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THE DIFFERENTIAL BLOOD CHANGES IN 
ASCARIASIS IN FOXES! 


By ARNoLD H. KENNEDY? AND RoNALD G. Law® 


Abstract 


Definite and characteristic changes occurred in the proportions of basophiles, 
neutrophiles and were cells in the blood of foxes which had received 
Ascaris eggs. The numbers of red blood cells and the amounts of haemoglobin 
remained within normal limits but the proportion of basophiles rose as high as 
69.5% of the total leucocyte count. This increase is apparently associated with 
the number of eggs given; the largest dose producing the highest proportions. 
The proportions of neutrophiles and leucocytes tended to fluctuate in opposite 
directions. An increase in the total numbers of white blood cells appears to be 
influenced by an increase in the proportion of basophiles. 


The cellular reaction of the blood to parasite infections indicates a definite 
attempt on the part of the body tissues to create a defence against them. 
The experiments described below were carried out in order to record the reac- 
tion of the various blood cells to Ascaris infection in foxes. 


Experiments 


Eight foxes, 10 to 12 weeks of age, giving a negative examination for one 
month by the D.C.F. technique, were infected with embryonated eggs from 
Toxascaris leonina. These and two control foxes were kept under identical 
conditions. The embryonated eggs were given orally in a small quantity 
of water by means of the stomach tube. The eggs were given to individual 
foxes in numbers approximating 500, 1,000 and 2,000. In one case, two 
doses of approximately 500 eggs were administered ten days apart, in order 
to record the cellular reactions when two small doses of eggs were given at a 
ten-day interval. After the eggs were administered the foxes were placed 
in individual wire cages so constructed as to allow the feces to pass through 
the wire to a drawer inserted some ten inches below. Previous to the adminis- 
tration of the eggs the blood of each fox was examined for normality on three 
occasions at intervals of three days. The standard for normal blood in foxes 
used in the investigation was that established by Kennedy (1). 

1 Original magus received July 4, 1934. 
Contribution from the Ontario Government Experimental Fur Farm, Kirkfield, Ontario. 


8 — Pathologist. 
. © Director. 
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Staining Reactions of Foxes’ Blood. 

According to Kennedy (1) the eosinophile is not found in foxes’ blood, but 
is replaced by a cell which is distinctly basophilic in staining reaction and 
which is classified as a basophile. Consequently, where investigations in- 
volving a rise in numbers of such cells are undertaken with foxes, eosinophiles 
are replaced by basophiles. 

The blood studies show the total number of red blood cells and leucocytes, 
the differential leucocyte count and the number of grams of haemoglobin 
per 100 cc. of blood. These counts were made on each fox at three-day 
intervals. The cell counts were estimated with a standard haemocytometer 
and the haemoglobin readings by a Sahli haemoglobinometer. The blood 
smears for the differential blood count were stained with Hasting’s stain, and 
two hundred leucocytes differentiated under the oil immersion lens. In all 
cases the foxes were bled in the morning, previous to feeding. Standardized 
pipettes were used throughout the study. . 


Red Blood Cells and Haemoglobin 


As stated above, the ten foxes, previous to infection, were all within the 
normal range of haemoglobin and of red blood cells per cu.mm. At no time 
during the experiment was there any trend towards abnormality in either 
red blood cells or haemoglobin. During the experiment the red blood count 
ranged from seven to twelve millions per cu. mm. of blood, the numbers 
increasing toward the conclusion of the experiment while the haemoglobin 
ranged from 11.5 to 12.9 gm. per 100 cc. of blood. These figures are within 
the limits established for the blood of normal foxes. 


Total Number of White Blood Cells 


In all of the foxes given embryonated eggs definite and characteristic changes 
in the proportions of basophile, neutrophile and lymphocyte cells occurred. 
The fluctuations of the basophiles and of the total numbers of white blood 
cells per cu. mm. appear very similar. The total number of white blood cells 
tends to increase or decrease according to the proportions of basophiles present 
in the blood stream. When the proportion of basophiles was within normal 
limits the total numbers of white blood cells were found to be normal also. 
Apparently the proportions of the various white blood cells are affected by 
Ascaris infection, the most striking change being shown in the proportion 
of basophiles. 

In Figs. 1-10 data from eight infected foxes and two controls are represented. 
In view of the similarity found in the counts of foxes receiving approximately 
the same number of Ascaris eggs, the graphs show the means for foxes given 
similar doses. 

Each graph was prepared from the data obtained on the total numbers of 
red and white cells per cu. mm. of blood. The proportions of the basophile, 
neutrophile and lymphocyte cells were also recorded at various dates through- 
out the experiment. The proportions of monocytes were so low that they 
were not plotted. The abscissas show the dates of the examination, while 
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the ordinates show the proportions or, in some cases, the numbers of basophiles, 
neutrophiles and lymphocytes found in the blood. The relations of the pro- 
portions of these three classes of leucocytes, as seen in foxes that have been 
given varying quantities of embryonated Ascaris eggs, can readily be studied 
from these graphs. 

Control Foxes G-31 and G-20 (Figs. 1 and 2) 

The proportion of basophiles fluctuated from 2 to 7%, which is within the 
limits for normal foxes’ blood. The general trend of the neutrophiles and 
lymphocytes was very similar to that of the basophiles, both being considered 
within the normal range. 
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Foxes G-21 and G-23 (Figs. 3 and 4) 

These foxes received approximately 500 embryonated Ascaris eggs. The 
proportion of basophiles reached the high point of 26% at 11 and 22 days 
after infection. After reaching the peak a decrease took place, the low point 
of 6.5% being reached on the thirty-ninth day. On the forty-eighth day an 
increase to 12.5% took place. An increase to 13.5% took place on the fifty- 
eighth day, after which the trend was downwards, the proportion of baso- 
philes ultimately returning to normal. 
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Foxes G-29 and G-30 (Figs. 5 and 6) 

Approximately 1,000 eggs were administered, and four days after infection 
the proportion of basophiles and the total numbers of white blood cells com- 
menced to increase. Fourteen days after infection the proportion of baso- 
philes reached a peak of 59.0% and the total number of white blood cells 


PER CENT 


rH Pe er 


ai ANE rr 
HIER eS 
GUL 





1 4 7111417 2123 = 29 32 3538 42 46 50 54 59 66 78 84 
pays 
Fic. 5. 


91 98 105 112 ung 


WHITE BLOOD 
CELLS /cu. 
IN THOUSANDS 





mom TTT Tit 
ee LULU Leet rl 
aE elt 
8 






oY NAR ile MAE ad 
HCE ee 
1470417 223 2REB DH H 66 78 8& a 9 105 ue ny 














282 CANADIAN JOURNAL OF RESEARCH 


reached a peak of 28,000 per cu. mm. of blood. On the thirty-second day the 
proportion of basophiles dropped to 22.5% and the total number of white 
blood cells to 8,900 per cu. mm. of blood. After this date the trends were 
more or less level, the basophile proportions varying from 19 to 31% and the 
total number of white blood cells from 8,600 to 13,700, until the ninety-first 
day, when a decrease occurred. On the one hundred and fifth day the pro- 
portion of basophiles had diminished to11%. After infestation the proportions 
of neutrophiles and lymphocytes decreased, the proportion of neutrophiles 
dropping from 44% to 18.5%, and the proportion of lymphocytes from 47 to 
22%, fourteen days after infestation. After this date the general trend of both 
neutrophiles and lymphocytes was upwards but their daily fluctuations were 
in opposite directions. In these foxes, the total number of white blood cells and 
the proportions of lymphocytes and basophiles had similar trends. With an 
increase in the total number of white blood cells, the proportion of basophiles 
also tends to increase. With an increase in the proportion of basophiles, the 
proportion of lymphocytes also increases over the proportion of neutrophiles 
and when the proportion of basophiles decreases the proportion of neutrophiles 
increases over the proportion of lymphocytes. 


Foxes G-27 and G-28 (Figs. 7 and 8) 


Approximately 2,000 embryonated eggs were given to these foxes and on 
the twenty-first day after infestation the proportion of basophiles reached 
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69%, decreasing to 15.5% forty-four days later. On the seventieth day 
there was an increase to 23%, and on the eighty-fourth day the level of 7.0% 
was reached. As the proportion of basophiles reached the highest peak the 
proportions of neutrophiles and lymphocytes decreased. After the basophilic 
low level was reached, the proportion of lymphocytes followed a trend similar 
to that of the proportion of basophiles and maintained higher levels through- 
out the experiment than did the proportion of neutrophiles. As the propor- 
tion of basophiles decreased, the proportion of neutrophiles increased, and 
their trends are almost opposite to each other. The trend of the numbers of 
white blood cells per cu. mm. followed very closely the trend of the basophiles. 
A peak of 37,000 was reached eighteen days after infestation. There was a 
decrease to 10,000 forty-four days after infestation. The decreases for both 
the total white blood cells and the proportion of basophiles occurred the same 
number of days after infestation. 
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Foxes G-37 and G-38 (Figs. 9 and 10) 

These foxes were given a preliminary dose of approximately 500 embryon- 
ated eggs and ten days later were given a similar dose. After the first dose a 
slight increase from 8.5 to 13% in the proportion of basophiles occurred, and 
in eleven days, when the second dose was given, the proportion was 3.5%. 
The proportion of basophiles then increased until twelve days after the second 
dose when it reached 49.5%. The proportion then decreased to 9.5% thirty 
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days after the second dose. The proportion increased to 18% at thirty-three 
days, and then showed a downward trend, and at eighty-one days after the 
second infestation reached 4.5% and remained within normal limits. 

As the proportion of basophiles reached the highest peak the proportions 
of neutrophiles and lymphocytes decreased, but these increased again when 
the proportion of basophiles began to decline. The proportion of lympho- 
cytes was greater than the proportion of neutrophiles thirty-three days after 
the second infestation and remained greater until the conclusion of the ex- 
periment. The fluctuations of the neutrophiles are opposite to those of the 
lymphocytes throughout. The total number of white blood cells per cu. mm. 
reached the peak of 23,000 twenty-three days after infestation, (which was 
the number of days in which the proportion of basophiles reached the high 
point), and declined to a low of 7,000 thirty-two days after infestation. Again 
the trends for both the total numbers of white blood cells and the proportion 
of basophiles were strikingly similar. 


Summary 


The proportion of basophiles rises to a marked degree after administration 
of embryonated Ascaris eggs to foxes. In some cases this proportion is as high 
as 69.5% of the total leucocyte count. The increase in the proportion of 
basophiles appears to be influenced by the number of eggs given; the largest 
dose producing the highest, and the smallest dose the lowest, proportion. 

There was a marked relation between the proportions of lymphocytes and 
neutrophiles. These tended to fluctuate in opposite directions. When the 
proportion of neutrophiles tended to rise the proportion of lymphocytes 
tended to decline. 

The red blood count and haemoglobin remained within normal limits 
throughout the experiment. 

The high proportion of basophiles tends to coincide with the marked 
increase in the total number of leucocytes. It demonstrates the fact that 
Ascaris infection can cause a marked rise to as high as 69% in the proportion 
of basophiles per cu. mm. of blood in the fox. The proportions of neutro- 
philes and lymphocytes are not involved to the same extent. Their actual 
numbers per cu. mm. of blood may be somewhat increased when the total 
numbers of white blood cells and proportion of basophiles reach a high peak 
but their proportions in the blood tend to be lowered. The increase in the 
total number of white blood cells over the normal appears to be greatly 
influenced by the increase in the proportion of basophiles; the lymphocytes 
and neutrophiles influence the increase to a minor extent only. 
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THE MICROBIOLOGY OF BUTTER 


I. THE YEAST AND MOLD COUNT OF BUTTER AS A MEASURE OF 
CREAMERY SANITATION! 


By F. W. Woop? ann H. R. THORNTON® 


Abstract 


Evidence is presented to show that the yeast and mold count of butter is 
inadequate as a sole measure of the microbial content of the butter and of 
creamery sanitation. 

It is recommended that for these purposes the bacterial count be used to 
supplement the yeast and mold count. 


Introduction 


The plate counts of yeasts and molds in butter are used extensively through- 
out Canada. Such counts may have a threefold purpose, viz., 

i. They may be used to measure the mold content of butter with a view 
to lowering the incidence of moldy butter. Consideration will be given to 
this question in a further communication from this department. 

ii. They may be used to increase the keeping quality of butter by maintain- 
ing its yeast and mold content at a minimum. A number of authors have 
reported that there is little or no relation between the yeast and mold counts 
and the keeping quality of butter. 

iii. They may be used as a measure of creamery sanitation. The present 
paper will be confined to a discussion of yeast and mold counts in this role. 


Literature 


Macy and Richie (1) review the literature on this subject and present a 
comprehensive bibliography. Reiteration of such readily available informa- 
tion will not be made here. This review serves to emphasize the need for 
more specific knowledge in this field. 


Procedure 


The butter samples were removed with sterile spatulas from two corners 
of packed butter boxes, thus insuring the inclusion of a surface portion, to 
sterile 4 oz. screw-top jars. After arrival at the laboratory the samples were 
melted in a water bath at 42-45° C. and were plated as soon as the butter 
reached a consistency permitting the use of a pipette. In almost all cases 
not more than 20 min. elapsed between the melting of the butter end the 
pouring of the plates with the medium. The usual 1 cc. pipettes were 
used and these were warmed to facilitate complete removal of the butter. 
The melted butter and the dilutions were shaken 25 times through a distance 


1 Manuscript received December 4, 1934. 
Contribution from the Department of Dairying, University of Alberta, Edmonton, Canada. 
The data contained in this paper are taken from a thesis presented by Mr. Wood to the Committee 
on Graduate Studies, University of Alberta, in partial fulfilment of the requirements for the degree 
of Master of Science. 
2 Graduate Student. 
* Professor of Dairying, University of Alberta. 
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of approximately 1 ft. The diluent was either physiological salt solution or 
tap water and losses due to sterilization were made up by the aseptic addition 
of sterile solution. 

The medium used for the plate counts of bacteria was standard nutrient 
agar. That for the yeast and mold counts was Bacto Dehydrated Malt 
Agar acidified to pH 3.5+0.1, since this acidity is rapidly becoming standard. 
pH measurements were made by the electrometric method* using a saturated 
calomel-potassium chloride half-cell and a quinhydrone electrode in solidified 
medium at room temperature. Just prior to plating, acidification with a 
1:10 aqueous dilution of 85% U.S.P. lactic acid was accomplished. The 
amount of acid specified by the manufacturers was not in all cases sufficient 
to lower the pH to 3.5. 

Before pouring the medium, the petri dishes containing the inoculums were 
placed upon a warm plate to melt the butter. This allowed more thorough 
mixing with the medium. The inoculums did not reach a temperature above 
45° C. 

The yeast and mold counts were made after incubation of the plates at 
25° C. for five days. This procedure was satisfactory except when a large 
number of molds were present, in which case it was found preferable to count 
at the end of three days and before the mold colonies lost their identity 
through merging. The bacterial plates were incubated at 25° C. for four days. 

The yeast and mold colonies were counted in transmitted light with the 
naked eye. The lumi-lens plate counter used for counting bacterial colonies 
proved unsatisfactory for yeast and mold work. When necessary a hand lens 
was used and, in the interests of accuracy, a hand tally was employed through- 
out. All counts are plate counts and are reported on a per cc. basis. 

Mycoderms in. colony form resemble molds of the Oidium type, while in 
stained preparations the individual cells are more nearly like yeasts. In the 
absence of any standard for the classification of these micro-organisms for 
counting purposes they are included with the yeasts in this investigation, 
except as noted. 





























Experimental 





Air Contamination during Plating 

Contamination by molds from the air seems to be frequent enough to 
have won for these micro-organisms the pseudonym, ‘‘Weeds of the Labor- 
atory’’. A surprisingly large proportion of modern butters are so low in 
mold content that air contamination of butter surfaces and plates during 
analytical procedures might be expected to introduce a percentage error of 
considerable magnitude. 

In order to obtain information on the extent of this error 93 control plates 
were interspersed with the routine plates on 65 days. The majority of these 
plates remained sterile. The largest number of colonies appearing on any 
one plate was two and this number occurred on a total of four plates. 
















*The son apparatus used in this work was purchased with funds granted by the 
National Research Council of Canada. 
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In another study of this problem 25 open plates containing solidified medium 
were exposed on the desk where plating was proceeding, on 12 days. The 
period of exposure varied from 20 to 30 min. Seven mold colonies developed 
which were confined to six plates and the highest number of colonies on any 
one plate was two. No yeast colonies were recognized in either study. 

Because of the location and construction of the laboratory in which this 
work was conducted, a high incidence of air contamination was expected. 
These data indicate that, contrary to expectation, the error introduced into 
yeast and mold counts of butter by air contamination during plating is likely 
to be so small as to be negligible. 


The Growth of Micro-organisms during Plating 


It is frequently convenient during routine butter analysis, and necessary 
during research, to hold the butter in a melted state for considerable time, 
although this practice is frowned upon in The Preliminary Butter Report (2). 
Five samples of butter were held at 43° C. for periods up to four hours, plates 
being poured from these every hour. In no case was any significant change 
observed in the counts of yeasts, molds or bacteria during this period. 


The Relations of the Numbers of Yeasts, Molds and Bacteria in Butter. 


The yeast, mold and bacterial counts of 46 samples of exhibition butter 
from 13 Alberta creameries are presented in Table I. It is apparent that 
there is no relation between the counts of these three types of organisms 


TABLE I 


THE PLATE COUNTS OF BACTERIA, MOLDS AND YEASTS IN 46 SAMPLES OF EXHIBITION 
BUTTER FROM 13 ALBERTA CREAMERIES 

































1 250 0 2 7,800 0 4 
2 400 0 2 8,700 7 0 
3 400 1 4 8,900 1 1 
4 400 6 7 10,600 11 8 
5 1,400 0 1 12,000 3 5 
6 1,400 0 0 15,800 1 47 
7 1,600 0 0 18,800 18 9 
8 1,800 0 30 19,700 1 93 
9 2,000 0 12 21,700 1 8 
10 2,100 28 0 22,100 1 8 
11 2,100 28 0 22,800 5 17 
12 3,000 0 0 25 ,000 0 140 
13 3,500 0 0 29 ,000 0 2,140 
14 3,800 16 21 31,800 2 5 
15 4,100 8 0 ¢ 38,000 1 19 
16 4,600 1 7 44,500 10 1,680 
17 4,600 1 157 55,000 1 1,900 
18 4,600 0 5 68 ,600 0 0 
19 5,600 1 8 75 ,400 4 26 
20 5,700 2 35 130,000 0 12,000 
21 6,100 0 0 136,000 34 28 
22 6,200 12 12 400 ,000 0 480 
23 7,200 0 16 ver 1,000,000 0 0 
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except that most of the high yeast counts occur in butter having high bac- 
terial counts, but it does not follow that a low yeast count is invariably accom- 
panied by a low bacterial count. The same type of relation is seen to exist 
in the 44 samples of fresh commercial butter from 16 Alberta creameries as 
listed in Table II. There can be no doubt that the yeast and mold counts 
constitute a very incomplete measure of the microbial content of many of 
these butters. 
TABLE II 


THE PLATE COUNTS OF BACTERIA, MOLDS AND YEASTS IN 44 SAMPLES OF FRESH 
COMMERCIAL BUTTER FROM 16 ALBERTA CREAMERIES 





—— Bacteria Molds | Yeasts Bacteria Molds | Yeasts 





a a a | 










1 650 1 12 23 27,250 0 
2 900 0 25 24 33,500 50 360 
3 2,000 0 45 25 42,000 4 4 
+ 2,400 0 0 26 44,300 35 
5 3,600 3 0 27 46,500 2 230 
6 5,600 0 0 28 48 ,500 4 
7 6,250 0 57 29 50,000 55 250 
8 6,500 1 3 30 55,000 30 180 
9 9,000 0 4 31 61,000 0 150 
10 9,600 1 0 32 79,000 1 131 
11 10,600 0 111 33 80,000 4 4 
12 10,900 2 2 34 93,500 0 980 
13 12,900 1 1 35 95,500 0 270 
14 15,000 45 485 36 106 , 500 2 990 
15 15,500 60 1,120 37 107 ,500 15 4,500 
16 15,500 2 3 38 113,000 80 310 
17 15,800 0 7 39 145,000 | 250 2,350 
18 17,000 40 350 40 196 ,000 1 22 
19 17,700 0 3 41 600 , 000 1 210 
20 18,200 0 41 42 1,000,000 0 0 
1,000,000 0 
1,000,000 4 


There appear to be at least three possible explanations for the general lack 
of correlation between numbers of yeasts, molds and bacteria in butter, viz:— 
i. Their comparative resistance to pasteurization temperatures. 
ii. Their comparative growth rates. 
iii. Their comparative resistance to the methods employed in the steriliza- 
tion of creamery equipment. 


A thorough study of the flora of cream pasteurized at high temperatures, 
as practised in Canadian creameries, has not yet been made and will be 
necessary before it can be said with certainty that pasteurization has no 
bearing upon the subject under discussion. However, surveys made on 
Alberta creameries by Dr. Linneboe, Dairy Analyst, Provincial Department 
of Agriculture, and by the writers, have not shown high yeast and mold 
counts or particularly high bacterial counts on the pasteurized cream even 
when high-count butter was encountered. In almost all cases the evidence 
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has pointed to insanitary equipment. This is well illustrated in Table III 
which shows the counts of five samples of cream before and after com- 
mercial pasteurization. No yeasts or molds were shown to survive the heat 
treatment and the highest bacterial count after pasteurization was 33,000. 


TABLE III 


THE PLATE COUNTS OF THE BACTERIA, MOLDS AND YEASTS IN CREAM BEFORE 
AND AFTER COMMERCIAL PASTEURIZATION 


Cc ao Raw Pasteurized 
oO. fn 
Bacteria Molds Yeasts Bacteria Molds Yeasts 
1 46,000,000 300 7,000 2,500 0 0 
2 12,000,000 10 6,600 600 0 0 
3 500 , 000 ,000 300 15,600 33,000 0 0 
4 49,000,000 10 640 1,900 0 G 
5 — a -- 16,800 0 0 


In a favorable environment the generation times of bacteria are much 
shorter than those of yeasts and molds. Bacteria have been known to 
reproduce as quickly as once in 17 min., but this generation rate in a mixed 
flora in butter seems unlikely. Yeasts reproduce much more slowly than 
bacteria, requiring, so it is reported, at least three hours in favorable cir- 
cumstances. The life cycle of the molds is much longer still. Because of 
the common practice of holding cream for some hours between pasteurization 
and churning and the probability of the holding temperature being high enough 
to allow growth in some cases, it is not difficult to visualize the growth differ- 
ences which may be effected in such cream. The tendency would be for the 
bacteria to reproduce much more rapidly than the yeasts or the molds, and 
the result would be a butter high in bacteria and with few or no yeasts or 
molds. The same condition may be produced when the butter itself is . 
exposed to high storage temperatures. 

To study the comparative resistances of bacteria, yeasts and molds to churn 
“sterilization”, a churn was filled approximately one-third full of water at 
80-90° C. and revolved from 20-30 min. It was then drained and filled with 

a definite quantity of cold water of 


TABLE IV very low bacterial content. The 


THE PLATE COUNTS OF THE BACTERIA, MOLDS plate counts of 1 cc. of such rinse 
AND YEASTS IN RINSE WATER FROM 


a, “uuemaane” ennenes water in five trials on widely sepa- 


rated days are reported in Table 
IV. It4s to be expected that butter 





Trial Bacteria Molds | Yeasts : ‘ 
|_| Sséamaadée in this churn from low count 
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YEAST AND MOLD COUNT OF BUTTER 


The Effect of External Butter Workers 


Some interesting observations were made upon butter from churns equipped 
with the external type of butter worker. An opportunity presented itself 
for the study of butter manufactured in such churns in seven creameries. 


TABLE V 


THE PLATE COUNTS OF THE BACTERIA, MYCODERMS AND YEASTS IN 83 SAMPLES OF BUTTER 
FROM CREAMERY A, EQUIPPED WITH A CHURN HAVING AN EXTERNAL BUTTER WORKER 
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This butter was characterized by a rather regular, high, yeast and mycoderm 
content. Table V shows the bacterial, yeast and mycoderm counts of 83 
samples of butter from a representative churn of this type, the counts being 
arranged in ascending order of the bacterial counts. It is seen that there is 
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no tendency for the mycoderm or yeast counts to increase with an increase in 
the bacterial counts. Similar counts of butter from consecutive churnings 
from two other creameries equipped with external butter workers are given 
in Table VI. Churning number 5 of Creamery C immediately followed the 
dismantling and thorough cleaning of the butter worker. 


The yeast and bacterial counts of 70 samples of butter manufactured in 
11 creameries known to be equipped only with churns having internal workers 
are set forth in Table VII. The mycoderm count of each sample was either 
zero or was so low that it was not considered separately from the yeast count. 
If the churns which we studied are representative then it seems probable that 
the external type of butter worker may contribute large numbers of yeasts 
and mycoderms to butter. One such worker was dissembled and decom- 
posing material taken from one of the joints gave the following analysis 
per gm.: bacteria 1,000,000,000; yeasts 10,000,000; mycoderms 1,300,000. This 
material does not show the same relation between the counts of the dif- 
ferent kinds of micro-organisms as is shown by the manufactured butter. 
Two possible explanations are: 


1. The inadequacies of the plate method of analysis. 


2. An equilibrium, differing from that above, which may exist at the very 
surface or edge of the material im sttu from which the contamination of the 
butter takes place. 


Discussion and Recommendation 


It is common—almost to the point of being universal—practice in Canadian 
creamery butter manufacture to pasteurize cream at not less than 170° F. 
for 10 min. and starters are not used for cream ripening after pasteurization. 
The bacterial content of butter so made, therefore, has a different significance 
than when low-temperature pasteurization and starters are used. In these 
circumstances a very high bacterial content in the finished butter denotes 
one, or a combination, of the following faults:— 

i. Large numbers of thermoduric bacteria in the raw cream. This seems 
unlikely to happen but is included here in the absence of more extensive 
proof to the contrary than is offered in this paper. In any case this would 
probably constitute a serious defect in butter. 


ii. Improper pasteurization. 
iii. Growth after pasteurization either in the cream or in the butter. 
iv. Contamination subsequent to pasteurization. 


It is believed that the evidence presented in this paper shows conclusively 
that the yeast and mold count of butter is an incomplete measure of its 
microbial content and is frequently misleading when so used. Therefore, it is 
recommended that, if the maximum reduction of the number of micro- 
organisms in Canadian butter is to be effected, the yeast and mold count be 
supplemented with the bacterial count. 
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TABLE VI 


THE PLATE COUNTS OF BACTERIA, MYCODERMS AND YEASTS IN CONSECUTIVE CHURNINGS FROM 
CREAMERIES B AND C, EACH EQUIPPED WITH CHURNS HAVING EXTERNAL BUTTER WORKERS 


Churning Creamery B Creamery C 
No. = $a | 
Bacteria Mycoderms | Yeasts Bacteria Mycoderms | Yeasts 
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TABLE VII 


THE PLATE COUNTS OF BACTERIA AND YEASTS IN BUTTER FROM 11 ALBERTA CREAMERIES 
EQUIPPED WITH CHURNS HAVING INTERNAL BUTTER WORKERS 
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Summary 


1. Air contamination of plates and butter surfaces with yeasts and molds 
during plating is likely to be small. 


2. No significant changes were observed in the bacterial, yeast and mold 
counts of five samples of butter held in the melted state for four hours. 
3. The only relation observed between the different kinds of micro- 


organisms in Lutter was that high yeast counts tended to occur in butter of 
high bacterial content. 


4. The possibility of the differing growth rates of bacteria, yeasts and molds 
and their resistance to equipment ‘“‘sterilizing’”’ procedures being responsible 
for the general lack of relation between the counts of these organisms in 
butter is discussed. 

5. Butter manufactured in seven churns equipped with external butter 
workers was characterized by a high yeast and mycoderm count. 


6. In view of the manufacturing practices common in creameries in Canada 
a very high bacterial content in Canadian creamery butter is evidence of a 
fault in the production or manufacture of the butter. 


7. The yeast and mold count of butter is an incomplete measure of its 
microbial content and consequently of creamery sanitation. 

8. It is recommended that consideration be given to supplementing the 
yeast and mold count of Canadian creamery butter with the bacterial count. 


References 


1. Macy, H. and Ricuie, H. B. The mold and yeast count as an index of the keeping quality 
of butter. J. Dairy Sci. 12 : 351-366. 1929. 


2. THE PRELIMINARY BuTTER REPORT. Committee on Bacteriological Methods of Analysing 


Dairy Products; American Dairy Science Association. J. Dairy Sci. 13 : 394-405. 
1930. 





| 
i 
| 
i 
' 
| 


295 


THE MICROBIOLOGY OF BUTTER 
II. THE GROWTH OF MOLDS IN AND UPON BUTTER! 


By. H. R. THORNTON? AND F. W. Woop? 


Abstract 


The improbability of the initiation of mold growth in the incorporated moisture 
of salted and unsalted butter and in the free moisture of salted butter is dis- 
cussed. An attempt is made to relate the initiation of mold growth to con- 
densation moisture on the various surfaces of butter. Experiments showing the 
comparative effect on the initiation of mold growth of four methods of parch- 
ment treatment are described. Cladosporium sp. was isolated from a butter 
showing the interior type of molding. 


Introduction 


Although the creamery butter of Alberta has been subjected to a yeast 
and mold control program of ten years’ duration, moldy butter continues 
to be a source of economic loss. Its sporadic and puzzling recurrence on 
prints and solids suggests that the control measures have failed to reach the 
vital spot of the trouble. Dependence has been placed almost solely on the 
lowering of the mold content of the butter although this is but one of many 
factors which may influence the growth of molds. 

Our knowledge of the physical and chemical structure of butter is yet too 
incomplete to permit prediction, or recognition with certainty, of all the 
conditions affecting the growth of molds in or upon butter. It is possible, 
even probable, that adequate control of the molding of butter will be deferred 
until more information relating to the chemistry of butter is available. In 
the meantime there is nothing to be lost in maintaining a viewpoint which 
takes cognizance of the few available facts regarding the chemical nature of 
butter on the one hand and the growth requirements of molds on the other. 
The present discussion is an attempt to correlate some hitherto apparently 
contradictory observations and to make some small contribution to our 
knowledge of the subject. 

Procedure 


The general plating procedure described in the first paper of this series 
was followed in the present investigation. 

One-pound prints of commercial creamery butter low in mold content were 
each divided into two approximately equal parts. One part was inoculated 
with mold spores while the other was left uninoculated. Each part was 
then subdivided into portions of approximately one-tenth pound each. These 
portions were wrapped as uniformly as possible with parchment treated as 
described later. Unless otherwise stated, reference will be made only to 
salted butter. . 

1 Manuscript received December 4, 1934. 
Contribution from the Department of Dairying, University of Alberta, Edmonton, Canada. 
The data contained in this paper are taken from a thesis presented by Mr. Wood to the Committee 


on Graduate Studies, University of Alberta, in partial fulfilment of the requirements for the degree 
of Master of Science. 


2 Professor of Dairying, University of Alberta. 
3 Graduate Student. 
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Two cultures of Penicillium, isolated from moldy butter and parchment, 
were the source of the inoculum. Slants of these strains, sufficiently aged 
to contain an abundance of spores, were employed and the inoculum was 
mixed as thoroughly as possible with the butter to effect even distribution 
of the spores. Owing to the fact that the action of specific molds was not 
being studied no precautions were taken to maintain aseptic conditions during 
the inoculation and wrapping of the butter. It is probable that many of 
these samples became inoculated with extraneous molds during the manipula- 
tions of preparation and observation. Therefore, while these results can be 
interpreted with certainty only for the two cultures of Penicillium used, it is 
possible that they are representative for many species. After the butter 
was wrapped mold spores were sprinkled on the outer surfaces of the parch- 
ment in the case of butters 1, 4 and 5 (Table [). 


The four methods of treating the parchment prior to wrapping were as 
follows :— 


i. Sterile parchment—dampened and autoclaved for 1 hr. at 15 Ib. pressure. 


ii. Hot brine —soaked for 1 hr. at 97° C. in brine which was saturated 
with sodium chloride at room temperature. 

iii. Cold brine —soaked in saturated sodium chloride solution for 
several hours at room temperature. 

iv. Dry parchment —no treatment. 


Incubating the samples together in one large humidifier proved unsatis- 
factory because of the collection of condensed moisture and the transference 
of this moisture with salt from one sample to another. Each sample was, 
therefore, placed in a sterile screw-top ointment jar in the bottom of which 
was 2 cc. of sterile water. These jars were then incubated at 12-15° C. in a 
chamber with a moisture saturated atmosphere. 


Types of Molding 


Two types of mold growth on butter are frequently confused both in 
the industry and in the literature. In the first type, mold growth occurs 
on the surface of the butter and—another instance of the same thing— 
in the cracks, crevices and creases of poorly packed butter. This type will 
be termed surface molding in the present discussion. It occurs very frequently 
on commercial butter but its presence does not always mean large monetary 
loss since the outer layers of butter can be removed leaving the freshly cut 
surface free from mold growth. 

The second type of molding occurs as smudgy, dark areas or patches 
throughout the interior of the butter mass. This type is encountered com- 
paratively rarely but its occurrence usually results in almost complete loss 
of the butter since it cannot be readily removed and the butter is unmarket- 
able. In the present investigation it will be referred to as interior molding, 
and discussion relating to it will be confined to the section entitled ‘Interior 
Molding”’. 
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Salt Relations 


Salt Concentration in the Butter 

After making an extensive review of the literature of this subject Macy (4) 
concludes that ‘‘the influence of salt varies depending upon the species and 
the substrate. Altho unsalted butter is most commonly affected, mold often 
appears on salted butter. This occurrence demands explanation.” 

In the present investigation the ease with which mold growth can be induced 
on the surface of unsalted butter was noted. At the same time no sample 
of salted butter was encountered which did not become moldy when stored 
under favorable conditions. The maximum salt concentration of the butters 
used in these experiments was 2.5%. 

While there is a degree of uniformity in the results of different investigators, 
the divergencies which appear are explainable, in part at least, by known 
factors. It is not to be expected that mold growth is initiated in the fat 
itself. Mold growth will, therefore, be influenced not by the salt concentration 
of the butter but by the salt concentration of the butter moisture. The salt 
concentration of the butter moisture is governed not alone by the percentage 
of salt and moisture in the butter but by the distribution of the salt and mois- 
ture as well. Microscopic examination of butter usually reveals undissolved salt 
crystals—proof that the salt is not homogeneously distributed. The travel 
of moisture in salted butter is evidence of a tendency for the salt to become 
homogeneously distributed throughout the moisture. But, owing to the 
slow diffusion of moisture through butter, an entirely homogeneous dis- 
tribution of the salt is very improbable. It is apparent that the percentage 
of salt in the butter is not a satisfactory measure of the inhibitory effect of 
salt even when consideration is given to the salt tolerances of different molds. 
But it is highly probable that the salt concentration of the butter moisture 
of commercial salted butters is sufficiently high to act as an important in- 
hibitory agent in view of the other unfavorable conditions of butter storage. 
This probability receives support in the observation that unsalted butter 
becomes moldy with ease and in the fact that the growth of mold on different 
salted butters of the same salt and moisture concentration is not more erratic 
than the growth of mold upon the butter in different boxes of the same 
churning. The fact that all the butters observed in this work became moldy 
when conditions other than salt concentration were optimum points to some- 
thing other than salt as the controlling factor in the case of salted butters. 


Salt on the Butter Wrappers 

Salt has proved to be an effective agent in the prevention of mold growth 
on the parchment liners and wrappers used in the packaging of butter. A 
number of investigators recommend boiling the parchment in saturated salt 
solution. Macy and Pulkrabek (5) present an excellent review of the liter- 
ature on this subject. 

In recent years some creameries have resorted to the use of dry parch- 
ment as received from the manufacturer for the wrapping of print butter. 
It has been reported to us that the adoption of this practice has been attended 
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with some success in reducing the incidence of moldy butter. It is claimed 
that this method of controlling mold growth on butter has been used with 
success in New Zealand. 


Data are presented in Table I to show the effect of parchment treatment 
upon subsequent mold growth both on the butter and on the parchment on 
incubation for periods varying up to 21 weeks in a moisture-saturated at- 
mosphere. It is seen that a wet sterile wrapper gives the least protection 
against moldiness, while dry wrappers were fairly effective for three weeks 
but not for eight weeks. Wrappers soaked in saturated salt solution, either 
hot or cold, gave the maximum protection. 


Oxygen Relations 


In reviewing the effect of oxygen upon mold growth on butter, Macy (4) 
concludes that ‘““The notion that molds are strict aerobes seems to be faulty, 
but it is evident that they do require some oxygen for their norfnal develop- 
ment, altho excessive quantities are inhibitory.”” He was able to obtain 
meagre growth on butter in partial pressures of oxygen and no growth in 
the absence of oxygen. In the present investigation initiation of growth 
of any mold studied was not observed on acidified malt agar in the complete 
absence of oxygen in Spray anaerobic culture dishes. The fact that mold 
growth is observed infrequently in butter at any place other than the surface 
seems to justify the conclusion that mold growth in or upon butter implies 
the presence of considerable quantities of oxygen. Therefore, any theory 
in explanation of the molding of butter will have to presuppose the presence 
of oxygen. 


The modern conception of organic complexes in solution is that they respire, 
consuming oxygen. The moisture of butter contains an organic complex 
in solution. Therefore—by analogy—it consumes oxygen. When this 
moisture is in situ in the butter the diffusion of oxygen must necessarily be 
very slow. Because of this it is highly probable that an anaerobic state is 
rapidly reached in this moisture in salted and unsalted butter after manu- 
facture. The theory that mold growth is initiated in the incorporated butter 
moisture then becomes untenable. 

When the butter moisture “‘leaks’’ onto the surface it is probable that the 
diffusion of oxygen from the atmosphere is more rapid than the consumption 
of oxygen by the butter moisture, and aerobic conditions result. This 
assumption is justified by the behavior of other solutions of organic complexes. 
The salt in the moisture of salted butter then tends to inhibit the growth of 
molds in the expressed moisture. 


Shutt (7) claims that “Alternaria, Spicaria and Oospora types can exist in 
the butter and spread rapidly when the mycelium reaches the surface, where 
more oxygen is available. Perfect union between butter and wrapper seems 
to be no hindrance to the growth of these molds.’’ Not only do these two 
statements seem to conflict but the author presents no data in support of his 
contentions and does not describe the technique by which he was able to make 
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such observations. Many workers do not differentiate between the growth 
of molds and the initiation of growth. The latter is the important considera- 
tion in the question under discussion. No theory of the structure of butter 
yet propounded suggests that the conditions immediately below the surface 
of butter differ from those in the farther interior. If the quoted statements 
are true, we are at a complete loss to explain the infrequent occurrence of the 
interior type of molding. 

Some authors have referred to molds growing through parchment. During 
the course of the present study the piercing of parchment by molds has not 
been observed and proof of this ability was not forthcoming with any of a 
number of techniques devised for the purpose. 


Moisture Relations 


The great lack of uniformity in the occurrence of moldiness on different 
samples of butter of the same brine concentration and even on different 
boxes of the same churning indicates that the key to control lies not in the 
butter moisture or salt concentration. Without exception, investigators 
who discuss this subject report the necessity for storage of the butter in 
atmospheres of high relative humidity in order to induce macroscopic evidence 
of mold growth (4). Shutt (7) was able to germinate and grow on the outer 
surfaces of butter wrappers the varieties of molds with which he worked 
when the relative humidity of the storage rooms was 100%, but was unable 
to do so when the relative humidity dropped to 90%. Tomkins (9), while 
studying the relation of mold growth on different surfaces to the relative 
humidity of the atmosphere, found that the spores of Allernaria citris germin- 
ated and the mycelium grew to a length of 7.7 yw in an atmosphere with a 
relative humidity of 100% and only 0.3 » when the relative humidity was 
90.8%. Skovholt and Bailey (8) failed to obtain macroscopic evidence of 
mold growth on bread crumbs held in atmospheres of less than 90% relative 
humidity. 

In the present investigation no sample of butter became moldy on incuba- 
tion for 60 days at ice-box or room temperature and humidity. Every 
sample quickly became moldy when removed to a moist chamber at the same 
temperatures. It is apparent that, apart from temperature, the humidity of 
the storage atmosphere is the greatest single factor influencing the molding 
of butter, and the evidence to date does not favor the theory that mold 
growth is initiated in the butter moisture, either free or incorporated. 

Atmospheric humidities are related to the size of moisture condensation 
films on surfaces coming into contact with the atmosphere. The condensation 
moisture on the surface of butter is likely to exist in the form of droplets 
rather than as continuous films due to surface tension effects. The work of 
Collins and Hammer (1) on the migration of bacteria through butter lends 
support to this theory. An exception would occur where the parchment 
forms a capillary space with the surface of the butter, at which point the mois- 
ture would form asafilm. In our experimental work the collection of moisture 
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in the folds and creases of the parchment during storage at high humidities 
has at times been so excessive as to cause the moisture to run off the parch- 
ment, carrying with it sufficient salt to change the results of the experiments. 
Many samples of butter in commercial storage were observed to be moldy at 
those places where moisture would be expected to collect in large amounts if 
the condensation moisture had been excessive. Condensation moisture would 
contain no salt while it would quickly dissolve nutrient material from the 
butter and the parchment. Mold could, therefore, grow in this solution. 
If the condensation moistuse came into contact with the salted butter moisture, 
the salt concentration would be lowered and its inhibitory effect lessened. 
Support is given to this theory in the investigations of Vernon (10) who 
found that free moisture on the surface of butter was superior to pure water 
as a pabulum for molds. 


If the foregoing is representative of the conditions existing on the surface 
of butter, then control measures will have to include methods of preventing, 
as far as possible, the formation of condensation moisture. Observations 
both in the laboratory and in the industry are to the effect that the storage 
of butter in atmospheres of low humidity tends to decrease the precipitation 
of moisture on the butter surfaces. The more carefully the butter is packed 
the less will be the tendency to leave microscopic or capillary spaces between 
the butter and the parchment and unnecessary creases or folds in the 
parchment. 

As has been pointed out by previous authors, there is another source of 
salt-free moisture in the manufacture of butter. This is rinse water, water 
from wet equipment, etc. which falls on the finished surface of the butter or 
is carried there by wet parchment. Like the condensation moisture this 
water is likely to be at oxygen pressure equilibrium with the atmosphere. 
It will dissolve nutrient material from the butter and parchment—thus per- 
mitting mold growth—and will dilute the free butter brine which may be 
expressed onto the surface of the butter. The control measure which imme- 
diately suggests itself is to keep this type of moisture away from butter. 


Temperature Relations 


The influence of temperature on the growth of molds on butter is too well 
known to justify enlargement here. Suffice it to emphasize that the tem- 
peratures of butter storages in this province cannot always solely be depended 
upon to inhibit mold growth but that the incidence of this butter defect is 
lessened as temperatures are lowered. 

There is, however, another relation dependent upon temperature which 
deserves consideration. The air which is incorporated into.butter is the air 
which is in the churn during manufacture. This air is incorporated at 
churning temperature and is likely to have a relative humidity approaching 
100%. The air which may be entrapped contiguous to the surface during 
the packing and wrapping of the butter will be that of the churn room and 
will likely be high in temperature and humidity. When the butter is stored, 
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the temperature of this air will fall and moisture will be precipitated from it 
on the butter surfaces, even on those of the incorporated air pockets. When 
the butter is again removed to higher temperatures, more moisture will be 
condensed from the warmer atmosphere on the cold surface of the butter and 
particularly on the parchment. Vernon (10) reports an increased incidence 
of moldy New Zealand butter in London during warm weather and relates 
the phenomenon to a greater condensation of atmospheric moisture on the 
cold surfaces of butter when it is removed from storage. 


Contamination Relations 


A number of Alberta buttermakers have reported to us their experience 
that the molding of commercial butters seems to bear no relation to the mold 
count of the butter. A few such butters were investigated by us and the 
results substantiate the opinion of the practical men. Our attention has 
been called to the fact that the incidence of moldy butter has decreased with 
the decrease in the mold content of the butter in one of the provinces of 
Canada. This, however, does not constitute proof that the two phenomena 
are directly related in the absence of information regarding the other variables. 

Combs and Eckles (2) were unable to find any relation between the mold 
count and the subsequent molding of the butter. Shutt (7) on the other 
hand states that ‘“‘molds can develop on the butter under the wrapper when 
the initial mold spore contamination of the butter stored is very high.” 
He presents no data in support of this statement. 

In the present investigation no sample of butter was encountered which 
became moldy when the humidity of the storage atmosphere was unfavorable. 
No sample failed to become moldy when the humidity of the storage atmosphere 
was favorable. Scrutiny of Table I reveals that there is no relation between 
the mold count and the subsequent development of moldiness either on the 
inoculated or uninoculated butters there reported. Our observations, there- 
fore, confirm those of Combs and Eckles. 

These results appear reasonable when consideration is given to all the 
factors affecting the growth of molds upon butter. Apart from temperature, 
the humidity of the atmosphere surrounding the butter has more influence 
upon the development of moldiness than any other single factor and, probably, 
than all other factors combined. In the commercial handling of butter this 
is likely to vary more than any other condition, the variations being sufficient 
to account for the erratic appearance of molding in different boxes of the same 
churning. 

Many of the incorporated mold spores and mfcelial portions are likely to 
be imbedded in the fat of the butter and cannot, therefore, initiate growth. 
There is no evidence that incorporated spores can initiate growth of the 
surface type in salted butter. Common observation suggests that, if this 
occurs at all, it does so infrequently and cannot constitute the major cause of 
this type of molding. It is highly probable that few or no samples of com- 
mercial butter are finished without mold spores or mycelium dropping on their 








GROWTH OF MOLDS IN BUTTER 303 


surfaces. It is possible that these unincorporated mold spores and mycelial 
portions constitute the major source of trouble. 


The results here cited should not be interpreted as meaning that the authors 
consider the mold count useless in the control of moldy butter. They do not 
constitute proof that there is no relation between the mold content of butter 
and the subsequent development of mold. But they can be taken as evidence 
that, if any relation exists, the methods of measurement at present available 
are not sufficiently exact to detect it. Moreover, an indirect relation between 
the mold count and the’ incidence of moldiness probably does exist. Any 
campaign which has for its purpose the lowering of the mold count of butter 
will undoubtedly eliminate many sources of external contamination. 


Interior Molding 


The fact that the interior type of molding is of comparatively rare occur- 
rence while surface molding is easily induced upon commercial butter imme- 
diately focuses attention upon the peculiar chemical and physical structure 
of butter. It has already been pointed out that the theory that mold growth 
is initiated in the incorporated butter moisture is untenable. Butter is known 
to contain appreciable quantities of air. Pickerill and Guthrie (6) found 
from 0.5 to 14.9% of the butter volume to be air, with an average of 4%. 
Although the nature of the distribution of this air is not well known, it is 
apparent that a considerable proportion of the microscopic droplets in butter— 
usually called moisture droplets—are in reality air globules. King (3) points 
out the difficulty of recognizing these air globules with certainty in the micro- 
scopic examination of butter. 


The incorporated air in butter comes from the air in the churn. This air is 
high in humidity and on the temperature being lowered during storage con- 
densation moisture will be deposited in minute droplets at the air-oil interface. 
The microscopic picture of many of the droplets in butter is such as to suggest 
that this has taken place. The conditions at the air-oil interface within the 
butter become not entirely unlike those at the outer surface of the butter and 
an environment favorable for mold growth is provided. 

Aneopportunity presented itself for a microscopic examination of a butter 
displaying this type of mold growth. When material from the dark smudgy 
areas in the butter was viewed under the microscope it was observed that all 
growth originated in droplets and it was assumed that these were air droplets. 
In many of these droplets the exact point of growth initiation was recognizable 
and was invariably at the air-oil interface. The mycelium then grew in the 
droplet or globule, frequently filling this with a coil ot mycelial growth which 
was sometimes unable to pierce the butter. At other times the mycelial 
growth penetrated the butter for considerable distances, often through other 
droplets. These observations emphasize the difference between the point 
of initiation of growth and the direction of subsequent mycelial growth. 

A mold was isolated from the moldy patches and was identified by three 
independent workers as Cladosporium sp. Morphologically it was indis- 
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tinguishable from the mold observed in situ in the butter. The mold count of 
this butter was low and no species of Cladosporium was observed on the plates. 
An explanation was provided for this in that the mold which was isolated 
would not grow into macroscopic colonies in the five-day incubation period. 
It is interesting to note that other workers have reported Cladosporium sp. 
as the cause of this type of molding (10). 


Summary 


1. Two types of butter molding are at present recognized, viz.,— surface 
molding, which occurs frequently and interior molding, which occurs in- 
frequently. 


2. Emphasis is placed upon sharp differentiation between the point of 
initiation of mold growth and the direction of subsequent mycelial growth. 


3. The evidence to date is interpreted as indicating that mold growth is 
not initiated in the butter moisture, either incorporated or expressed, of 
salted butter because, 

(a) Salt is an important inhibitory agent, especially in the presence of other 
unfavorable conditions inherent in commercial butter storage. 

(6) Incorporated butter moisture is likely to become anaerobic within a 
short time after manufacture. 


4. No sample of butter was encountered which became moldy at ice-box 
temperature when the humidity of the surrounding atmosphere was low. 
No sample failed to become moldy at the same temperature when the humidity 
of the surrounding atmosphere was high. 


5. Condensation moisture is suggested as the most probable major factor, 
apart from temperature, controlling the molding of butter. Butter should 
be handled and stored, therefore, in such a manner that the minimum of 
condensation moisture will be formed. 


6. No relation was found between the mold count of butter and the sub- 
sequent development of moldiness. 


7. Of four methods of treatment, wet sterilization of the parchment was 
the least satisfactory, while parchment soaked in hot or cold saturated brine 
gave the maximum protection against subsequent mold growth. 


8. Cladosporium sp. was isolated from a moldy patch in butter showing 
interior molding. The point of initiation of growth was observed and is 


described. 


“ 
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WEATHER AND WHEAT YIELD IN WESTERN CANADA 


I. INFLUENCE OF RAINFALL AND TEMPERATURE DURING 
THE GROWING SEASON ON PLOT YIELDS! 


By J. W. Hopkins? 


Abstract 


A statistical study of plot yields recorded at a number of agricultural ex- 
periment stations in central and southern Saskatchewan and Alberta has demon- 
strated a significant correlation between yield and the amount and distribution 
of seasonal rainfall. On the whole, above-average rainfall is associated with 
higher yield, but the result of a given increment of rain at different times is 
partly dependent on soil conditions. On fertile soil, rainfall prior to harvesting 
results in a reduction of yield, probably owing to lodging. The maximum in- 
fluence of precipitation upon yield appears to be exerted during the month of 
June. The average summer rainfall sequence is very similar in each of the above 
four districts. There is a moderate degree of correlation between the amounts 
of rain recorded in different districts during the same season, but the simul- 
taneous occurrence of extremely wet or dry seasons over the whole area seems 
infrequent. 

Temperature conditions during the growing season seem to be secondary to 
rainfall in influencing yield. Above-average temperatures are beneficial at the 
time of sowing, detrimental during mid-summer and again beneficial prior to 
ripening, but as in the case of rainfall, the effect produced is influenced by soil 
conditions. No consistent relation is evident between either rainfall or tem- 
perature and the relative yield of early and late maturing varieties. 

It is apparent that the yields secured are influenced by factors other than 
those considered, and the precipitation during the autumn, winter and spring 
months prior to sowing is being studied in this connection. 


I. Introduction 


It is generally believed that the large fluctuations in the annual yield of 
wheat and other Canadian field crops are to a considerable extent attributable 
to variations in the weather conditions prevailing in different seasons. Apart 
from certain broad generalizations, however, little can be said to be definitely 
known at present concerning the actual quantitative relations involved and 
the extent to which these may be modified by soil conditions or cultural 
practices. 

Estimation of theseeffects iscomplicated by the fact that, during the growing 
season at least, equal increments of rainfall, heat, etc., occurring at different 
times cannot be assumed to exert a similar influence upon yield; for in addition 
to meteorological considerations (13), the requirements of the plant do not 
remain constant throughout the growth cycle (e.g., 6, 14, p. 67). Several 
agricultural meteorologists have therefore been led to postulate the existence 
of ‘‘critical periods” in crop development. At these times (which may be of 
relatively short duration) the plant is supposed to be particularly sensitive 
to an excess or defect of some meteorological element, of which fluctuations 


in the then available amount exercise a major influence upon yield (1, 2, 3, 
4, 7, 12, 17). 


1 Manuscript received February 7, 1935. 

Contribution from the National Research Laboratories, Ottawa. Abstracted from a thesis 
approved for the degree of Ph.D. in the Faculty of Science, University of London. Published as 
Paper No. 73 of the Associate Committee on Grain Research. 

2 Biologist, National Research Laboratories, Ottawa. 
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In this connection it was observed during the course of wheat varietal 
tests at the University of Alberta that in some years marked differences 
in the relative order of yield of early and late maturing varieties occurred, 
and it appeared possible that these fluctuations in relative yield were partly 
attributable to some differential effect of weather conditions upon the varieties. 
Knowledge of any such differential effects upon varietal yields would be of 
considerable practical importance. 

The present exploratory investigation was, therefore, undertaken with the 
object of studying the relation between rainfall and temperature conditions 
during the growing season and the yield of a specified variety of wheat, 
and also the relation between these weather factors and the relative yield of 
early and late maturing varieties. The method employed is the statistical 
one, in which numerical measures of the weather conditions are correlated 
with the associated crop yields with the object of deducing the quantitative 
relations between them. 


2. Observational Data 


Owing to the nature of the relations which it is sought to investigate, 
a relatively large number of statistical parameters must be calculated in 
order to specify them. This in turn necessitates an extensive series of observa- 
tions, if reliable results are to be obtained. Fisher, for example, investigating 
the influence of rainfall upon the yield of wheat at Rothamsted (10), utilized 
the records of 60 harvests from Broadbalk field. 

In Western Canada continuous crop records of this extent are naturally 
nowhere available. At the various Dominion Experimental Stations through- 
out the prairie provinces, however, and also at the agricultural colleges of 
the provincial universities, plot yields of the more commonly grown wheat 
varieties have been recorded for a number of years. An attempt has, there- 
fore, been made to combine the experience of a number of these stations in 
order to arrive at an estimate of the average relation in the area which they 
represent. This, being based on a more extensive body of observations than 
is available at any single station, may be expected to be correspondingly more 
reliable, though it must be recognized that the relations thus deduced are 
probably not exactly representative of the situation at any individual station. 

In Tables I and II the annual yields in the variety test experiments referred 
to are recorded, the values in each case representing the average performance 
of replicate plots. These were obtained from the annual reports of the 
directors of the various Dominion Experimental Stations and from records 
supplied by the Department of Field Husbandry, University of Saskatchewan, 
and Department of Field Crops, University of Alberta. Various members 
of the Experimental Farms Branch also assisted by supplying information 
not available in published form. 

These trials are commonly carried out on well prepared summerfallow land, 
the fertility of which is maintained by the rotation of leguminous or other 
crops. The investigation has, therefore, been extended to a further series 
of crop returns, obtained under less favorable conditions. This comprises 
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TABLE I 
YIELD, IN BUSHELS PER ACRE, OF MARQUIS WHEAT GROWN IN VARIETY TEST PLOTS 












54.6 47.2 17.0 
33.0 46.7 20.8 
29.0 43.7 27.8 
46.4 36.9 42.5 
59.9 40.5 26.7 
35.6 46.4 38.1 
44.6 55.9 25.1 
22.6 24.3 26.0 
33.6 54.3 12.2 
58.8 46.3 2.6 





Results in 1923 and 1924 are omitted, as yields at several stations were adversely affected by 
hail storms. 


1 University of Alberia. 
2 University of Saskatchewan. 
TABLE II 
YIELD, IN BUSHELS PER ACRE, OF WHEAT VARIETIES GROWN IN TEST PLOTS 








Station Variety 1931 | 1930 | 1929 | 1928 | 1927 | 1926 | 1925 
Edmonton Garnet 48.5 | 28.2 | 24.3 | 47.4 | 51.9 | 50.4 | 56.7 
Reward 40.3 | 26.5 | 18.3 | 45.6 | 46.9 | 46.3 | 45.8 

Red Fife 57.1 | 41.4 | 36.9 | 39.6 | 47.7 | 41.7 | 43.2 

Lacombe Garnet 67.2 | 38.8 | 41.1 | 41.3 | 48.0 | 48.8 | 42.0 
Reward 64.6 | 34.7 | 35.8 | 42.5 | 41.0 | 45.0 | 35.6 

Kitchener* 47.2 | 49.8 | 45.7 | 29.6 | 46.0 | 55.0 | 56.0 

Lethbridge Garnet 20.7 | 20.3 | 28.3 | 37.9 | 18.8 | 24.4] 25.0 
Reward 19.1 | 17.6 | 27.0 | 35.5 | 19.8 | 23.9 | 25.0 

Red Fife 15.8 | 25.0 | 22.5 | 36.1 | 29.6 | 27.2 | 34.9 

Indian Head Garnet 7.3 | 40.8 | 28.3 | 34.7 | 30.0] 57.5 | 48.8 
Reward 7.5 | 35.6 | 25.4 | 32.9 | 19.1 | 45.8 | 41.9 

Red Fife 19.8 | 29.8 | 32.5 | 37.3 | 36.1 | 49.1 | 41.2 

Swift Current** Garnet 6.9 | 25.7 | 20.1 | 45.7 | 39.1 | 31.1 | 26.5 
Reward 6.7 | 20.0 | 22.1 | 34.6 | 34.4 | 29.5 | 22.0 

Marquis 9.4 | 24.4 | 18.5 | 42.0 | 39.3 | 31.0 | 23.3 

Red Fife 9.4 | 24.2] 17.7 | 44.8 | 37.4 | 34.0 | 28.0 

Scott Garnet 29.4 | 43.9 | 15.7 | 42.8 | 44.0 | 30.9 | 38.4 
Reward 23.5 | 38.9 | 18.1 | 43.3 | 39.0 | 28.0 | 36.4 

Red Fife 30.4 | 38.4 | 17.3 | 33.1 | 42.3 | 22.7 | 33.7 

Rosthern Garnet 13.8 | 29.2 | 44.8 | 27.5 | 34.0] 20.0 | 28.3 
Reward 12.7 | 26.6 | 34.9 | 23.3 | 33.7 | 20.8 | 24.0 

Red Fife 18.3 | 29.9 | 50.9 | 28.8 | 28.0 | 18.3 | 33.3 





*A late maturing variety; Red Fife not grown. 
** Station established in 1921. 
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the annual yields (usually of Marquis, but in some cases of Reward or Garnet), 
listed in Table III, of the summerfallowed and stubble plots of a summer- 
fallow-wheat-wheat rotation. Owing to its superior cultivation, the soil of 
the variety test plots will undoubtedly differ both chemically and physically 
from that of these grain rotation plots. In particular it may be expected 
to have a higher moisture-holding capacity, an increased proportion of 
fibrous constituents and greater amounts of readily nitrifiable plant residues 
and of plant nutrients of all kinds. 


TABLE III 
YIELD OF WHEAT, IN BUSHELS PER ACRE, FROM SUMMERFALLOWED AND STUBBLE LAND 





Lethbridge | Indian Head |Swift Current Scott Rosthern 
Year Fal- |Stub-| Fal- |Stub-| Fal- |Stub-| Fal- | Stub-| Fal- | Stub- 

low ble low ble low ble low ble low ble 
1931 7 t S.E PUES) GS -& 0.2 | 33.0 | 12.3 | 34.0} 15.2 
1930 28.5 | 12.1 | 24.2 | 23.1] 19.0] 5.1 | 35.8 | 26.1 | 26.4] 23.1 
1929 26.2 | 25.6} 18.3 | 0.01 9.5} 4.51 7.77 §.8 1 22.1 | 17.3 
1928 §2.1 | 27.1 | 32.7 | 19.5 | 32.1 | 18.9 | 37.0 | 23.8 | 20.8 | 19.3 
1927 38.4 | 32.3 | 33.6 | 18.3 | 42.7 | 34.3 | 36.0 | 27.2 | 30.8 | 31.2 
1926 28.2 | 29.7 | 30.8 | 25.4] 20.5 | 12.5 | 18.8 | 10.9 | 30.4 | 20.0 
1925 25.4 | 12.5 | 29.2 | 25.1 | 20.5 | 18.5 | 32.9 | 29.9 | 34.8 | 32.8 
1922 34.5 | 18.6 13.3 | 12.2 
1921 31.6 | 24.3 24.3 | 24.0 
1920 14.6! 5.2 25.0 | 18.0 


The meteorological data used, consisting of daily observations of precipita- 
tion and temperature, were extracted from the Monthly Record published 
by the Meteorological Service of Canada. Certain gaps occurred in the 
printed records but the information represented by these was supplied by 
the Director of the Meteorological Service. 


3. Method of Analysis 


The considerations outlined in Section 1 necessitate the determination of the 
effect on yield of any weather factor, at different periods of the growth cycle; 
and if a search is to be made for “‘critical periods’’ of possibly short duration, 
the season must be fairly finely subdivided. Five-day intervals have been 
suggested for this purpose. The amount of precipitation, etc., in each period 
may then be regarded as constituting an individual variable, the effect of 
which upon yield may be determined by the methods of partial correlation. 

This direct analysis suffers, however, from serious disadvantages (12). 
The arithmetical work incidental to the formulation and solution of the 
necessary sets of simultaneous equations in 25 or more unknowns is extremely 
laborious. More important, however, is the consideration that owing to the 
high proportion of independent variables to observations, the standard errors 
of the regression coefficients will in all probability be so large as to destroy 
their significance. 

These difficulties were overcome by Fisher (10), who introduced the concept 
of the regression integral, involving a regression function varying continuously 
with the time. The practical application of this method of analysis consists 
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in expressing the amount and distribution of, for example, rainfall in each 
year by means of the six numerical coefficients of a fifth-degree polynomial 
function of the time; 

polo + pili + poT2 + pals + pals + psTs 
the six terms 7, of order 0, 1, 2. . . 5, being mutually orthogonal. The 
six coefficients p (found by Least Squares) corresponding to each year are 
then employed as independent variables with which the crop is to be cor- 
related, in place of the 25 or more actual values of rainfall by periods. Then 
by virtue of the orthogonal properties of the functions 7, employed, the 
regression integral of yield upon rainfall may itself be expanded as a fifth- 
degree polynomial function of the time involving the partial regression 
coefficients a of yield upon the distribution coefficients p, and the effect on 
the crop of an extra inch of rainfall at any period displayed. 

This method has been employed in the present study. But, whereas 
Fisher, dealing with an autumn sown crop, took into consideration the rainfall 
sequence throughout a period of 366 days ending on the 31st of August of 
the year in which the crop was harvested, attention has here been confined 
to a period of 125 days commencing 15 days before the date of sowing each 
year. Instead of the 61 six-day intervals into which he grouped the daily 
records of rainfall for the determination of the six coefficients p of each season, 
25 five-day intervals have been employed. Otherwise, the procedure is the 
same. The time of sowing, rather than any fixed calendar date, was used 
as a reference point because the former may vary considerably from year to 
year, and from district to district in the same year, in exceptional cases by 
as much as one month. 

In Table IV will be found the 462 rainfall distribution coefficients approp- 
riate to the various stations and seasons for which crop records are available. 


TABLE IV 
RAINFALL DISTRIBUTION COEFFICIENTS 


Station and year Po pl p2 ps pa ps 
Edmonton 1931 2.772 1.138 — .740 — .543 —.190 — .208 
1930 1.398 .273 — .588 — .413 .042 074 
1929 1.500 .450 —.119 — .784 — .153 090 
1928 2.170 . 504 — .849 — .649 . 134 660 
1927 1-996 — .204 —1.026 — .650 1.041 668 
1926 1.962 605 — .656 .304 .595 — .392 
1925 1.924 — .039 1.084 — .443 1.380 — .609 
1922 1.492 463 653 .420 476 225 
1921 1.894 1.087 — .358 — .065 — .038 —.227 
1920 2.246 — .250 — .862 .165 — .278 — .666 
Lacombe 1931 2:592 1.634 —.741 — .868 577 1.391 
1930 1.810 575 168 |“ —.270 — .542 — .451 
1929 1.012 098 019 . 208 261 275 
1928 2.368 577 —1.228 . 284 1.667 —.176 
1927 2.584 160 — .867 — .438 390 722 
1926 2.250 803 —.111 .962 328 — .285 
1925 1.072 140 — .361 .029 319 043 
1922 1.422 291 115 — .107 477 362 
1921 1.644 .387 —.128 — .243 — .443 . 188 
1920 1.018 —.251 —.149 .008 . 132 . 202 
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TABLE I1V—Concluded 


RAINFALL DISTRIBUTION COEFFICIENTS 
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The six numbers constituting each horizontal line specify the amount and 
distribution of rain in each season. As it was necessary to calculate 77 sets 
of six coefficients using the same series of values of the polynomial functions, 
six series of 25 multipliers, consisting of the values of the terms Jo, 71, . . . 7s 
for the successive periods of time were constructed as suggested by Tippett (20) 
so that multiplication of the amount of rainfall in each period by the corres- 
ponding term of the appropriate series and addition of the resulting products 
gave the distribution coefficient directly. This procedure facilitated the 
computation of the coefficients p, but even so the arithmetical work involved 
was considerable. 


4. Regression of Yield on Rainfall During the Growing Season 
Rainfall and Yield of Marquis Wheat 


The coefficients of regression ao, a1, . . . as of yield y on the rainfall dis- 
tribution coefficients po . . . . ps are to be obtained from the solution of the 
six simultaneous equations 

@eS(0g) + O,S(Pehi) + MSPs) +... ee ee + &S(PoPs) = S(Poy) 
oS (PoP1) + %S(P:*) + HS(PiP2) + 2. ee ee ee + %S(PiPs) = S(rry) 
oS (PoPs) + %S(P1Ps) + OS(s) +... - we ee + %S(e#) = Sosy) (1) 


the quantities p and y being measured from their respective means and the 
summations extending over the whole of the observational series. In order 
to eliminate any bias resulting from differences in average soil fertility or 
rainfall conditions between stations, the quantities to be correlated will not 
be the values of y, Po, P1, etc., actually observed, but the amount by which 
these are in excess or defect of their mean value at the appropriate station 
during the period under consideration. The sums of squares and products 
in Equation (1) will therefore be of the form 
SS(0rp — Prep) (Pop — Pap) 
where p,, is the mean value of p, at station p, rather than 
S(Pr — Pr) (Pe — Be) 
where p, is the mean of all values of p,, and will correspond to 63 of the total 
69 degrees of freedom. 
Computation of these quantities yields the following numerical values: 
16.916160a, + .704273a, — 7.374067a2. — -153974a3 + 3.366082a,4 
+ 2.677181a, = 129.51124 
— 5.327921a3 + .2387600.4 
+ 2.090910a; = —46.54435 
++ 4.591087a; — 3.442486a4 
— 2.087548a5 = —130.18403 
— .153974a9 — 5.327921a, + 4.591087a, + 16.260811a; + .035470a4 
+ 
+ 
+ 


.704273a9 +13.445682a, — .966389a2 


— 7.374067a) — .966389a, + 15.382760a, 


3.460086a, = 4.45425 
.035470a3 + 18.600138a4 
.981074a5 = 32.98832 

3.460086a; + .981074c%4 

13 .885332a5 


3.366082a0 + .238760a, — 3.442486a, 


2.67718lao + 2.090910a, — 2.087548, 


3.59866 (2) 
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giving: 
ao = 5.16397 a3; = .80204 
a, = —3.74079 ag = — .67708 
a = —6.74085 a = — .93891 


The significance of the correlation between yield and the rainfall coefficients 
may be investigated by determining the distribution of the total variance 
of yield, corresponding to 63 degrees of freedom, between the six degrees 
of freedom of the regression formula and the 57 degrees of freedom appropriate 
to deviations of the actual yield values from the regression formula. It is 


apparent from Table V TABLE V 
that the regression for- DISTRIBUTION OF SEASONAL VARIANCE OF YIELD OF 
mula has in fact account- MARQUIS WHEAT 


ed for more than a pro- 
portionate amount of the 
variance and a definite, 


Variance due to Degrees of} Sum of 
freedom squares 


Regression formula 6 1698.31 
though moderate degree Deviations from regres- 57 6370.15 
of relationship between sion formula 
the variables in question Total 63 8068.46 





is indicated. 


In certain circumstances it is conceivable that the pooling of the seven 
sets of sums of squares and products derived from the different stations 
might not be valid. The incidence of rainfall at the various stations during 
the period under review is not entirely independent; there is, as might be 
expected, a certain intra-annual correlation between stations with respect 
to both its amount and distribution. If then any extraneous source of varia- 
tion in yield is common to different stations in a given year, an error will be 
introduced into the regression formula deduced in the preceding paragraphs. 

This point was investi- 
gated by analyzing the 
total seasonal variance and 
covariance of yield and 
the rainfall coefficients (63 
degrees of freedom) into two 
portions, of which the first 
represented the average 
variance and covariance 
between years common to 
all stations (nine degrees of 
freedom) and the second 
represented the residual 
variance and covariance 
at each station each year 
(54 degrees of freedom). eisty ace eh tai al 


Two sets of the sums of Fic 1. Average effect of rainfall on yield of Marquis 
squares and products form- wheat (variety test plot data). 


QusHeLs PER ACRE) 





AVERAGE EFFECT OF EACH ADDITIONAL INCH OF RAIN 
. ‘ e 
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ing the numerical coefficients of Equations (1) were thus obtained, 
yielding two series of values of the coefficients a. The differences between 
these were, however, insignificant when considered in relation to the residual 
deviations; consequently, the original regression formula, derived from the 
total 63 degrees of freedom, is employed in subsequent calculations. 


As explained in Section 3, the values of a enable the average effect on the 
crop of an extra inch of rainfall at any time within the 125-day period to be 
estimated. Fig. 1 shows the course of the effect throughout that period. 


Rainfall and Differential Yield of Wheat Varieties 


For this study the crop yields of Table II are available, together with the 
corresponding yields of the variety Marquis from Table I. The quantities 
whose seasonal variation in excess or defect of their respective station averages 
is to be correlated with the rainfall coefficients are not the absolute yields, 
but the differences in yield between the varieties Garnet and Reward (both 
early maturing sorts), Garnet and Marquis (the latter medium-early maturing), 
Garnet and Red Fife (the latter late maturing), and Marquis and Red Fife. 
The resulting elimination of irregular variation due to differences in the 
fertility of the soil on which the trials were conducted in different years, 
common to both plots of any comparison, may be expected partially to offset 
the reduction in number of the observational series. 


Owing to the smaller number of crop records, and the inclusion of results 
from the Swift Current station, the following new series of numerical values 
of the sums of squares and products of rainfall coefficients in Equations (1) 


must be computed: 


12.441809 .447835 —4.232000 1.855544 2.699285 5.111341 
.447835 9.186527 —.422194 —3.748035 .565771 .077764 
—4.232000 —.422194 10.456264 .260742 —4.267647 —3.509379 
1.855544 —3.748035 . 260742 7.905838 —.955070 —1.870688 
2.699285 .565771 —4.267647 —.955070 16.362681 1.468324 
$.111341 .077764 —3.509379 —1.870688 1.468324 10.407955 (3) 


As it will eventually be necessary to solve equations involving the present 
set of values of these independent variates in conjunction with ten different 
dependent variates, it will be advantageous to calculate at this stage the 
values of 

Cre = ey 
where 4,, is the co-factor of the element in the r” row and s” column of the 
determinant formed from the matrix of the numerical coefficients of Equations 
(3). These c values, determined in the manner outlined by Fisher (11, Sec. 
29) are given in Table VI. The partial regression coefficient a, of any depend- 
ent variate y on p, may then be rapidly obtained from the relation: 


Se = BAA) FGPG 4H ook cece ew eee ss + 6,5(,y) (4) 








.127657 
— .028570 
.026392 
— .060570 
— .011081 
— .062904 
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TABLE VI 
VALUES OF Cre 
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— .028570 


. 143207 
.002755 
.081472 
.002704 
.028120 


C.2 c.3 
.026392 — .060570 
.002755 .081472 
. 125698 .000207 
.000207 . 195697 
.026035 .013039 
.025766 .062541 
TABLE VII 


006697 
"145746 


SUMS OF PRODUCTS OF VARIETAL YIELD DIFFERENCES AND RAINFALL COEFFICIENTS 


Garnet-Reward | Garnet-Marquis 


Garnet-Red Fife 


Marquis-Red Fife 


S(psy) 
S(pay) 
S(psy) 


15.6465 


—8.3329 


12.3720 
4.4511 


—1.7251 


1.2907 


22.7739 
7.2119 

-- 24.0684 
11.9291 
42.6996 
8.3294 


5 
—3 


8 
3 


2.2249 
4.0452 
7.8503 
4.9571 
1.0705 
1.2984 


32.1449 
—3.7745 
—41.7687 
—12.4828 
33.2225 
33.2491 


Utilizing the differences in yield of the above-mentioned wheat varieties 
as successive dependent variates, the sums of products in Table VII and the 
following regression coefficients result: 


Difference in yield of Garnet and Reward 


Difference in 


Difference in 


Difference in 


ao 
ay 
‘Qe 


ao 
a@ 
a2 


ao 
a) 
a2 


ao 
ay 
a2 


2.23029 a; = 
— 1.21200 “a = 
1.93438 as = 

yield of Garnet and Marquis 
. 34679 a; = 
1.63877 a = 
— 1.07564 a = 

yield of Garnet and Red Fife 
2.98575 a = 
— .31845 a = 
— .45215 a = 

yield of Marquis and Red Fife 

1.40547 a3 = 
—1.56621 a = 
— 2.69321 a = 


— .69475 
.07248 


= — .44483 


2.61538 
2.33236 
. 39608 


= —.79706 


4.25054 
.64787 


— 2.19339 
.92497 
1.08334 


The distribution of the seasonal variance of the varietal yield differences 
between the six degrees of freedom of the appropriate regression formula and 
the remaining 36 degrees of freedom is shown in Table VIII. 
the proportion of the total variance which can be expressed by means of the 


In each case 





BUSHELS PER ACRE 
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regression formula is low, and in only two instances does the regression mean 
square exceed the mean square deviation from the regression formula. Both 
of these involve the late-maturing variety Red Fife, but the actual ratios 
of the mean squares correspond to z values of only .216 and .279, respectively, 
whereas a value of .430 or more might be expected to occur on the average 
once in 20 times as the result of chance differences alone. It must be con- 
cluded therefore that the observed fluctuations in relative yield are not asso- 
ciated in any consistent manner with the rainfall conditions prevailing during 
the growing season. 
TABLE VIII 


DISTRIBUTION OF SEASONAL VARIANCE OF VARIETAL YIELD DIFFERENCES 









Variance Garnet- 
Reward 


Garnet- | Garnet- | Marquis- 
Marquis | Red Fife | Red Fife 


5, a (6d f freedom) 65.14 
ue to regression legrees o om ‘ 
Deviations from regression (36 degrees of | 518.18 


freedom) 
Total (42 degrees of freedom) 583.32 





179.69 540.58 257.71 
1721.79 | 2104.05 884.94 


1901.48 | 2644.63 | 1142.65 


Mean square: 
Due to regression 10.86 


90.10 42.95 
Deviations from regression 14.39 


58.45 24.58 


Rainfall and Yield of Wheat from Summerfallowed and Stubble Land 


The crop data to be used in this study are those of Table III. Inspection 
reveals that the level of yield is definitely higher on the fallowed than on the 
stubble plots, but that both are inferior in this respect to the variety test plots. 
This point is illustrated in Fig. 2, in which the mean yield and seasonal standard 
deviation of both series are 
shown, together with the corres- 
ponding values deduced from 
the Marquis yields of Table I. 
Inspite of the increased average 
YY), yield, the annual fluctuations 
of the summerfallowed are 
. g UY no less, but rather slightly 

YY greater, than those exhibited 
by the stubble land crops. 
The extra moisture and plant 
nutrients accumulated as a 

pug Y Us result of the fallow are re- 
Yj Yj Wha flected in the higher level of 

GRAIN ROTATION «GRAIN ROTATION «= VARIETY TEST yield, but apparently have not 

ne re re ere sufficed to exert any absolute 
Fic. 2. Average yield nh -hatched column) and stabilizing effect, though the 


annual variation (cross-hatc in bushels per acre of oe : 
three series of wheat plots. percentage variation is lower. 
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The numerical values of the sums of squares and products of rainfall 
coefficients in Equations (1) corresponding to the years and stations for which 
crop yields are available in Table III are: 


12.843475 — .335858 —4.365599 1.826327 1.794788 4.455203 
— .335858 8.950589 —.005577  —4.254751 1.254894 . 280159 
—4.365599 — 005577 9.856559 2.416812 —4.985248 — .724501 
1.826327 —4.254751 2.416812 10. 123139 — .217015 — .367921 
1.794788 1.254894 —4.985248 — .217015 15.041716 1.515575 
4.455203 . 280159 —.724501 — .367921 1.515575 10.644241 (5) 
From these, the quantities c,, shown in Table IX may be obtained in the 
manner already indicated. Multiplication of the six values in each column 
TABLE IX 
VALUES OF Crs 


. 123494 — .018813 .070191 ‘ 5 .050149 
— .018813 . 149510 — .037348 F 4 .007155 
.070191 — .037348 . 175304 ‘ ‘ .026646 
— .048457 .074932 — .070012 : ; .022308 
-014451 — .022246 .054516 ‘ .014792 
— .050149 -007155 — .026646 -022308 — .014792 - 115813 


of Table IX by the corresponding six sums of products from Table X, as in 
Equation (4) page 314, yields the following regression coefficients of yield 
on the rainfall coefficients: 


Yield from stubble land: 
a= 9.09916 . 83645 
a = —2.34024 — 1.61300 
ae — 2.37545 = .03891 


Yield from summerfallowed land: 


9.80077 as .00891 
— .48082 a, = —.70476 
— 3.87536 a= .19378 


Difference in yield from summerfallowed and stubble land: 


a = . 70161 a3 = —.82754 
a 1.85942 a= .90824 
a = —1.49991 a= .15488 


It will be observed from TABLE X 
: SUMS OF PRODUCTS OF YIELD AND 
Table XI that an appreciable qastienes commeennienn 
proportion of the annual vari- © ——————________ 
ance of yield of both the Stubble ——- ae ae 
stubble and summerfallowed yields yields yield 


plots is accounted for by the 


: ; 126.7605 | 142.5702 15.8097 
appropriate regression for- S(p,) —29.6080 —8.4421 21.1659 
mula. The ratio of the mean —52.9048 | —77.5864 | —24.6816 


’ 29.6453 10.7511 — 18.8942 
square due to the regression 2946 25.9970 25.7024 


to the mean square deviation 5( 39.2906 | 47.3286 8.1190 





AVERAGE EFFECT OF EACH ADDITIONAL INCH OF RAIN (US PER acAc) 
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TABLE XI 
DISTRIBUTION OF SEASONAL VARIANCE OF YIELD FROM SUMMERFALLOWED AND STUBBLE LAND 


Summer- Difference 
Variance Stubble fallow i 
yields yields 
Sum of squares: 
Due to regression (6 degrees of freedom) 1374.22 1692.98 
Deviations from regression (36 degrees of | 2331.27 2498 .29 
freedom) 
Total (42 degrees of freedom) 3705.49 4191.27 
Mean square: 
Due to regression 229.04 282.16 
Deviations from regression 64.76 69.40 





corresponds to a.z value of 
.632 and .701 in the case of 
the stubble and fallow series, 
respectively, both being in ex- 
cess of the 1% point, .623. 
The regression of the differ- 
ence in yield in the various 
seasons between the summer- 
fallowed and stubble plots 
on the rainfall coefficients is, 
however, totally insignificant. 
Fig. 3 shows the course of 


Fic. 3. Average effect of rainfall on yield of wheat the two regression functions 


from stubble plot (solid line) and summerfallowed plot throughout the season. 
(broken line) of grain rotation. 


Discussion of Results 


The results as a whole may be held to reflect clearly the sub-optimal moisture 
conditions prevailing during the average western Canadian season. Certainly 
the curves in Figs. 1 and 3 differ markedly from the corresponding portions 
of those obtained by Fisher (10) for winter wheat and by Wishart and Mac- 
kenzie (21) for spring-sown barley under English conditions, with which they 
may be compared by taking the average date of seeding to be May 1. 

The detrimental effect of rain later than 90 days after sowing, a pronounced 
feature of Fig. 1, is not reproduced in either of the curves in Fig. 3 and may 
probably be ascribed to loss of crop through lodging of the longer-strawed 
plants produced on the more fertile variety test plots. It is interesting 
to note that Fisher found the dunged plot 26 to differ from all the other 
Broadbalk plots in the severe damage caused by rain immediately preceding 
harvest. A similar adverse effect manifests itself on all the Hoos field barley 


plots receiving nitrogenous dressings (21), but is biggest and earliest in the 
case of the farmyard manure plot. 
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It should be emphasized in any discussion of Figs. 1 and 3 that the curves 
indicate not the direct effect of so much rain, but the resultant of the entire 
combination of weather conditions, such as low temperature, high relative 
humidity, lack of sunshine, etc., associated with rainfall at any given time. 
Thus the detrimental effect on the varietal plot yields of rainfall prior to 
sowing may be due largely to the accompanying reduction in temperature. 


The lack of correlation between varietal differences and the rainfall coeffi- 
cients does not necessarily mean that the former are unaffected by weather 
conditions, but does indicate that there is no consistent relation between the 
observed differences (which are sometimes large) and the rainfall sequence 
during the growing season. Such inconsistencies may be attributable to 
any or all of a variety of factors, e.g., other meteorological phenomena not 
associated with rainfall, weather conditions (including precipitation) during 
the preceding summer, autumn and winter, and soil differences both within 
and between stations. 


Considering the season as a whole, both the summerfallowed and stubble 
plots of the grain rotation show a greater response to rainfall than do the 
variety test plots. This is particularly the case in the later stages of growth. 
The enhanced moisture-holding capacity of the more fertile soil of the variety 
plots may result in a larger reserve supply of moisture still available at this 
time to meet the requirements of the crop, which, as already noted, was found 
by Briggs and Shantz (6) to attain its maximum rate of transpiration prior to 
the onset of ripening. It is well known, moreover, (e.g., 15, 18) that the 
degree of fertility of the soil has a marked effect on the efficiency of trans- 
piration, with the result that the water requirement on a poor soil is con- 
siderably higher than on a fertile one. There is some indication (Fig. 3) 
that the stubble plots of the grain rotation are more responsive to wet weather 
in the early part of the season and that the fallowed plots are able to utilize 
later rains more advantageously, perhaps owing in the latter case to a more 
extensive plant development in the preceding period made possible by greater 
reserves of moisture and nutrients. As, however, these differences are not 
at any time pronounced, and are statistically quite insignificant, such con- 
siderations cannot be more than conjectural. 


The maximum benefit from additional rain on the variety plots is obtained 
approximately one month after the date of sowing. About this time tillering 
is taking place and a period of active growth is being initiated. It is possible, 
therefore, that additional moisture in the surface level induces a more prolific 
vegetative growth, which the fertile soil is capable of reflecting in increased 
yields. Smith (19), for example, found from 8000 observations on 64 varieties 
of spring wheat, oats and barley grown in a cereal nursery’at Dickinson, N. 
Dakota, a high correlation between rainfall from May 16 to July 15, extent 
of tillering, and yield in different seasons, and concluded that under local 
conditions there was a close relation between rainfall and tillering and between 
tillering and yield. 
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The degree of correlation between seasonal yield and the rainfall coefficients 
is only moderate, being represented by values of r = 0.46, 0.60 and 0.62 in 
the case of the variety test (Marquis), stubble and summerfallow series, 
respectively. It should be recognized that the residual variance, not accounted 
for by the regression formulas, may arise not only from the operation of all 
other meteorological factors not associated with rainfall, and of weather 
conditions antecedent to the growing season, but also from variations intro- 
duced by growing the crop on different land each season (particularly in the 
case of the variety plots) and from dissimilarities between the true values of 
the regression function at different stations. Such variance may thus in part 
reflect the imperfections of the data available rather than an actual absence 
of causal relation. 


5. Characteristics of Seasonal Rainfall 


The extent to which the annual yields are influenced by the tainfall con- 
ditions in any particular period will depend not only on the effect of a given 
increment of rain at that time, but also on the variability of the amount of rain 
actually falling during this period from year to year. A study has therefore been 
made of the seasonal rainfall at several meteorological stations in central and 
southern Saskatchewan and Alberta, the results of which are presented graphi- 
cally in Fig. 4. With the exception of Lacombe (no records prior to 1908), data 
for the years 1900-1931, inclusive, were available from the following stations: 

Central Alberta: Edmonton, Lacombe. 

Central Saskatchewan: Battleford, Saskatoon. 

Southern Alberta: Calgary, Lethbridge. 

Southern Saskatchewan: Indian Head, Swift Current. 


SOUTHERN ALBERTA 





Fic. 4. Upper portion: Average precipitation (inches of rainfall) recorded at meteorological 
stations in central and southern Alberta and Saskatchewan. Below: Annual variation (standard 
deviation in inches) in amount of rainfall. 
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The upper portion of Fig. 4 portrays the average rainfall in each of the five 
months April-August, inclusive, recorded at the stations in these districts. 
The average rainfall sequence is similar in all districts, progressing from the 
low value of one inch or less in April to a maximum of approximately three 
inches in June. Seasonal fluctuations in the amount of rain received in the 
various months are measured by the standard deviation of the monthly 
totals, shown in the lower part of the diagram. These also follow a similar 
course in all districts except southern Alberta, where the variability of May 
rainfall exceeds that of June and July, otherwise uniformly the most variable 
months. The uncertainty of the rainfall in these regions is well illustrated 
by the high ratio of the standard deviations to the corresponding average 
monthly totals. 


Referring to Figs. 1 and 3 and taking the average date of sowing to be May 1, 
the beneficial effect of additional rain on both the variety test plots and the 
fallowed plots of the grain rotation is seen to be greatest during the month 
of June. As this month also exhibits the greatest variability in amount 
of rain, June precipitation appears to exert a greater influence upon the yield 
of these plots than that of any other month. Seasonal differences in July 
and August rainfall may also affect appreciably the yield secured from the 
variety test and summerfallowed grain rotation plots, respectively. May 
rainfall is in three out of the four cases more constant from year to year but, 
if deficient, may be expected to reduce the yield of the stubble plots and of 
the variety test plots. 


The seasonal incidence of precipitation in the four districts may be investi- 
gated by averaging the total amounts recorded during the five months at 
the two meteorological stations in each district. The mean values and 
standard deviations of these averages are: 


Central Southern Central Southern 
Alberta Alberta Saskatchewan Saskatchewan 
Mean (inches) 12.3 10.6 9.2 10.3 


2.76 4.03 2.57 2.71 


Table XII gives the coefficients of correlation calculated from the amounts 
of rain recorded in different districts during the same season. Five of the 
six coefficients are signifi- TABLE XII 

cant, showing that the rain- Copeeiation OF SEASONAL RAINFALL IN FOUR DISTRICTS 
fall in the four districts Score 
during a given season is not 

independent. The degree 

of association indicated, 

however, is not as high as oe 
might perhaps have been Southern Alberta 
anticipated between the Central Saskatchewan 


precipitation in adjacent res We bed of aathnane = 2 
parts of a comparatively we of t for \ 1 level of significance = .45 


Standard deviation 
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unbroken area such as the Canadian plains, with the result that the 
simultaneous occurrence of extremely dry or wet seasons over the whole 
area may be expected to be infrequent. 


6. Additional Effect of Temperature 

The regression of wheat yield on rainfall having been found, the partial 
regression of yield on temperature, eliminating all associated rainfall effect, 
may be determined in the manner outlined by Tippett (20). The tem- 
perature sequence of each season is first expressed as a continuous function 
of the time, using the orthogonal polynomial functions of zero, first, second, 
etc., order in time already employed in the reduction of the rainfall data. 
In this way a series of sets of six numerical coefficients 00, 6:;, . . , 9s, rep- 
resenting the temperature conditions of different seasons, is obtained. The 
regression of each of the temperature coefficients on the six rainfall coefficients 
is next found, and the residual deviations of the temperature. coefficients 
are then correlated with the corresponding residual deviations of wheat yield 
from the rainfall regression formulas of Section 4. 

As an index of the prevailing temperature conditions the daily means 
estimated from maximum and minimum thermometer readings have been 
used. The imperfections of such estimates in comparison with those to be 
obtained from a continuous record, or from thermometer readings at regular 
intervals, is well known, but more precise observations are not available. 

From the observed daily temperatures the numerical values of the coeffi- 
cients 8 were computed by the procedure used in obtaining the rainfall 
coefficients of Section 3, and are listed in Table XIII. Owing to the inequality 
of the number of seasons in the three series of crop yields, the regression of 
the temperature coefficients on the rainfall coefficients of each series must be 
determined separately. 


TABLE XIII 
TEMPERATURE DISTRIBUTION COEFFICIENTS 


Station and year 


Edmonton 1931 
1930 
1929 
1928 
1927 
1926 
1925 
1922 
1921 
1920 


Lacombe 1931 
1930 
1929 
1928 
1927 
1926 
1925 
1922 
1921 
1920 








Station and year 


Lethbridge 


Indian Head 


Rosthern 


Swift Current 


1931 
1930 
1929 
1928 
1927 
1926 
1925 
1922 
1921 
1920 


1931 
1930 
1929 
1928 
1927 
1926 
1925 
1922 
1921 
1920 


1931 
1930 
1929 
1928 
1927 
1926 
1925 
1922 
1921 
1920 


1931 
1930 


1929 . 


1928 
1927 
1926 
1925 
1922 
1921 
1920 


1931 
1930 
1929 
1928 
1927 
1926 
1925 
1922 
1921 
1920 


1931 
1930 
1929 
1928 
1927 
1926 
1925 
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TABLE XIII—Concluded 
TEMPERATURE DISTRIBUTION COEFFICIENTS 
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Temperature and Yield of Marquis Wheat 


In order to determine the regression of the temperature coefficients on the 
rainfall coefficients, the sums of products of the seasonal deviations of these 
two sets of quantities from their average values at the various stations must 
be found. These sums of products are given in Table XIV. Together with 
the corresponding sums of squares and products of the rainfall coefficients in 
the left-hand side of Equations (2), they provide the necessary data for the 
calculation of the required regression coefficients, which have the values 


indicated in Table XV. 


TABLE XIV 
SUMS OF PRODUCTS OF SEASONAL DEVIATIONS OF RAINFALL AND TEMPERATURE COEFFICIENTS 
(SER1Es I) 
6o a 62 3 6 4s 
Po 428.7280 |—1239.4226 | —650.8052 366. 5496 —11.0012 164.9752 
pi —177.4032 135.4570 291.2175 156.0703 | —217.5229 12.7797 
p2 — 150.8797 746.1933 307.1580 160.0163 84.3347 — 63.3566 
pa 78.4404 | —535.4481 129.0684 | —124.1476 | —150.7129 14.8424 
pa 145.7326 | —322.5385 | —510.7768 275.8483 | —525.8077 | —226.1613 
ps — 167.5631 248 .8022 316.1284 —89.4370 625.4040 | —500.9116 
TABLE XV 
REGRESSION OF TEMPERATURE COEFFICIENTS ON RAINFALL COEFFICIENTS 
(SER1Es 1) 
9 A 6, Os 64 6 
Po 29.83306 | —63.22830 | —48.55154 39.03065 | — 4.73207 14.47272 
pi —13.30091 —3.53529 28 .07211 7.77900 | —27.85799 3.67032 
p2 1.59389 33.03970 . 79859 29. 12029 18.11712 2.96908 
ps —3.12282 | —38.48140 23. 26893 | —15.80847 | —13.48704 —7.536429 
ry) — 3.62326 . 22028 21.92718 | —14.95925 27 .58732 15.00287 
Ps —16.09770 26.00405 32.27009 | —13.64222 47.55967 | —41.90966 


It is not necessary to calculate the actual values of the individual deviations 
in order to determine the regression of the yield residuals on the temperature 
residuals. The required sums of squares and products of these quantities 
may be obtained from the corresponding values prior to the fitting of the 
rainfall regression by means of the relations: 


S(6, —BrPo —B,P, WF es ig Se Oe rsPs)* = S(6,?) — Bw S(6,P,) ~6,5(6,P,) Re 
° *- —B,,S(6,P,) (6) 
S(6, — BrP. — Bani .- . —BrsPs } (6, — BP —BuP, Pe oP ee —BuPs) 
= S(6,9,) —BroS (9,0) —B,,S(O,p6) ee 6S. aoe ae — BrsS(6,p5) (6a) 


where 6,, is the regression coefficient of 6, on p,. 
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The seasonal variance and covariance of yield and the temperature coeff- 
cients provide the following sums of squares and products: 
60 Ch 02 93 % 
421227.6 —296405.0 —138159.2 —1075.5 —18083.1 
—296405.0 524177.7 85035.5 87427.9 30509. 1 
—138159.2 85035.5 295441.8 —112294.4 42250.4 
9 a 
1 8 
7 


Yield 

8887 .49 
— 23447 .17 

2117.49 
— 1418.98 
— 1673.04 
— 4472.18 


5 
29596. 2 
—31833.7 
—45360.9 
— 3234.2 
—58200.9 
160786. 1 


—1075.5 87427.9 —112294.4 305773.5 —106631.8 
— 18083. 1 30509. 42250.4 —106631. 260322 .8 


29596.2 —31833. —45360.9 —3234.2 —58200.9 


Adjusting these by means of (6) and (6a) gives the following equations to 
determine the regression coefficients y of the yield on the temperature residuals: 


404393. 8yo 
— 257652. 5 yo 
—111718. 570 
— 11017. 970 
— 13090. 875 


16108 . 8yo 


— 257652. 57: 
+394632.4y; 
+ 31626.07; 
+ 103356. 74; 
+ 4310.67 


+ 5617.4) 


~111718. 542 
— 13090.8y,4 
+ 31626.0y-2 
+ 4310.67, 
+ 253418 . 92 
+ 42515 .4y, 
— 99390.7y2 
—~ 108490. Ov, 
+ 42515.4y2 
+ 234522 Ov, 
— 16038. 272 


— 11017.9y; 
+ 16108. 8y; 
+103356.773 
+ 5617.45 
— 99390. 77s 
— 16038. 245 
+ 286537. 073 
— 18409. 5ys 
— 108490. 073 
— 25653.8ys 
— 18409. 


4871.30 


— 11065 .47 


8485 .58 


— 1446.99 


— 1509.39 


Sy3 


— 25653.8y7, +141051.0y; — 6338.72 


The values of ¥ satisfying these equations are: 
Yo = .0064312 
‘1 — .0306698 
v2 .0445184 


.0150464 
— .0108923 
— .0394075 


vs = 
vA 
% = 


The effect on the crop of a given increment of temperature at different 
times throughout the season may be deduced from these coefficients as ex- 
plained in Section 3, and is illustrated in Fig. 5. Table XVI shows the 
proportion of the variance of the yield residuals expressible by the temperature 
regression formula. This is seen to be only moderate, and even though the 
greater part of it is in fact attributable to the terms of first and second degree, 
the significance of the estimated effects is in some doubt. 


TABLE XVI 
DISTRIBUTION OF SEASONAL VARIANCE OF YIELD OF MARQUIS WHEAT 


Sum of 
squares 


| Degrees of | 
_freedom | reedom 


Variance due to 


Rainfall regression formula 
Residual temperature regression formula 
Deviations from segression formulas 


1698.31 
992.93 
5377.22 


Total 


{ 
} 
= 
| 
| 
oa I lia 
' 
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Temperature and Differential Yield of Wheat Varieties 


The seasonal covariance of the rainfall and temperature coefficients during 
the period for which varietal comparisons are available is productive of the 
sums of products listed in Table XVII. Multiplication of these by the 
successive columns of ¢ values of Table VI as indicated in Equation (4), 
page 314, enables the regression of the seasonal fluctuations of the temperature 
coefficients on those of the rainfall coefficients to be readily determined, the 
numerical values being given in Table XVIII. Tables XVII and XVIII 
enable the sums of squares and products of the seasonal deviations of the 
temperature coefficients to be adjusted to give the corresponding functions 
of the residuals, required for the determination of the regression of the yield 
residuals on the temperature residuals. The unadjusted values are: 


9 A 62 Os M% 8s 
8 221797.8 —124098.4 -—51633.4 —31824.8  —14355.6 12472.0 
A — 124098 .4 269063 .8 30872 .4 79153.2 55718.1 —36307.3 
Os —51633.4 30872 .4 184827.8  —98700.8 32392.8  -—12016.8 
63 —31824.8 79153.2 —98700.8 271426.7 —149209.5 25215.5 
04 —14355.6 55718.1 32392.8 —149209.5 252250.8  —54108.7 
6s 12472.0 -—36307.3 —12016.8 25215.5 —54108.7 99198 .7 


The residual sums of squares and products obtained from these by means of 
(6) and (6a), page 324, are: 


60 4, 62 43 % 4s 
Oo 204048 .9 —96491.7 —27884.0 —45544.9 6142.4 — 3898.0 
6 —96491.7 191653.0 10222.5 74448 .6 16047 .2 —5523.7 
02 — 27884 .0 10222.5 131552.5 —48547.1 —9461.7 7645.8 
6; —45544.9 74448 .6 — 48547 .1 174538 .0 —75222.7 9732.0 
% 6142.4 16047 .2 —9461.7 —75222.7 128662.1 —22180.9 
5 —3898.0 —5523.7 7645.8 9732.0 —22180.9 72142.3 
TABLE XVII 


SUMS OF PRODUCTS OF SEASONAL DEVIATIONS OF RAINFALL AND TEMPERATURE COEFFICIENTS 
(SErtEs IT) 


T) A 02 3 % 4s 
po 375.9597 | —633.5518 | —548.8587 220.6140 —55.3437 237.2455 
pi 86.8920 | —279.2813 163.4911 31.9991 | —207.7983 45.3690 
pe —53.2085 366.5013 198.1474 | —318.5033 478.7752 43.7068 
pa 92.6571 | —422.6514 —29.3746 | —409.6782 —10.3239 230.1808 
ms 216.3768 | —227.0621 | —454.6179 279.1454 | —665.9853 —36.7730 
ps —27.5618 14,9081 56.4651 | —320.2172 689.6705 | —320.8584 


TABLE XVIII 


REGRESSION OF TEMPERATURE COEFFICIENTS ON RAINFALL COEFFICIENTS 
(Series IT) 


80 
Po 37 .83086 
pL 8.91464 
p2 8.41584 
ps 3.52677 
ps 10.81043 


ps — 19. 34992 


—36.04711 
—55.51375 
22 .96345 
— 69 .04345 
—5.45320 
15.66161 


—66.24214 
37 .60539 
48301 

38 . 46038 
—20.05465 
47.57602 


60. 70666 
—44.22227 
—35.19257 

—107.38116 
1.34768 
—91. 60699 


—23.87033 
—10.10621 
58 .57674 
18. 95003 
—29 48537 
105 . 38554 


36.79157 
9.47086 
2.70329 

13 .83466 

—3.08344 
—45 .13542 
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In the same way, the sums of products of yield and temperature residuals, 
shown in Table XX, may be obtained from the unadjusted quantities of 
Table XIX. The regression of the four series of yield-difference residuals on 
the temperature residuals may then be deduced, the regression coefficients 
y taking the values shown below: 

TABLE XIX 
SUMS OF PRODUCTS OF VARIETAL YIELD DIFFERENCES AND TEMPERATURE COEFFICIENTS 












Garnet-Reward | Garnet-Marquis | Garnet-Red Fife | Marquis-Red Fife 









1516.18 












S(Ay) —1226.75 530.76 —4553.94 —4322.18 
S(6xy) —2176.92 —212.75 —582.83 —492.84 
S(O) 549.25 1150.77 800.51 420.69 
S(y) 1574.68 —4817 .03 — 8068 .49 — 3050.12 

1063.25 — 1808.47 —2956.67 — 663.46 










TABLE XX 
SUMS OF PRODUCTS OF VARIETAL YIELD DIFFERENCE RESIDUALS AND TEMPERATURE RESIDUALS 









Garnet-Reward | Garnet-Marquis | Garnet-Red Fife | Marquis-Red Fife 
















—3975.51 
S(hy) —1131.74 3231.43 — 1947.62 —3615.25 
S (Oey) — 1100.29 1037.59 3069.95 1363.18 
S(6sy) 264.81 1226.46 557.25 — 1506.94 
S(0y) 868.01 — 2635.16 —5377.19 —2162.12 
S(6sy) 524.42 — 2307.24 — 3083.17 78.36 










Difference in yield of Garnet and Reward 







yo =  .0090922 vs = .0143016 
v1 = —.0080484 vs = .0174913 
v2 . 0001092 vs = .0105813 







Difference in yield of Garnet and Marquis 
yo = —.0151360 ys = —.0254974 
yun = .0221254 ys = —.0449486 






v2 — .0073164 ys = —.0407105 
yield of Garnet and Red Fife 









Difference in 






yo = —.0048636 vs = —.0114844 
v1 = —.0061265 ys = —.0567342 
v2 =  .0180241 ys = —.0612737 





Difference in’ yield of Marquis and Red Fife 

, vo = 0133333 ys = —.0079812 
v1 = —.0079323 ve = —.0213125 

¥2 .0096345 ¥ = — .0052979 
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The proportion of the yield variance accounted for by the respective regres- 
sion formulas is shown in Table XXI. This is uniformly low except in the 
case of the difference Garnet-Reward. Even this, however, by no means 
attains the level of significance. It cannot be said, therefore, that these 
data reveal any consistent association between temperature conditions during 
the growing season and varietal differences in yield. 


TABLE XXI 
DISTRIBUTION OF SEASONAL VARIANCE OF VARIETAL YIELD DIFFERENCES 





: Garnet — | Garnet— | Garnet— | Marquis— 
Variance Reward | Marquis | Red Fife | Red Fife 
Sum of squares: ° 
Due to rainfall regression (6 degrees of 65.14 179.69 540.58 257.71 


freedom) 


Due to temperature regression (6 degrees 52.09 305.18 557.24 145.78 
of freedom) 


Deviations from regression (30 degrees of | 466.09 1416.61 1546.81 739.16 
freedom) 


Total (42 degrees of freedom) 583.32 1901.48 2644.63 1142.65 


Mean square: 


Due to rainfall regression 10.86 29.95 90.10 42.95 
Due to temperature regression 8.68 50.86 92.87 24.30 
Deviations from regression 15.54 47.22 51.56 24.64 


Temperature and Yield of Wheat from Summerfallowed and Stubble Land 


The sums of products due to the seasonal covariance of the rainfall and 
temperature coefficients will be found in Table XXII, and the regression of 
the temperature coefficients on the rainfall coefficients, calculated from these 
sums of products by means of the c values of Table X, in Table XXIII. 
Tables XXIV and XXV show the sums of products of yield and temperature 
coefficients, and of the yield and temperature residuals, respectively. 


TABLE XXII 


SUMS OF PRODUCTS OF SEASONAL DEVIATIONS OF RAINFALL AND TEMPERATURE COEFFICIENTS 
(Series III) 

























6 6, cP % 
¢ 

Po 441.8201 | —788.6703 | —701.4976 195.1780 144.0828 184.3352 
pr —422.1688 40.4009 233.1662 30.9984 | —217.2446 56.8434 
p: — 161.5960 315.0244 276.9642 | —309.6751 527.6116 80.6117 
ps 36.2710 | —444.0487 39.4245 | —235.5485 86.5130 51.6520 
ps 8.5924 | —289.1353 | —429.1750 357.2893 | —553.9455 | —156.0151 

9.1924 43.2061 181.6013 | —273.6773 590.4029 | —302.0484 
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TABLE XXIII 


REGRESSION OF TEMPERATURE COEFFICIENTS ON RAINFALL COEFFICIENTS 
(Series III) 









| af ff 














po 49.06755 | —60.87200 | —88.79662 32.08566 17. 10857 37 .74296 
pr —62.80296 | —17.41988 $1.51523 | —15.02830 | —31.86609 7.20007 
p2 16. 13444 12.53382 | —41.02141 1.51637 58. 73163 20.87404 
ps —36.06837 | —42.29236 52.34389 | —36.77078 | —20.41775 —5.30237 
ps 6.70989 | —10.06720 | —40.68117 25.53857 | —22.64785 —4.35172 






—19.12481 30.82053 57.72701 | —43.54940 55.66114 | —42.50655 










TABLE XXIV TABLE XXV 


SUMS OF PRODUCTS OF YIELD AND SUMS OF PRODUCTS OF YIELD AND 
TEMPERATURE COEFFICIENTS TEMPERATURE RESIDUALS 




















‘ Difference, Difference, 
ummer- | summer- Summer- | summer- 
Stubble fallow fallow — Stubble fallow fallow — 









stubble stubble 











S(oy) 11092.36} 12419.15) 1326.79 S(Boy) 5683.49} 7263.69) 1580.20 












S(Ay) |—15506.33) —17499. 16} —1992.83 S(@y) | —7583.85) —8737.51) —1153.66 
S(@xy) | —7682.48) —7037.33 645.15 S(02y) —828.16 685 . 33 1513.49 
S(@sy) | —1904.54) —305.71 1598 . 83 S(@sy) | —3559.59| —3096.87 462.72 
S(Hy) 1635.25] —1633.74| —3268.99 S((y) 80.28) —1611.21; —1691.49 

990.28} 2880.99 —645.60} 1362.21 2007 . 81 








The sums of squares and products of the seasonal deviations of the tem- 
perature coefficients are: 













Oo. 4 02 4s o 9s 
60 313126.8 —209260.8 —82955.1 —23973.3 26458 .9 4804.4 
6 — 209260 .8 346463 .7 51448 .4 81359.4 6196.9 —8256.4 
62 —82955.1 51448 .4 209450.0 —98929.8 1136.6 —24013.9 
63 — 23973 .3 81359.4 —98929.8 272087.4 —152747.6 28170.2 
Os 26458 .9 6196.9 1136.6 --152747.6 245352.8  —44700.7 
4804.4 —8256.4 —24013.9 28170.2 —44700.7 115230.8 








and the corresponding values required in the estimation of the regression of 
the yield residuals on the temperature residuals are: 






4 6; 62 4; % 6s 
268968.0 —186357.8 —30584.6  —42734.2 15361.4 —48338.0 
— 186357 .8 272188 .8 1086.9 99731.5 —15544.2 12867 .4 
—30584.6 1086.9 116502.9 —53235.8 —14380.3 1118.8 





—42734.2 99731.5 —53235.8 237056.0 —118395.6 15717.4 
15361.4 —15544.2 -—14380.3 —118395.6 161335.8 36444.0 
— 4838.0 12867 .4 1118.8 15717.4 36444.0 92937 .4 







F2FFPFP oF 







resulting in the values of + listed below: 
Yield from stubble land: 






yo =  .0039530 ys = — .0214618 
v1 = —.0187830 va = — .0208558 
v2 = —.0183526 ys =  .0078885 
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Yield from summerfallowed land: 


yo = .0129911 ¥3 = — .0338827 

m1 = —.0154951 vs = —.0482360 

2 = —.0123744 ~ = .0422734 
Difference in yield from summerfallowed and stubble land: 

vo = .0090381 vs = —.0124289 

v1 = .0032879 v1 = — .0273802 

v2 = .0059782 vs =  .0343849 


The limited extent to which the residual variance of either the observed 
yields or the difference in yield between the stubble and fallow plots can be 
represented by means of these regression coefficients is shown in Table XXVI. 
However, there is a significant correlation between temperature conditions 
and the difference in yield of the Marquis plot from the variety series and 
the fallow plot of the present series. It is possible, therefore, that the relations 
illustrated in Figs. 5 and 6, are not wholly chance effects. 


TABLE XXVI 
DISTRIBUTION OF SEASONAL VARIANCE OF YIELD FROM SUMMERFALLOWED AND STUBBLE LAND 


Summer- Difference 
Variance Stubble fallow in 
yields yields yield 
Sum of squares: 
Due to rainfall regression (6 degrees of freedom) 1374.22 1692.98 127.71 
Due to temperature regression (6 degrees of 249.74 461.51 129.14 
freedom) 
Deviations from regression (30 degrees of | 2081.53 2036.78 1566.65 
freedom) 
Total (42 degrees of freedom) 3705 .49 4191.27 1694.36 
Mean square: 
Due to rainfall regression 229.04 282.16 21.28 
Due to temperature regression 41.62 76.92 21.52 
Deviations from regression 69.38 67.89 47.92 


Discussion of Results 


As the significance of the estimated temperature effects is not definitely 
established, and as the coefficients y are, owing to the high residual variance, 
subject to relatively large errors, an extensive discussion of these results 
would be unprofitable. 


The extent to which the observed yields are affected by temperature 
conditions other than those associated with rainfall or the absence of rainfall 
may be judged from Tables XVI, XXI and XXVI._ It has been suggested 
by Blair (5) that a considerable part of the apparent effect of either precip- 
itation or temperature upon yield is really due to the accompanying effect 
of the other. Contrasting the numerical values of a and y, however, and 
remembering that the standard deviation of the temperature coefficients 
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6 is from 110 to 150 times that of the corresponding rainfall coefficients p, it 
is seen that the total effect on yield of temperature variations (as portrayed 
by the available measurements) is to be judged secondary to that of fluctua- 
tions in rainfall. The value of R, the coefficient of multiple correlation 
between yield, and the rain- 
fall and temperature co- 
efficients, is 0.58 for the 
Marquis plots of the varietal 
test series and 0.66 and 
0.72 for the stubble and 
fallow plots of the grain 
rotation respectively. 

The curves of Figs. 5 and 
6 agree in suggesting that 
above-average temperatures 
are beneficial at the com- 
mencement of the growing 
season, detrimental during ict” eee Mais 
midsummer, and again de- _ Fic. 5. Average effect of temperature (eliminating in- 
sirable during the period of a — on yield of Marquis wheat 
maturation; these various 


WERAGE EFFECT OF LACK ADDITIONAL IF MEAN DALY TEMA (us PER Acre) 


phases succeeding each other at 
somewhat shorter intervals of 
time on the less fertile soil of 
the grain rotation (Fig. 6) than 
on the more fertile soil of the 
variety plots (Fig. 5). Such 
effects are in reasonable accord- 
ance with expectation. The 
reduction in yield accompany- 
ing high temperatures still later 
in the season might be due 
partly to the acceleration of 
respiration in the ripe grain, 
with increased losses of carbo- 
ilies °Tensiatiak: < Deetal hydrate material, and partly to 


Fic. 6. Average effect of temperature aon reduction of the moisture con- 
influence of associated rainfall) on yield of wheat from i i 
stubble plot (solid line) and summerfallow plot (broken eae of the grain * 5 he time of 
line) of grain rotation. threshing and weighing. 

Reference has already been made to the greater moisture-holding capacity 
of the variety test plot soils, and the higher specific heat of moist soil will of 
course lead to slower warming in the spring. Above-average air temperatures 
before seeding may be desirable therefore in order to raise the temperature 
of these relatively moist soils to levels more favorable to the germination 
of the seed and the early growth of the plant. This consideration will apply 
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with diminished force in the case of the grain-rotation plots; and on the 
other hand, owing to its lower fibre content, their soil is liable to blow in 
the spring. The annual reports of the superintendents of the experimental 
stations contain frequent references to such occurrences, which must be 
aggravated by the drying out of the surface layer, a process which will un- 
doubtedly proceed more rapidly at higher temperatures. Even in these 
circumstances, however, additional warmth during the period of germination 
and early growth appears advantageous. 


As the temperature conditions in the spring vary more than do those of 
the succeeding summer months (Fig. 7) it seems likely that some of the most 


Fic. 7. Upper portion: Average mean monthly temperature (degrees Fahrenheit) recorded at 
meteorological stations in central and southern Alberta and Saskatchewan. Lower portion: Annual 
variation (standard deviation in degrees Fahrenheit) in mean monthly temperature. 


pronounced effects of temperature upon yield may be exerted at this time. 
At all stations the month of July is characterized by the maximum mean 
temperature. That both the average values and the annual variation (stan- 
dard deviation) of the corresponding temperatures are in general higher at 
the southern than at the northern, and at the Saskatchewan than at the 
Alberta stations, is also to be observed from Fig. 7. 


7. General Conclusions 


Cole and Matthews (9), studying the use of water by spring wheat on 
the Great Plains, found a correlation of 0.76 between the total quantity of 
water removed from the soil and the yield of grain. The supply of moisture 
upon which the crop may draw is derived from the amount available in the 
soil at the time growth begins, supplemented by the rain falling during the 
growing period. The former they found quite inadequate to carry the crop 
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to maturity, and under western conditions both combined are seldom in 
excess of the optimum requirements of the crop. In these circumstances the 
yield secured might certainly be expected to be in some degree dependent on 
summer rainfall; and the results of the present investigation, as well as the 
significant correlation between yield and seasonal precipitation found by the 
above authors, indicate that this is the case. The nature of the relation 
seems to be influenced by the soil conditions, and even under similar soil 
conditions is doubtless affected by the amount of moisture originally present. 
If this is so, the various curves of the preceding sections may be expected to 
indicate the average relation in the case of the soils, stations and seasons 
studied. 


In the same way the practical importance of any “critical periods” in plant 
development will depend not on the absolute moisture and temperature 
requirements of the plant, but on the deviation of these requirements from 
the environmental conditions experienced at that time. Disparities between 
the rainfall response curves of different plots, as in Figs. 1 and 3, are not 
therefore, necessarily incompatible with the existence of relatively ‘‘critical” 
periods of plant growth in the sense referred to in Section 1; but they do 
suggest that the extent to which yields are actually limited by the rainfall 
during such periods may be subject to considerable modification by soil 
factors. Temperature conditions affect the plant more directly, and this 
circumstance perhaps accounts for the greater degree of similarity between 
the effects illustrated in Figs. 5 and 6. 


Whatever may be the expectations on physiological grounds, the results 
of this study indicate that in the region to which they apply, the influence 
of weather conditions on wheat yield is not largely exerted in a few relatively 
short periods of time, but extends at least throughout the growing season. 
This conclusion is of importance in the formulation of prediction equations 
designed to forecast commercial yields from the known weather conditions. 
The limited effect upon yield of temperature conditions other than those 
associated with rainfall also has an obvious bearing in this connection, though 
in view of the tendency of cultivators to commence seeding operations as 
early as possible each spring, the observation that low temperatures during 
the initial stages of growth are associated with reduced yields may not be 
without practical significance. 


It might be thought that the moderate degree of correlation between weather 
conditions and the plot yields studied, renders unlikely the successful fore- 
casting of commercial production; but this is not necessarily so. The crop 
returns considered, being obtained from a different small area of land at 
each station each year, were undoubtedly affected by soil differences. If, 
however, the total production of a district is considered, soil differences 
above and below the average may be expected to a considerable extent to 
offset each other, resulting in a greater relative importance of crop variation 
due to seasonal weather conditions. 
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The results of Section 6, indicating that the yield secured is influenced less 
by temperature conditions than by rainfall during the growing season, are of 
interest in view of Chilcott’s conclusion (8) that seasonal precipitation is 
not the dominant factor determining crop yields in the Great Plains area. 
Whilst the “inhibitory factors” such as hot winds, diseases, and insects, to 
which this author rightly directs attention are undoubtedly of great im- 
portance in some seasons, it seems probable that the statistical method em- 
ployed by him, namely the calculation of the ratio of the yield of grain to the 
seasonal precipitation, as well as the fact that the total precipitation only, 
regardless of its distribution, was considered, has led to some under-estimation 
of the actual association between yield and rainfall. It is obvious, however, 
that the yields secured are markedly influenced by factors other than rainfall 
and temperature during the growing season, and the amount of precipitation 
during the autumn, winter and spring months prior to sowing is being studied 
in this connection. 
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SOME FACTORS IN THE DIFFERENT CHROMOSOME 
SETS OF COMMON WHEAT! 


By W. P. THompson?, T. J. ARNASON® AND R. M. Love? 


Abstract 


In crosses between 21-chromosome (vulgare) and 14-chromosome (emmer) 
wheats many hybrid lines which had only 14 pairs of chromosomes were iden- 
tified and were then studied from the standpoint of their emmer and vulgare 
characters. Any vulgare character which appears in such lines must be deter- 
mined by genes in the primary set of 14 chromosomes which mate in F; with 
those of emmer. Vulgare characters not appearing in such lines are determined 
by genes in the secondary set of 7 ak tenia unmated in F;. In this 
way it has been determined that 7 of the 19 characters studied in durum 
crosses and 9 of the 21 studied in dicoccum crosses are determined by genes in 
the primary set of 14 (for list see Table II). For several other vulgare characters 
it has been shown that necessary genes are in both the primary 14 and the 
secondary 7 (Table III). There is little to indicate that genes determining im- 
portant species-distinguishing characters are concentrated in the secondary set 
of 7, as is implied in the theory that vulgare wheat has resulted from the crossing 
of an emmer with Aegilops and that the set of 7 is derived from Aegilops and 
homologous with a set in that genus. Evidence in regard to the mode of in- 
heritance of certain characters is given. 


Introduction 


All of the stable lines which result from crossing common wheats with emmer 
types have the parental chromosome numbers (” = 21 or 14). The offspring 
of hybrids with intermediate numbers tend to revert to those of the parents 
(1, 2, 4) and eventually the reversion becomes complete in all lines. The 
inheritance of any character of common wheat will therefore depend, among 
other things, on whether it is determined by genes in the two sets of 7 which 
mate with those of emmer (hereinafter called the primary 14 or sets A and B), 


or by genes in the 7 which remain unmated in F, (hereinafter called secondary 
7 or set C). 


If a vulgare character appears in a hybrid segregate which has only 14 pairs 
of chromosomes, all the genes necessary for its production must be located 
in the primary 14 (Aand B). A hybrid line with 14 pairs only cannot show a 
common wheat character whose genes are in the secondary 7 (C); nor can a 
line with 21 pairs show an emmer character whose vulgare counterpart is deter- 
mined by genes in the secondary 7, unless the emmer character is epistatic to 
that of vulgare. The substitution of a particular emmer character for a 
vulgare one in a hybrid of the general vulgare type will therefore depend on 
whether the vulgare genes involved are in the primary 14 or in the secondary 7. 


The purpose of this investigation has been to determine in which set of 
chromosomes (A and B on the one hand or C on the other) the genes for many 
characters are located. This has been done by determining whether the 
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characters in question occur in 14-chromosome* hybrid lines or whether their 
emmer counterparts occur in 21-chromosome hybrid lines. The results 
obtained with the 14-chromosome lines are given in the following pages. 
The results with 21-chromosome lines which are more difficult to obtain and 
study in the necessary large numbers, are not as yet sufficiently extensive 
to report. Supporting evidence has been obtained by determining whether 
vulgare characters which are not found in 14-chromosome lines occur in plants 
which have 15 (14 bivalents and one of the C chromosomes). 

This problem is given additional significance by the widely held view, 
originally suggested by Percival (3) on taxonomic grounds, that common 
wheat originated through the crossing of an emmer type with Aegilops, and 
that the characters which distinguish common wheat from emmers were 
contributed by that genus. This conception has been given a precise cyto- 
logical significance by recent observations, which are taken.to indicate that 
the secondary set of 7 chromosomes (C) present in vulgare wheat has come 
from Aegilops and therefore should carry the distinctively vulgare genes. 
The cytological observations on which this conclusion is based have been 
reviewed by Thompson (5). 

Before cytological work had been done on Aegilops or its hybrids with 
wheat, Sax (4) had reached the conclusion from his studies on vulgare-emmer 
hybrids that the distinguishing characters of common wheat are deter- 
mined by genes in the C chromosomes. He wrote ‘Apparently the 7 addi- 
tional chromosomes of vulgare varieties determine the distinguishing charac- 
ters of common wheats due to the reduplication of hereditary factors, or to 
specific factors in these chromosomes .... The cytological and genetic 
data would indicate that many of the desirable characteristic properties of 
the emmer and vulgare wheats cannot be combined in a homozygous condition.” 

In the present work no general attempt has been made to analyze the 
genetics of the different characters, to determine how many genes are in- 
volved or their mode of interaction. In most cases the results reported 
show only whether all the genes which are necessary for the production of 
the character in question are in the primary 14 chromosomes. In some cases 
it is shown that necessary genes are in both the primary 14 and the secondary 7. 


Materials and Methods 


In extensive earlier work on the cytology of wheat species crosses a large 
number of hybrid segregates had been identified as having only 14 bivalent 
chromosomes at the first meiotic division. Others had 14 bivalents and 
1 univalent. Some of them were F; plants resulting from selfing F; and F, 
of vulgare Xemmer. Many were the results of backcrossing F,; to the emmer 
parent. Offspring of them were grown and studied both genetically and 
cytologically. 


* For the sake of clearness and ease in deciding which set is involved in each case haploid 
numbers are used throughout. ‘A 14-chromosome plant" means one which shows only 14 bivalent 
chromosomes at the first meiotic division and therefore has 2 of the primary sets of 14 or 28 somatic 
chromosomes. ‘‘A 15-chromosome plant"’ means one with 14 bivalents and 1 univalent, the latter 
belonging to the secondary 7. 
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With respect to some of the distinguishing parental characters exact classi- 
fication was difficult. Various intermediate conditions appeared. Several 
characters described as units in taxonomic work broke up into sub-characters 
in the hybrids. Entirely new conditions not found in either parent occurred. 
The exact condition was recorded as accurately as possible for each plant. 

The chromosome number and mating at the first division were determined 
for several plants in each line by means of smear preparations. This not 
only checked the parental count but made it possible to state explicitly that a 
particular vulgare character had been found in a particular 14— (or 14+1) 
chromosome plant. Since all plants had 14 bivalent chromosomes plus 
0 or 1 univalent it was only necessary to determine whether or not a univalent 
was present at the meiotic division. At least 20 pollen mother cells of each 
plant were examined. 

The amount of material studied is shown in Table I. The total for 14- 
chromosome lines is 83 families and 757 individuals, and for 15-chromosome 
lines 34 families and 274 individuals. All lines are descended from different 
F, or backcross plants. The numbers should be large enough to make the 
results reliable. In addition six families including 38 individuals were examined 
from 14-chromosome backcrosses of vulgare X persicum F, to persicum. 


TABLE I 


NUMBER OF LINES AND PLANTS STUDIED (EACH LINE FROM A 
DIFFERENT F; OR BACKCROSS PLANT) 


F, by selfing Fi, F2, and Fy Backcrosses 








Chromosomes*| Durum crosses | Dicoccum crosses With durum | With dicoccum 


Lines Plants 


Lines | Plants Lines Plants Lines | Plants 


iciiideenniniesibaak a 
| 





. 141 6 72 14 149 30 253 33 283 


1411 + 11 2 m4 8 77 13 99 11 72 


| 
*] = univalent; II = bivalent. 
The parental varieties used were as follows: vulgare—Marquis; durum— 
lumillo; dicoccum—Vernal; persicum—Black Persian. 


The Characters Studied 


There are very few characters which distinguish all varieties of vulgare 
from all varieties of durum or dicoccum. But there are several characters 
which are confined to one species although they are not found in all varieties 
of that species. And there are several characters which do not differ in 
kind as between two species but which are usually developed to a greater 
degree in one species than in the other. Moreover there are numerous 
characters which distinguish the great majority of individuals and forms of 
one species from those of the others and which may therefore be regarded as 
typical of the species in question. Watkins (7) has discussed at length the 
differences between and interrelationships of the species of Triticum. 
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The characters studied in the present work are in some degree species- 
distinguishing, although, as indicated above, few of them distinguish all forms 
of vulgare from all forms of durum or dicoccum. For example there is some 
overlapping in regard to the diameter of the stem but vulgare forms as a rule 
have thicker stems than durum or dicoccum. The absence of beards is a 
feature of many vulgare types (including the one used in this study), but not 
of all; beards are never absent from durum or dicoccum types. The des- 
criptions which follow apply to the particular forms used in this study. 


1. Diameter of stem—measured at 2 cm. below the spike; emmer (durum 
and dicoccum; hereafter designated “‘e’’) 1.1 to 1.5 mm.; vulgare (hereafter 
designated “‘v’’) 1.6 to 2.2 mm. 

. Cavity of stem—observed at 2 cm. below the spike; e solid, v hollow with 
thin walls. 

. Leaf hairs—of young leaves; durum none or very few, short and scattered; 
dicoccum very numerous, long, evenly distributed; vulgare long on tops 
of ridges, short on sides, sparse. 

. Collar—a collar-like structure surrounding the stem at the base of the 
spike; in e a complete circle; in v open at one side. 

. Density of spike—calculated by dividing the length of the spike in milli- 
metres by the number of spikelets. 

The value so obtained is slightly greater than the average internode 
length. e 3 to 4; v 4.4 to 5.6. 

. Form of spike—width of one-ranked side divided by width of two-ranked 
side. In dicoccum 0.6 to 0.8; durum 0.7 to 1.0; vulgare 1.3 to 1.7. 

. Beards—e 8 to 12 cm.; v short tip awns. 

. Glume shape—empty glume, e, long, tapering to both ends; v, short, 
wide, blunt. 

. Glume cross section—e, V-shaped; v, U-shaped. 

. Glume tip—durum long, sharp; dicoccum short, pointed; vulgare short, 
blunt. 

. Glume shoulder—the portion at the side of the tip; e narrow, sloping; 
v wide, more horizontal. 

. Glume adherence—dicoccum tight, difficult to remove; durum and vulgare 
loose, easily removed. 

. Keel—ridge on back of empty glume; e prominent, sharp, extending full 
length of glume; v prominent above, weak below. 

. Keel-teeth—durum fine, hair-like, numerous from tip to near base; 
dicoccum few, short, blunt, near the tip; vulgare short, blunt, sparse, on 
upper half. f 

. Shape of segment of rachis between successive spikelets—e slender, top 
twice as wide as base; v stout, top little wider than base, margin curved. 

. Rachis hairs—durum, a marginal tuft at the top of each segment tapering 
downwards; dicoccum, a central tuft between the glume bases; vulgare, 
a uniform fringe along the sides and across the top. 
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17. Rachis articulation—dicoccum fragile, the spike readily breaking into 
segments; durum and vulgare tough, remaining intact during threshing. 

18. Seed shape—e long, narrow, tapering to ends; v short, stout, blunt. 

19. Seed cross section—e somewhat triangular with ridge on back; shallow 
furrow; v rounded back and sides, deep furrow. 

20. Seed end—embryo end; durum one projection; dicoccum a longitudinal 
furrow; vulgare two blunt projections. 

21. Seed hairs—on stigma end; durum few, short, dicoccum many short; 
vulgare many long. 


Experimental Results 
A. Vulgare CHARACTERS IN SEGREGATES WITH 14 PAIRS OF CHROMOSOMES 


The characters of the 21-chromosome species (vulgare) which have been 
found in 14-chromosome hybrid lines are shown in Table II. In regard to 
leaf hairs only a single plant with 14 chromosomes in durum crosses and back- 
crosses was recorded 


TABLE II 

as having the vulgare 

oe Vulgare CHARACTERS FOUND IN 14-CHROMOSOME LINES AND 
condition. The same NUMBER OF LINES IN WHICH THEY OCCURRED 
is true with respect 
to the cross section of F, of crosses with |F;, of backcrosses with 
the glume. Since the Character pee En nana ee 

ape “ durum | dicoccum | durum | dicoccum 
Ses «GE Ue ain eee pine nen 


errors in the records 


; Stem diameter 2 2 6 1 

cannot be entirely Leaf hairs 1 6 0 9 

. Beards 5 11 9 8 

ruled out, it may be Glume section 1 2 0 2 

questionable whether Glume ti 3 4 9 7 

the vulgare condition — adherence ; : ; es 

of these two charac- Keel teeth 0 5 0 7 

Seed hairs 2 0 4 0 

ters has actually been oka ; 3 2 5 
found in 14-chromo- 

6 30 33 





some durum segre- Lines studied 
gates. There is no 
doubt that they occur in 14-chromosome dicoccum segregates. In regard 
to all the other characters more than one family and several plants in a 
family showed the vulgare condition. 

It will be observed that seven (or possibly only five) of the 19 vulgare 
characters have been found in 14-chromosome vulgare-durum hybrids, and 
nine of the 21 vulgare characters in 14-chromosome vulgare-dicoccum hybrids. 
The total number of characters studied is different in the two crosses because 
in regard to two of them, glume adherence and rachis articulation, vulgare 
is similar to durum and quite different from dicoccum. The reason for the 
greater ease with which vulgare characters can be transferred to 14-chromo- 
some plants in vulgare-dicoccum than in vulgare-durum crosses will be dis- 
cussed later. But the evidence from both crosses shows that by no means all 
vulgare characters are determined by the secondary seven chromosomes. 
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Furthermore the vulgare characters which have been transferred include 
some important species-distinguishing characters such as the condition of 
the leaf hairs, cross section of the glume, diameter of stem, and the keel. 
But they also include some which do not distinguish all vulgare varieties 
from all emmer varieties, such as absence of beards, or blunt glume tip. 


It is of interest to note that with respect to ten of the 21 characters studied 
in vulgare-dicoccum, the F, is intermediate; with respect to nine others the 
vulgare condition is partially or completely dominant; with respect to the 
remaining two it is recessive. Of the nine characters transferred four are 
intermediate in F;, four are partially or completely dominant and one recessive. 
In vulgare-durum four of the five characters transferred are dominant in some 
degree, while one is intermediate in F,. 


In view of the large amount of sterility and chromosome irregularity in F, 
it is not considered advisable to draw any conclusions with respect to the 
mode of inheritance from the proportion of families in which a character 
appeared. 


Some characters which have not been found in the fully developed vulgare 
condition in 14-chromosome plants, have been found in a near-vulgare con- 
dition. These are shown in Table III. Presumably important genes 
influencing them are located in the primary 14 chromosomes but 

TABLE III for their complete ex- 
pression genes located 
CHARACTERS FOUND IN NEAR-vulgare CONDITION IN 


14-CHROMOSOME LINES AND NUMBER OF FAMILIES in the secondary 7 are 
IN WHICH THEY OCCURRED necessary in addition. 


It may be, however, 

F, of crosses with |F,ofbackcrosseswith that all the genes 
aoeasem | dicoccum  mecessary for the ex- 
pression of some of 
Spike density them are in the prim- 
Glume shape ary 14 and _ that 
— — emmer genes which 
Rachis segment inhibit their full de- 
— velopment have not 
been entirely elimin- 

Lines studied ated from any of the 
plants studied. 


Character 


durum | dicoccum | durum | dicoccum 


Anne Sw I 


is 
aos 


Most of those characters which have been found in the near-vulgare con- 
dition in 14-chromosome plants have been found in the completely vulgare 
condition in plants which have a single C chromosome in addition to the 
14 pairs. This is shown in Table IV. It will be observed that the characters 
listed in this table are almost the same as those in Table III. Presumably 
the extra chromosome from the secondary 7 carries the gene which is necessary, 
in addition to those in the primary 14, for the full expression of the vulgare 
condition. (The same vulgare characters which are found in plants which 
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have no univalent chromosomes have, of course, been found in those plants 
which have one univalent in addition, but they have not been shown in 
the table.) 


New characters TABLE IV 
: CHARACTERS FOUND IN vulgare CONDITION IN PLANTS WITH 
found in Table IV 14 BIVALENT AND 1 UNIVALENT CHROMOSOMES 
may of course be due 


entirely togenes in the f F, of crosses with | F,of backcrosses with 

additional chromo- Character pace aera ene eee 
: Uru 1cO0 

somes. Since there are | ee) ee 


seven extra chromo- 








: Spike form 1 - 2 1 

somes and since the Spike density = 1 2 2 
number of families Coles : 1 : 3 a 
with anextrachromo- ¢ — eats 3 7 3 
some was small, it is Rocks hairs 1 i 1 - 
achis segment =- = - 

probable that further Rachis auieaiaahiis = 2 ~ - 
workon such lineswill Seed shape - 1 1 - 
add more characters | ines studied 2 8 13 11 


to this category. It 

should add all those 

vulgare characters which are determined by a single kind of gene in the 
secondary 7, except those for which the gene must be represented twice as 
in pure vulgare. 


In addition to the durum and dicoccum crosses which have just been des- 
cribed, the progenies of six 14-chromosome plants resulting from backcrossing 
vulgare X persicum F, to persicum were studied. T. persicum differs from the 
variety of vulgare which was used, in having black pubescent glumes, lighter 
seed color and rugose caryopsis, in addition to many of the characters which 
distinguish durum and dicoccum from vulgare. The vulgare seed and glume 
color and absence of pubescence were found among the six families. A near- 
vulgare condition of the surface of the caryopsis was also found. 


B. ADDITIONAL RESULTS AND THE MODE OF INHERITANCE 
OF CERTAIN CHARACTERS 


In regard to certain characters additional evidence was obtained which 
cannot be given in tabular form. Some of this evidence bears on the mode 
of inheritance of certain characters. Although a complete analysis of the 
mode of inheritance was not attempted, and would be very difficult owing 
to sterility and chromosome irregularity in F, certain conclusions may be 
pointed out. - 


1. Cavity of Stem 


No 14-chromosome hybrid has been obtained with the hollow, thin-walled 
stem characteristic of vulgare. Many have been obtained with small cavities 
but they all had thick walls. It would appear, therefore, that a factor for 
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hollow stems is present in the primary 14 chromosomes of vulgare and that 
in order to produce the large cavity with thin walls at least one additional 
factor located in the secondary 7 must be present. 


2. Leaf Hairs 

T. vulgare has long hairs on the crests of the ridges of the young leaves and 
short ones on the sides of the ridges; in the variety used, both kinds are 
sparse. 7. dicoccum has numerous long hairs uniformly distributed, and 
durum has very few short ones or none at all. The vulgare condition was 
transferred without difficulty to 14-chromosome dicoccum hybrids and is 
therefore determined by genes in the primary 14 chromosomes. Consequently 
it should be transferred to 14-chromosome durum hybrids. But this occurred 
in only one case. It may be that durum genes inhibit the vulgare condition 
and only one plant was obtained which possessed the necessary vulgare 
chromosomes and not the durum chromosomes which carry the inhibiting 
genes. It is of interest in this connection that in vulgare X durum F, the 
durum condition is nearly dominant, the hairs being few and short, and that 
in vulgare Xdicoccum the hairs are numerous as in dicoccum, mostly rather 
long, with some still longer on the crests of the ridges as in vulgare. 

Possibly the fact that the vulgare condition of the leaf hairs is readily trans- 
ferred to dicoccum but not to durum hybrids may be due to the absence of 
genes for long hairs in durum and the primary 14 of vulgare and their presence 
in dicoccum and the secondary 7 of vulgare. In that case there would be no genes 
for long hairs in 14-chromosome durum hybrids; the long hairs of dicoccum 
hybrids would be due to genes from dicoccum, and the difference in length 
between the top and the side of ridges would be due to genes in the primary 
14 of vulgare. 

According to Percival (3) the arrangement, number, and kind of hairs in 
the young leaves of wheat are among the most constant species-distinguishing 
characters. In the hybrids these are very variable. It appears that different 
genes affect the characters since various combinations of the numbers, lengths, 
and arrangements of the parental species were observed. But all the genes 
necessary to produce the vulgare condition are apparently in the primary 
14 chromosomes. 


3. Density of Spike 

In several 14-chromosome lines the spikes were almost, but not quite, as 
lax as in pure vulgare. Apparently there is a spike-lengthening gene (or 
genes) in the primary 14 chromosomes but an additional gene (or genes) is 
required from the C chromosomes in order to praduce the full vulgare laxity. 
There is a similar situation with respect to several other characters. It 
may be mentioned also that some of the backcrosses to vulgare (not reported 
in this paper) had even laxer spikes than vulgare. This indicates either that 
durum, though very dense, contains spike-lengthening factors not present 
in vulgare, or that vulgare contains spike-shortening factors which have been 
eliminated in the hybrids. 
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4. Glume Tip 


The short blunt tip on the empty glume which is characteristic of vulgare, 
appeared in many 14-chromosome hybrids. In addition several lines had a 
much longer, sharper tip even than durum; it was in fact a short awn. One 
of the parents must carry, in the primary 14 chromosomes, genes for a short 
awn on the empty glume and also genes which suppress the awn. In some 
hybrids the awn-producing gene is present without the suppressor. 


5. Keel 


The weakly developed keel characteristic of vulgare was present in many 
14-chromosome vulgare Xdicoccum hybrids. Some had an even weaker keel 
than vulgare. But the vulgare condition was not found in durum hybrids. 
The nearest approach to it was a low but distinct keel extending the full 
length of the glume. It appears that factors are present in durum and not 
in dicoccum which promote keel development and that these have not been 
eliminated in any of the durum hybrids studied. There is a similar difference 
between durum and dicoccum crosses in the difficulty of transferring the 
vulgare condition of keel teeth and of leaf hairs. With respect to the latter 
two characters, however, dicoccum resembles vulgare much more closely than 
does durum. The fact that three vulgare characters can be readily trans- 
ferred in dicoccum but not in durum crosses therefore has no general sig- 
nificance. One character, seed hairs, was transferred in durum crosses but 
not in dicoccum crosses. 


In a cross between vulgare and turgidum Watkins (6) found plants without 
keels or with keels even less developed than in vulgare. His interpretation is 
that turgidum has the factors KK for the prominent keel and that vulgare 
has kk for the absence of keel in the primary 14 and K,K, for weak keel in 
the secondary 7. Keel-less extracted hybrids are then kk; 21-chromosome 
plants of the genotype kkK,K;, have the vulgare keel, and those of the genotype 
KKK,K, have an even more strongly developed keel than turgidum. The 
facts mentioned with respect to the vulgare-dicoccum cross support this inter- 
pretation, and the inheritance of other characters, such as glume tip, may 
be interpreted in a similar way. 


6. Shape of Rachis Segment 


In durum and dicoccum the segment of the rachis between successive spike- 
lets is narrow, slender and somewhat triangular, the top being nearly twice 
as wide as the base; in vulgare it is thick and stout, with a characteristic mar- 
ginal curve and with the top little wider than the base. Such a variety of 
conditions was found in the hybrids as to indicate that several gene differences 
are involved. It may be that all the genes necessary to produce the vulgare 
condition are in the primary 14 chromosomes but not all have been obtained 
in one plant. Since certain elements of the vulgare condition are readily ob- 
tained in 14-chromosome plants it appears that some of the necessary genes 
are in the primary 14. 
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7. Rachis Hairs 


Durum has a marginal tuft of hairs at the top of each segment, decreasing 
in length in a downward direction; dicoccum has a central tuft between the 
bases of the glumes; vulgare has a uniform fringe of hairs along the full length 
of the margin and across the top. Various combinations of these conditions 
were noted in the hybrids. For example some had tufts in both the durum 
and dicoccum positions; some had a near-vulgare fringe with a dicoccum tuft; 
some had no hairs at all. The character is apparently influenced by several 
genes and some of these are in the secondary 7 chromosomes. 


Discussion 


It has been pointed out that if a vulgare character occurs in a hybrid segregate 
which has only 14 bivalent chromosomes, all the genes necessary for its 
expression must be in the primary 14 chromosomes which mate with those 
of emmer. Five (or possible seven) of the characters studied fall in this 
category in durum crosses, and nine in dicoccum crosses. If such a character 
appears in any considerable proportion of 14-chromosome plants, genes in 
only one or a few chromosomes must be sufficient for its expression. Several 
characters appear to be of this kind. 


If a vulgare character is partly but never fully expressed in 14-chromosome 
segregates two explanations are possible: (1) Some of the vulgare genes in- 
fluencing the character are in the primary 14 and some in the secondary 7 
chromosomes, or (i) all are in the primary 14 but the presence of emmer 
genes, whether allelomorphic or not, inhibits the full development of the 
character. The latter intrepretation assumes that in all the plants examined 
the emmer chromosomes carrying the modifying genes are still present. 
The question then is whether the number of plants examined is large enough to 
include some which have only vulgare chromosomes carrying the genes in- 
fluencing the character. The number examined in the present study appears 
to be large enough, unless the character is influenced by genes in several 
chromosomes. Therefore the first interpretation appears in general to be 
correct—the vulgare characters in question require for their full development, 
genes in both the primary 14 and the secondary 7. 


This conclusion is supported by the fact that nearly all the characters in 
this category appear in the full vulgare condition among plants which have 
an extra C chromosome. 


It is clear from the results that the primary 14 chromosomes carry all the 
genes necessary for the production of several vulgare characters and part of 
those necessary for the production of others. There is little to indicate that 
the important species-distinguishing genes are concentrated in one set of 
chromosomes rather than the other. This conclusion is not in agreement 
with that which was reached when these interspecific crosses were first 
studied (4). It was thought that the specifically vulgare genes are concen- 
trated in the secondary 7 chromosomes and that it would therefore be im- 
possible to obtain 21-chromosome plants with emmer characters. 
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Some of the conclusions reached in this work are supported by the pre- 
liminary study of the extracted hybrids which have 21 chromosomes. These 
plants show the emmer counterparts of some of the vulgare characters which 
have been found in 14-chromosome segregates. The presence of the secondary 
7 chromosomes has not caused the development of the vulgare characters in 
question. 


According to a theory which has received wide recognition, vulgare wheat 
originated through the crossing of a member of the genus Aegilops with an 
emmer wheat. Aegilops contributed the secondary set of 7 chromosomes 
which are still homologous with a set in Aegilops. It also contributed (in 
those chromosomes) the genes for the distinctively vulgare characters. If 
that theory is correct the genes must have been extensively shifted into the 
other sets by translocations. But there is no evidence of such translocations 
in the numerous reports of the behavior of the chromosomes at meiosis in 
hybrids. 
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FRACTIONATION OF GLUTEN DISPERSED IN SODIUM 
SALICYLATE SOLUTION! 


By A. G. McCALta? AnD R. C. Rose? 


Abstract 


The solubility of gluten in sodium salicylate solution increases regularly with 
increasing concentrations of the solvent. When gluten is completely dispersed 
and then precipitated by the addition of magnesium sulphate, the relation 
between the amount of gluten nitrogen precipitated and the concentration of 
magnesium sulphate is also regular. Successive fractions of the precipitated 
pao protein (except the last 10-15%) contained progressively more amide and 
less arginine nitrogen. None of the fractions were similar to gluten, but when 
they were redispersed, combined and oe as a whole, a gluten was 
obtained. The most soluble 10-15% of the gluten protein is probably distinct, 
but the remainder is a single protein complex which can be progressively 
fractionated. 


Introduction 


Recently the view that gluten is composed of two distinct proteins, glutenin 
and gliadin, has been seriously challenged. ‘These proteins, in common with 
most others, are defined chiefly on the basis of their solubilities, and there is 
considerable evidence that such definitions are only arbitrary. The fraction 
of gluten which is soluble in 70% alcohol has been called gliadin, and the 
remainder, which is soluble in dilute acid or alkali, glutenin. 

Gortner, Hoffman and Sinclair (8), however, have shown that practically 
any amount from 15 to 60% of the protein in flour can be extracted by choosing 
suitable salt solutions. Doubtless for each salt solution there is a concen- 
tration which extracts the maximum amount of protein, and a further increase 
in concentration either does not affect, or else decreases, the amount extracted. 
Recently Cook and Rose (4) have shown that a sodium salicylate solution of 
sufficient concentration will disperse practically all of the protein in gluten. 
It appears that solubility in the salt solutions used by Gortner, Hoffman and 
Sinclair, or even in sodium salicylate solution, might be used in defining 
gluten proteins, with as much reason as solubility in alcohol, acid and alkali. 

There is also other evidence that the definitions of glutenin and gliadin are 
only arbitrary. Thus Herd (9) has shown that the amount of gliadin obtained 
from samples of flour depends on the method of extraction, and believes that 
so-called ‘‘gliadin”’ is not an individual protein. Cook (2) fractionally preci- 
pitated gluten from a dispersion in 30% urea solution and concluded that, 
on the basis of solubility, there was no clear-cut distinction between glutenin 
and gliadin. 

Chemical analyses have shown that so-called “glutenin” and “gliadin” 
differ in their amide and arginine content. Different workers have obtained 
gliadin samples of approximately constant amide and arginine nitrogen content 
and so, from the standpoint of consistent chemical composition, gliadin may 


1 Manuscript received February 13, 1935. 
Contribution from the Department of Field Crops, University of Alberta, with financial 
assistance from the National Research Counci of Canada. Issued as paper No. 68 of the Associate 
ittee on Grain Research. 
2 Research Assistant, Associate Committee on Grain Research, University of Alberta. 
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lay some claim to individuality. This, however, is not true of so-called 
“glutenin”’, as different workers have obtained widely divergent results (3). 
Sandstedt and Blish (14) fractionated the gluten protein by precipitation 
from dispersions in acetic acid, and found that the fractions produced differed 
gradually in chemical composition. They suggested that gluten is composed 
of three main proteins, the third being intermediate between glutenin and 
gliadin in its properties. Cook’s results (2), however, show that the relation 
between solubility and arginine content of fractions was essentially linear. 

In order to fractionate the gluten protein it is necessary to have a suitable 
agent which will disperse the gluten completely. Ideally, the reagent would 
be such that the protein could be recovered without the loss of any of its 
original properties. The known solvents are alkali, acid, urea, and sodium 
salicylate solutions. Sandstedt and Blish (14) showed that dilute acetic 
acid is preferable to dilute alkali, and the recent work of Cook and Rose 
(5, 6, 13) showed that the neutral solvents are superior to either alkali or acid. 

Sodium salicylate was chosen for use in the present study since the use of 
urea involves the laborious procedure of removing the urea nitrogen. 


Materials 


Four flours,were used in this study, namely: a commercial pastry flour, a 
low protein flour milled from Marquis wheat, a medium protein flour milled 
from a Reward wheat and a high protein flour milled from another Marquis 
wheat. These are designated respectively as Flours 1, 2,3 and 4. The last 
three were experimentally milled and free from bleaching agents and chemical 
improvers. All but Flour 1 were of high baking quality. Analytical data 
concerning these flours are given in Table I. The methods used in deter- 
mining these values are given in the following sections. 


TABLE I 
NITROGEN IN FLOUR AND WASHED GLUTEN 


Nitrogen in washed gluten 
Total nitrogen 
in flour Total, as Amide, as Arginine, as 
(dry basis), percentage of total | percentage of total | percentage of total 
0 nitrogen in flour | in washed gluten | in washed gluten 


Washed Gluten Methods 


All glutens were washed from the flour, using the phosphate buffer solution 
of pH 6.8 recommended by Dill and Alsberg (7). The amount of flour weighed 
out depended on the experiment for which the gluten was used, and the rate 
of flow of the washing solution and the time of washing were varied according 
to the sample size. 
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Extraction of Gluten 


Concentrations of sodium salicylate varying from 0 to 10% were employed, 
and the method was essentially that of Cook and Rose (4) in which the gluten 
washed from 5 gm. of flour was shaken at frequent intervals for six hours with 
100 ml. of solvent. Cook and Rose found that increasing the amount of 
solvent, or the time of extraction, had little effect on the amount dispersed, 
and although this was verified with the higher concentrations of sodium sali- 
cylate as used by Cook and Rose, repeated extraction with the lower concen- 
trations did slightly increase the amount extracted. A six hour extraction 
period, however, was used with all solutions. When extraction was carried 
out with water it was first saturated with toluene to limit bacteriological 
action. In the first experiments the extract was filtered through glass wool 
as centrifuging left floating particles, but in subsequent experiments the foam 
produced by shaking was broken by the addition of caprylic alcohol and the 
extract centrifuged. The results obtained by the two methods were slightly 
different as shown later. 


Preparation of Dispersions 


Dispersion was accomplished after the method of Rose and Cook (13) by 
placing the wet gluten in 10% solution salicylate solution and shaking vigor- 
ously at frequent intervals. Sufficient wet gluten was used to give a con- 
centration of approximately 5 mg. of nitrogen per ml. At the end of five 
hours when dispersion was complete the dispersion was placed in the refriger- 
ator at 1° C. and left overnight. The following morning it was centrifuged 
at approximately 2500 r.p.m. for 20 min. and the liquid decanted into a 
volumetric flask of appropriate size, and made up to volume. 


Precipitation of Dispersed Protein 


The protein was fractionally precipitated by the successive addition of 
calculated volumes of a magnesium sulphate solution, saturated at 20°C. An 
attempt was made to prepare a saturated magnesium sulphate solution in 
10% sodium salicylate but this was impossible, owing probably to the form- 
ation of magnesium salicylate. When, however, equal volumes of saturated 
magnesium sulphate and 10% sodium salicylate were mixed, no precipitate 
formed even when the solution was stored at 1° C. for several days. After 
the addition of the salt the dispersion was mixed thoroughly by shaking, and 
stored at 1° C. for approximately 22 hr. It was then centrifuged at about 
2500 r.p.m. for 20 min., the liquid decanted and the precipitate washed once 
with a sodium salicylate-magnesium sulphate solution of the same concen- 
tration as that in which precipitation occurred.” 


Total Nitrogen Determination 


This determination was made by the Kjeldahl method, using mercuric 
oxide as the catalyst. 
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Nitrogen Distribution 

This determination was restricted to amide and arginine analyses as these 
are the main groups that are known to differ quantitatively in so-called 
“glutenin’’ and “gliadin”. The wet gluten or freshly precipitated protein 
was hydrolyzed for 24 hr. with 20% hydrochloric acid, the hydrolysate filtered 
and the residue washed. After making the filtrate up to a volume of 100 ml., 
aliquots were taken for total nitrogen and amide nitrogen determinations. 
The amide nitrogen was determined by the aeration method of Plimmer and 
Rosedale (11) except that magnesium oxide was used instead of calcium oxide. 
The arginine nitrogen was determined on the residue from the amide nitrogen 
determination, again using the method recommended by Plimmer and Rose- 
dale (12). 

Experiments 

Extraction of Gluten 

This experiment was conducted to determine the relation between the per- 
centage of nitrogen extracted from washed gluten and the concentration of 
the sodium salicylate solution employed. The results are shown in Fig. 1. 
Since with the lower concen- 
trations of sodium salicylate 
the gluten ball could not be 
broken up thoroughly, it is 
possible that the percentage 
extracted by these concentra- 
tions is too low. There is a 
difference between the results 
obtained by filtering the ex- 
tract and by centrifuging it, ® 
and the percentage of nitrogen Ca 
extracted by a given concen- = : 
tration of sodium salicylate 2 
obviously depends, to some 2 
extent, on the flour from $ 
which the gluten was washed. 
In all cases, however, the 
relation is regular, there be- 
ing no deviations from the 
smooth curves that cannot 
be attributed to experimental 


error. ~-0--FILTERED 
. —O— CENTRIFUGED 


GLUTEN NIT 


Precipitation of Dispersed 
Protein 
The object of this experi- 


tion between the percentage Fic. 1. Extraction of gluten with sodium sdiedia solutions. 
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of the total nitrogen precipitated and the concentration of the precipitat- 
ing reagent. Only Flours 1 and 3 were used, as sufficient amounts of 
the other two were not available. Dispersions were prepared as already 
described and suitable aliquots pipetted into centrifuge tubes. The protein 
was then precipitated, washed, transferred to Kjeldahl flasks and the total 
nitrogen determined directly. 

The results are shown in Fig. 2. Since the magnesium sulphate solution 
was made up in water, its addition would decrease the concentration of sodium 
salicylate. This secondary factor probably affected the shape of the curves 


5 ——e—- Fiour! 
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CONC. OF MAGNESIUM SULPHATE IN TERMS OF SATURATION 
Fic. 2. Precipitation of dispersed gluten with magnesium sulphate solution. 
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in Fig. 2, but it could in no way affect the continuity. As in the preceding 
experiments, there are minor differences between the results obtained with 
the glutens from the two flours, but the relation between the percentage of 
total nitrogen precipitated and the magnesium sulphate concentration is 
regular, there being no indication of discontinuity. 


Fractionation and Analysis of the Fractions 

Having shown by the preceding experiments that there is apparently no 
discontinuity in the solubility behavior of gluten protein, an attempt was 
made to fractionate the dispersed protein and to determine what chemical 
differences existed between successive fractions. Dispersions were prepared 
and made up to 200 ml. Two 5.0 ml. aliquots were analyzed for total nit- 
rogen, and the remainder transferred quantitatively to centrifuge bottles. 
The dispersions were then adjusted to 0.02 of saturation with magnesium 
sulphate solution and the precipitated protein separated and washed as 
described earlier. The supernatant liquid was decanted into a second cen- 
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trifuge bottle and the washings added. The entire process of adding mag- 
nesium sulphate solution, separating and washing the protein precipitate 
was repeated, and fractions obtained at salt concentrations of approximately 
0.025, 0.035, 0.05, 0.07, 0.1, 0.2 and 0.5 of saturation. The nitrogen 
distribution of these fractions was determined, and the results are plotted 
in Fig. 3. The value represented by the point on the extreme right of each 
of the curves was obtained 
by difference, as explained 
later. The amide and argi- 
nine nitrogen content are 
expressed as percentages of 
the total nitrogen in the 
respective fractions. Each 
of these values is plotted at 
the cumulative mid-point 
of the fraction for which 
it was obtained, the mid- 
points being expressed as 
percentages of the total 
nitrogen in the original 
dispersion. 

In this experiment each 
fraction was exposed to 
sodium salicylate solution 
for one day more than the 
preceding fraction. Conse- 
quently, before drawing any 
conclusions from the results, 
it was necessary to deter- 
mine whether exposure to 
sodium salicylate affected 
the chemical composition of 
the protein. Dispersions were MIDPOINT OF FRACTION AS ‘% OF DISPERSED N 
prepared using gluten from _ Fic. 3. Amide and arginine nitrogen content of gluten 
the four floursand allowed to /79#0"5- 
stand overnight in the refrigerator. After 24 hr. exposure to sodium salicylate 
two aliquots of each dispersion were pipetted into centrifuge tubes. One of 
each pair was immediately brought to 0.2 of saturation with respect to 
magnesium sulphate, and the precipitated protein separated and analyzed 
as already described. The other aliquot of each pair was stored for a further 
nine days, a period corresponding to the longest time any of the fractions 
were exposed, and the protein then precipitated and analyzed as before. 
The results are given in Table II, each value being an average of duplicates, 
and show that neither the percentage precipitated nor the chemical com- 
position of the protein was changed significantly by the prolonged exposure 


Amipe ano arcinive N as %oF-ToTAL N OF FRACTION 
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to sodium salicylate. In the ‘F’’ columns of Table II are given the corres- 
ponding values obtained from the results of the preceding fractionation study, 
by summing the total, amide and arginine nitrogen of the successive fractions 
obtained up to 0.2 of saturation with magnesium sulphate. The amide 
nitrogen results show excellent agreement with those obtained by direct 
precipitation, as do the arginine nitrogen results for Flours 3 and 4. The 
arginine nitrogen values for the other flours appear to be low, a discrepancy 
that can probably be attributed to error. 


TABLE II 


EFFECT ON GLUTEN OF EXPOSURE TO, AND FRACTIONATION FROM, 
10% SODIUM SALICYLATE SOLUTION 


Total nitrogen Amide nitrogen, Arginine nitrogen, 
precipitated, as as percentage of as percentage of 
percentage of total nitrogen total nitrogen 
total nitrogen precipitated precipitated 
in dispersion 


10 


1 Precipitated after exposure for one day to sodium salicylate solution. 
10 Precipitated after exposure for 10 days to sodium salzcylate solution. 
F Fractionally precipitated. 


It is impossible to precipitate all the protein in dispersions in sodium sali- 
cylate by the addition of magnesium sulphate and unfortunately the usual 
protein precipitants, which are acidic, cannot be used on account of the in- 
solubility of salicylic acid. The nitrogen distribution of the protein remaining 
in solution was calculated, however, from the data given in Tables I and II 
and the average of the corresponding values after 1 and 10 days’ exposure of 
the protein to sodium salicylate was plotted, as previously stated, in Fig. 3. 
Except for the percentage of arginine nitrogen calculated for the non-preci- 
pitated fraction of the gluten washed from Flour 4, these results appear to 
fit on the respective curves very well. It is probable that the value 7.6% 
reported in Table I for the arginine nitrogen content of the washed gluten 
is too low, and a slight error in this value would greatly affect the results 
calculated for the non-precipitated residue. 


In order to check these calculated values by direct analysis the soluble 
nitrogenous substances removed in washing the gluten from Flours 2 and 4 
were hydrolyzed and analyzed as before. It is obvious from the results in 
Table I that washing would remove about 10% of the total nitrogen in these 
flours. The amide and arginine nitrogen content, expressed as percentages 
of the total nitrogen washed from the gluten, were as follows: Flour 2, amide 
nitrogen, 7.7%, arginine nitrogen, 16.7%; Flour 4, amide nitrogen, 8.1%, and 
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arginine nitrogen, 17.4%. In each case the amide nitrogen is lower and the 
arginine nitrogen higher than the value represented by the extreme right- 
hand point in the respective curves in Fig. 3. This is in direct agreement 
with the trend of the curves in Fig. 3, and the value 17.4% for the arginine 
nitrogen content o{ the nitrogenous material washed from Flour 4 lends 
support to the contention that the value 5.8% calculated for the non- 
precipitated fraction is in error. 

It is obvious from Fig. 3 that, except for the most soluble 10-15% of the 
protein in washed gluten, the amide nitrogen content increases and the 
arginine nitrogen content decreases with increasing solubility of the protein. 
There are minor differences among the results obtained with the different 
flours but, apart from small deviations which can be attributed to error, the 
curves are smooth and regular with no evidence of discontinuity. The most 
soluble portion of the protein is discussed in a later section. 


Examination of Gluten Fractions and Reversibility of Fractionation 


Having demonstrated in the preceding experiment that gluten protein 
could be fractionated into a number of progressively different fractions, an 
additional experiment was performed in which the fractions were examined 
and then redispersed, the resulting dispersions mixed, and the protein again 
salted out and examined. The object was to determine whether the fractions 
resembled gluten and if not, whether they could be recombined to give a gluten. 

In this experiment a dispersion of gluten from Flour 3 was used. Fractions 
were obtained at 0.025, 0.04, 0.07 and 0.2 of saturation with magnesium 
sulphate, as previously described, and each fraction was washed with dis- 
tilled water, those that were not tenacious being repeatedly rinsed and cen- 
trifuged, until the wash liquid gave a negative test for sulphate ion. 

Fraction 1 contained about 30% of the total nitrogen originally in the 
dispersion and, on washing, it matted into a tough, tenacious, smooth mass. 
Fraction 2 came down as a soft gel which contained about 25% of the total 
nitrogen originally in the dispersion. On washing, the protein exhibited a 
slight tendency to hold together but could not be matted into a tenacious 
ball. Fraction 3 contained about 20% of the total nitrogen originally in the 
dispersion, and precipitated as small particles which showed little tendency 
to cohere. Fraction 4 contained approximately 15% of the total nitrogen 
originally in the dispersion and precipitated as a finely divided powder that 
was not changed during washing. 

After examination, the fractions were redispersed in 10% sodium salicylate. 
Some difficulty was experienced in dispersing Fraction 1, as the tenacious mass 
would not break up and the dispersing agent appeared to attack the surface 
only. On transferring the protein mass to a mortar and grinding with sodium 
salicylate solution, complete dispersion was finally effected. Fraction 2 dis- 
persed with little difficulty and Fractions 3 and 4 seemed quite soluble. The 
resulting dispersions differed considerably in turbidity which varied progres- 
sively from that of Fraction 1, which was quite opaque, to that of Fraction 4, 
which was almost clear. 
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The dispersions of the four fractions were combined, mixed thoroughly 
and the protein then precipitated by the addition of magnesium sulphate 
solution to 0.2 of saturation. The protein came down as a tenacious mass 
which after washing was indistinguishable from one which had not been 
previously fractionated, and distinguishable from freshly washed gluten 
only by its freedom from starch and a slightly increased toughness and 
decreased elasticity. The precipitated gluten closely resembled the residue 
of a washed gluten after extraction with water and, as extracting with water 
dissolves about 10% (Fig. 1) and precipitation at 0.2 of saturation with 
magnesium sulphate leaves about 10% (Table II) of the total nitrogen of 
washed gluten, this similarity is readily explained. 


Discussion 


Before discussing the results in general, the break in the amide and arginine 
curves (Fig. 3) requires explanation. This discontinuity, and the results 
obtained by analyzing the nitrogenous substances removed in the preparation 
of washed gluten, indicate that the most soluble protein differs from the 
remainder. This fraction represents 20-30% of the protein in flour, only 
about half of which is included in washed gluten as prepared in these ex- 
periments. The solubility curves (Figs. 1 and 2) do not exhibit any breaks, 
which may suggest that this most soluble fraction is not distinct from the 
rest of the protein on the basis of solubility. The curves in Fig. 2, however, 
are flattening out in the region where 85-90% of the protein is precipitated. 
In extracting washed gluten it is impossible, as previously stated, to break 
up the gluten ball with the lower concentrations of sodium salicylate, so the 
lower points on the curves in Fig. 1 may not be sufficiently accurate to exhibit 
a slight break if one did exist. 


It cannot be stated from the results of the present study whether this most 
soluble 20-30% of the protein in flour is one distinct protein, or several, or a 
protein complex, but this portion appears to be distinctly different from the 
rest of the protein. The terms “‘non-gluten” and “‘salt soluble” are applied 
to the more soluble protein in wheat flour, and this fraction is commonly, 
though arbitrarily, as Blish (1) points out, estimated by extraction with 5% 
potassium sulphate. This salt solution extracted 13.4 and 11.4% of the 
total nitrogen of Flours 2 and 4, respectively, which is much less than the 
amount which appears, from the results of the present study, to be distinctly 
different from the rest. 


None of the results of this investigation lend any support to the classical 
view that gluten is composed of glutenin and gliadin. The results in Fig. 1 
show that, by choosing a solution of sodium salicyJate of suitable concentration, 
practically any desired percentage of the gluten can be extracted and that 
this percentage is approximately constant with the different flours used. 
The curves in Fig. 1 are regular and if glutenin and gliadin are distinct proteins 
the curves might reasonably be expected to exhibit a break in the region of 
50% extraction. The curves in Fig. 2 are also regular, so on the basis of 
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solubility in sodium salicylate there is no evidence of the existence of two 
distinct gluten proteins. This is in agreement with Cook’s conclusion (2) 
based on the precipitation of gluten protein from a dispersion in urea solution. 

It is also well known that so-called ‘‘glutenin’’ and ‘“‘gliadin”’ differ in 
chemical composition. The results in Fig. 3 show, however, that by suitably 
adjusting the magnesium sulphate concentration a number of fractions can be 
obtained which vary progressively and uniformly in chemical composition. 
The most soluble 10-15% of the dispersed protein appears to be different and 
has already been discussed, but there is no evidence of any break in the 
curves that could be taken to indicate that gluten is composed of glutenin and 
gliadin or of glutenin, ‘‘mesonin’’ and gliadin as Sandstedt and Blish (14) 
suggested. 


The fact that no break which could be taken to indicate that gluten is 
composed of so-called “‘glutenin” and ‘‘gliadin”’ was found in any of the curves 
in Figs. 1, 2 and 3, might be explained on the basis that one protein is adsorbed 
on the other to such an extent that extraction with, or precipitation from, 
sodium salicylate solutions, merely removes fractions containing progressively 
less of one and more of the other protein. Were such an explanation valid, 
however, then that fraction which contains the two components in the same 
ratio as they are present in gluten should, of necessity, have the same physical 
and chemical properties as gluten. None of the fractions obtained in the 
preceding experiment exhibited a close similarity to washed gluten. Bearing 
in mind the fact that Fraction 1 contained the most insoluble 30% and Frac- 
tion 2 the next 25% of the total nitrogen in the original dispersion, reference 
to Fig. 3 and Table I will show that Fraction 2 had approximately the same 
chemical composition, as determined by amide and arginine analyses, as the 
original washed gluten. Its physical properties, however, were not at all 
similar to those of gluten since, as previously stated, it could not be matted 
into a tenacious mass. The difference in physical properties between Frac- 
tion 2 and gluten cannot be attributed to the action of sodium salicylate 
solution in the former case since, when the fractions were redispersed, com- 
bined and the protein precipitated, it came down as a gluten. 


A view which will explain the results obtained in this investigation is that 
gluten, or at least the major portion of it, is a single protein complex which 
may be separated into a great many fractions which differ progressively and 
systematically in both physical and chemical properties. The solubility of 
the fractions and the turbidity of the resulting dispersions indicate that the 
protein particles decrease in size with increasing solubility. Sorensen (15) 
regards soluble proteins such as egg albumin and so-called ‘‘gliadin’’ as 
“reversibly dissociable component systems’. The results of the present 
investigation show that gluten is reversibly fractionable, and in terms of 
Sorensen’s nomenclature, gluten may be regarded as a “highly dissociable 
component system.” 

This view is in accord with the results obtained by Gortner, Hoffman and 
Sinclair (8) in their extraction studies, and will account for some of the wide 
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deviations in amide and arginine content of so-called ‘‘glutenin’’ preparations, 
as different methods doubtless yield samples containing different fractions 
of the gluten complex. For example, the average amide and arginine content 
of glutenin prepared by Cook and Alsberg (3) from dispersions in urea were 
respectively 19.5% and 9.8%, while the corresponding values obtained by 
Larmour and Sallans (10) for glutenin prepared from a dispersion in 0.007 N 
acetic acid were 21.0% and 8.5%. The average amide and arginine content 
of the most insoluble 40% of the gluten in the present study (Fig. 3) were 
19.6% and 9.6%, which values are in good agreement with those obtained 
by Cook and Alsberg. Larmour and Sallans, however, failed to secure 
complete dispersion of the gluten in the reagent used, and the undispersed 
portion was undoubtedly that portion which was lowest in amide and highest 
in arginine nitrogen. The failure to include this in the prepared glutenin 
would result in a higher amide and a lower arginine content than would 
otherwise have been obtained. The differences in the results’ presented by 
Cook and Alsberg and Larmour and Sallans are, therefore, explained. Doubt- 
less this view will also explain the variability in the chemical composition of 
glutenin prepared by other workers. 


The results of the extraction studies in this and previous investigations (8) 
show that definitions of flour proteins based on solubility in a solution such 
as 70% alcohol are purely arbitrary, as the changes in physical and chemical 
properties are not sudden but gradual and regular. The precipitation ex- 


periment in the present study appears to give a much better picture of the 
nature and composition of gluten, and it therefore seems preferable that the 
terms ‘“‘glutenin’’ and ‘‘gliadin’’ should be discarded. These suggest that 
gluten is composed of two distinct proteins, whereas the results of the present 
study, and those of earlier investigations (2, 8), substantiate the view that 
it is a single protein complex which can be divided into a great many progres- 
sively different fractions. 
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DENSITY AND ADSORPTION STUDIES 
IN THE REGION OF THE CRITICAL TEMPERATURE: 
SYSTEM DIMETHYL-ETHER-ALUMINA! 


By J. Epwarps? AND O. Maass? 


Abstract 


An apparatus and an experimental procedure for simultaneously determining 
the density and adsorption of liquids, vapors and gases in the region of the 
critical temperature are described. Measurements of adsorption and density 
of dimethyl ether through the temperature range 124-135° C. are given for the 
liquid, vapor and gas. Adsorption of the liquid does not occur below 124° C. but 
beginning at this temperature, it rapidly increases to a maximum, beyond which 
the adsorption curve resembles a normal isobar. It is shown that the apparent 
adsorption cannot be ascribed to a lag in the critical temperature of liquid held 
in the pores, and further that an absurd value for the density of such liquid must 
exist to account for the total adsorption observed. Hysteresis effects of a complex 
nature are noticed, and these are discussed from the point of view of the hypo- 
thesis that: ‘there is an essential difference, other than concentration, between 
a liquid and a highly compressed gas.” 


Introduction 


The series of researches that has been made in this laboratory on the 
behavior of liquids and gases in their critical regions was recently extended 
by Morris and Maass (8) to the study of adsorption for the system propylene— 
alumina. Inasmuch as the results of this investigation were of sufficient 
interest to warrant a more exhaustive survey of the phenomenon of adsorption 
at high pressures, and particularly in the region of the critical temperature 
and pressure, adsorption isotherms from —78 to 130° C. through a pressure 
range 0.1-52.6 atm. were determined for the system methyl-ether—alumina. 
The results of this investigation will be published by the authors in a later 
paper. 

On account of the large error involved in the estimation of the density of 
the methyl ether surrounding the gel above the critical temperature, the 
buoyancy corrections applied to the apparent weight of the gel yielded values 
for the adsorption which were obviously greatly in error. An apparatus 
was therefore constructed which would allow simultaneous determinations of 
density and adsorption to be made under all the experimental conditions 
imposed. Quite apart from the significance which this method gave to the 
adsorption data, the variation of density with temperature under a variety 
of conditions proved to be of considerable interest in the further elucidation 
of the nature of the liquid-gas equilibrium at the critical temperature. In 
this paper it will be shown that results have been obtained which cannot be 
satisfactorily explained on the basis of continuity of state. 


1 Manuscript received January 4, 1935. 
e Contribution from the Physical Chemistry Laboratory, McGill University, Montreal, 
‘anada. 
2 Demonstrator in Physical Chemistry, McGill University. 
8% Macdonald Professor of Physical Chemistry, McGill University. 
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Experimental 
Description of the Apparatus 

The bomb (Fig. 1) consisted of specially selected heavy-walled Pyrex 
tubing, with the arm D of 14 mm. inside diameter and wall thickness 4 mm., 
while A and B were of 11.5 mm. inside diameter and 3.5 mm. wall thickness. 
Connection among the three arms was 
made with capillary tubing of 1.2 mm. 
bore. 

The bomb A contained a quartz spiral 
suspended by a hook from a glass collar 
which was held snugly in place by a 
small indentation in the bomb. The 
alumina gel, contained in an aluminium 
basket perforated with small holes, was 
hung from the bottom ring of the quartz 
spiral by means of a fine aluminium 
wire fastened to the basket. The gel 
was prepared according to the method 
of Munro and Johnson (9) and subse- 
quently heated in a muffle furnace to 
600° C. The quartz spiral had a sensi- 
tivity of 0.00366 gm. per mm. extension, 
and the maximum load was 0.50 gm. 

The bomb B was of the same dimen- 
sions as A, and differed from A only by 
the substitution of a float for the sor- 
bent. The sensitivity of the spiral was 
0.00517 gm. per mm. extension. The 
float occupied a volume of 0.8971 cc. 
and weighed 0.3769 gm.: the density of the float was thus 0.4201. If a 
liquid of density greater than this surrounded the float, it would rise through 
and break the spiral; measurements of the density and adsorption of liquid 
methyl ether were therefore confined to the region above 121° C. Since the 
adsorption values up to the critical temperature had previously been obtained 
very accurately by applying the density values of Cardoso and Coppola (2), 
Cardoso and Bruno (1), and Tapp, Steacie, and Maass (13), this lack of flexi- 
bility in the apparatus was of no real concern in the present work. 

The connecting capillary tubing was covered with a thin layer of asbestos 
paper, about which was wound a length of resistance wire. The temperature 
of the capillary tubing could thus be maintained from 10 to 20° C. above 
that of either arm. 

The total volume of the apparatus was approximately 96 cc., distributed 
so that the volume of D was nearly equal to the combined volumes of A and B. 

Connection was made through C to a high-vacuum line and to a fractional 
distillation unit. The complete lengths of D, and A and B up to the top of 








Fic. 1. The bomb. 
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the gel and float were maintained at a temperature of 360° C. while the system 
was alternately evacuated and flushed with air dried over phosphorus pent- 
oxide. The dry, evacuated system was then flushed a number of times with 
methyl ether, and an amount sufficient to give a meniscus disappearing near 
the top of D, or, alternatively, of A and B, was condensed into D. The 
methyl ether was a thrice fractionated portion of 99.95% pure stock supplied 
in a steel cylinder by Ohio Chemicals Ltd. The bomb was sealed off at C 
and the arm D immersed in liquid air while A and B were held at a temperature 
of 140° C. Readings of the spiral extensions in both of these bombs then 
gave a check on the weights of the float and gel determined previous to assem- 
bly of the apparatus. Though the weights of the float as determined by both 
methods agreed within less than 0.05%, the weight of the gel as determined 
by spiral extension was appreciably lower than the original weight. This 
loss in weight was attributed to desorption of water vapor during the 
process of alternate evacuation and flushing with dry air described above. 
The actual weight of sorbent as determined after assembly of the apparatus 
was 0.2829 gm., which included 0.0473 gm. of metallic aluminium comprising 
the basket and wire. The volume of the complete unit was 0.0842 cc. This 
value was obtained from the known density of aluminium and the density of the 
gel measured by water displacement. The latter value was 3.44 gm. per cc. 


The apparatus was immersed in two thermostats so that the temperatures 
of D and of AB could be separately varied and controlled. (In future the side 
of the apparatus containing the gel and float will be designated ‘‘AB”’, while 
the other side will be referred to as ‘‘D’’). The thermostats consisted of 
Pyrex jars 10 in. in diameter and 18 in. in height, free from striations in 
the walls, and filled with dibutyl phthalate to a level about one inch below 
the horizontal capillary tube connecting the arms. This liquid possessed 
the advantages of transparency and low volatility at the temperatures 
required. Suitable heaters, copper cooling coils and stirrers permitted tem- 
perature control to within 0.1° C., the limit of accuracy of the calibrated 
thermometers which were read through a telescope situated alongside the 
cathetometer. Properly adjusted lighting permitted a clear vision of the 
ends of the spiral. 


The complete apparatus was encased in a wooden box reinforced with strap 
iron and containing three built-in shatter-proof glass windows. The possi- 
bility of damage to the cathetometer and telescope was minimized by placing 
a plate glass window directly in front of these instruments. 


Experimental Procedure 


In general, two procedures were adopted. The arm D was held at some 
arbitrarily fixed temperature above the critical temperature (126.9° C.), 
e.g., 130.0° C., and AB heated slowly from 122 to 135° C. and cooled to 
122° C. Readings were taken at approximately 1° intervals. That tem- 
perature equilibrium between the thermostat and the contents of the bomb 
had been attained was assured by following the change in the position of the 
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float after the thermostat had reached the required temperature. When, 
after a 10 min. interval, no change in the position of the float was detectable, 
the spiral extensions in both bombs were determined. In so far as the attain- 
ment of true equilibrium in adsorption was concerned, test experiments showed 
that no change in the position of the basket was detectable after 30 min. 
at true temperature equilibrium. 

The second procedure adopted was the reverse of the one described above. 
The arm AB was held at a fixed temperature, e.g., 130° C., and the temper- 
ature of arm D varied through the range 122-135° C. 

Values for the apparent adsorption were obtained by means of the equation 

W, = (A—N)S+dmnv,—w 
where w = the true weight of the gel, d,, = density of the medium surrounding 
the gel, v, = volume of the gel, A = length of the spiral under experimental 
conditions, NV = normal length of the spiral under no load, and S = sensitivity 
of the spiral. It is seen that the term d,,v, is the buoyancy correction for 
the medium displaced by the gel, and that (A—N)S is the apparent weight 
of the gel. 

The expression for W,, the apparent adsorption, does not take into account 
the correction due to the buoyancy effect of the adsorbed material itself, 
so that the true adsorption, W = W.+d,c, where v, is the volume of the 
adsorbed phase. Since v = W/d,, where d, is the density of the adsorbed 


TABLE | phase, 
Wode . 
RESULTS OF DENSITY DETERMINATIONS We= see 


I sree. A full discussion of the 


Temp. Temp. | 
of AB, Density || of AB, | Density | of AB, Density _use of this equation in esti- 
" . it * . * . 


| 
; : 
mating the true adsorption 
Arm D at 128° C. will be presented later. 
124.0 | 0.3740 II ’ ' The total error in the de- 
125.4 | 0.3599 : : ; termination of W, did not 


27:2 | 0.3070 tis 21> exceed 0.0008 gm. 


128.3 | 0.2280 || 130.6 | 0. 410. The density of the 
A 5 D at 130° C medium surrounding the 
sans : float was estimated to with- 


123.7 | 0.3805 in +0.0005 gm. per cc. 
123.9 | 0.3762 
125.2 | 0.3559 
125.4 | 0.3573 
125.8 | 0.3499 


126.8 | 0.3345 The results obtained for 
127.4 | 0.3274 


128.2 | 0.3074 ‘1973 ~ the densities are given in 
129.5 | 0.2744 2194 Tables I and II, and in 
ArmDati27’C. Figs. 2 and 3; and those for 


the adsorption, in Tables 


125.0 | 0.3560 0.2465 || 129.0 | 0.1842 
126.4 | 0.3300 128.0 | 0: 1960 | III and IV and Figs. 4 
and 5. 
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TABLE II 


Adsorption is expressed 
RESULTS OF DENSITY DETERMINATIONS 


in the usual manner as 























x/m, where x is the weight . . 

: ‘emp. emp. emp. 
adsorbed, and m is the of D, | Density|} of D, | Density || of D, | Density 
weight of adsorbent. <. % <. 

In all the curves, the Gel and float (Arm AB) at 128° C. 
solid lines represent valaes © @_———___£__ A __—__ 
123.9 | 0.1493 129.7 | 0.3150 128.8 | 0.2933 


obtained in the first half = 495"5 | 0'1646 | 131.2 | 0.3277 || 127.7 | 0.2239 


of the experiment, i.e., 126.7 | 0.1870 134.8 | 0.3460 125.6 | 0.1665 
h 127.6 | 0.2172 || 132.7 | 0.3378 
when temperature was ap- 128.7 | 0.2875 || 130.8 | 0.3249 | 


proached from below, while . 
the dotted lines indicate Ge! and float (Arm 4B) at 130° C. 


the values obtained dur- m2 o-1080 130.4 9.2008 429.5 0.2399 
‘ ‘ , 1 131. . 127.6 | 0.1788 
ing the cooling, or second = 497'3 | 9/1762 || 134.7 | 0.3272 || 125.6 | 0.1558 


part of the experiment. 128.2 | 0.1945 || 132.0 | 0.2966 
The hysteresis obtained 129.5 | 0.2288 130.5 | 0.2756 

in the density curves 

cannot be explained on the basis of a temperature lag, although such a lag 
would tend to bring the curves together. In order to demonstrate this, 
a procedure such as the following was adopted: If in the early part of the 
experiment, the arm AB was at, say, 127.0° C., and a measurement was to 
be taken at 128.0°, the arm would be heated to 128.2° and held at this tem- 
perature until the float was stationary for a 10 min. interval, and then 
allowed to cool to 128.0° C. When the float was again stationary, temperature, 
density and adsorption measurements would be made. In the second half of 
the experiment the arm would be cooled from a higher temperature, 
e.g., 129.0° C., to 127.8°, and subsequently reheated to 128.0°. In this 
manner any temperature lag was in a direction tending to bring the two curves 


128 
TEMPERATURE, “Cc. 


Variation with temperature of the 








Fic. 2. 


density of material in AB. Solid curves: tem- Fic. 3. Variation, with temperature of 
perature approached from 124° C. Broken arm D, of the density of the material in arm 
curves: temperature approached from 134.7°C. AB. Solid curve: temperature approached 


I. Arm D held at 127.0° C. II. Arm D from below. Broken curve: temperature 
held at 128.0° C. III. Arm D held at approached from above. I. Arm AB held at 
130.0° C. 128.0° C. II. Arm AB held at 130.0° C. 
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together, so that the recorded values indicate the minimum possible hysteresis. 
Since the error in temperature measurement was 0.1° C., the heating and 
cooling curves could be made to coincide in the regions of small hysteresis 
by choosing values within the limits of error. However, the tabulated 
results were chosen from a large number of observations for which in no case 
did the hysteresis disappear. The experiments were all repeated from three 
to five times. 


In Fig. 2 the density-temperature relations of arm AB are given for three 
fixed temperatures of D. The incomplete and unchecked Curve I has been 
included to indicate the approximate density behavior in AB when D is 
held at the critical temperature. An important consideration arising out of an 
inspection of Fig. 2 concerns the hysteresis effects in Curves II and III. Tapp, 
Steacie, and Maass (13), Winkler and Maass (16), and Morris (7) found 
that the densities existing above and below the level of disappearance of the 
meniscus differed by a considerable amount at the critical temperature; the 
density in the region previously occupied by liquid being greater than that in 
the region previously occupied by vapor. It was further found that on 
increasing the temperature of the bomb this difference in density rapidly 
decreased, finally disappearing about 14° C. above the temperature at which 
the meniscus disappeared, and that on re-cooling, even to the critical tem- 
perature, the density difference did not reappear, that is, on re-cooling, the 
density of material below the level of meniscus disappearance remained 
steady down to the critical temperature. In the apparatus used in the 
present investigation the process is complicated by the temperature variations 
between the two arms, so that a large variation in the distribution of material 
results. Nevertheless, the hysteresis shown in Fig. 2, Curves II and III, 
represents an exact parallel to the results of other investigators, and serves to 
corroborate their data. 


It will be seen that the procedure adopted to yield the data in Table II and 
illustrated in Fig. 3 gives a picture of the behavior of the density of the material 
in D corresponding to Fig. 1 for AB. Thus the hysteresis in density in Fig. 2, 
showing a lower density upon re-cooling, is reflected by a hysteresis in D 
such that a higher density is obtained on re-cooling. This hysteresis is not 
considered important except in that it corresponds to the different distribution 
of material between the two arms, and that it gives further evidence for the 
existence of the hysteresis in Fig. 2. 


In Curve II, Fig. 2, it will be seen that the density (solid line) at 128 .0° C. is 
0.2505, while the corresponding density in D, also at 128.0° C., is 0.2420. 
Under these conditions of equal temperature, the apparatus may be con- 
sidered as a straight bomb, AB representing the region below the point of 
disappearance of the meniscus. Thus the density difference obtained by 
other investigators is again confirmed. That this density difference decreases 
rapidly with increasing temperature may be seen by comparing Curve III, 
Fig. 2 and Curve II, Fig. 3. For instance, the density of material in AB at 
130° is 0.251, while the density of material in D, also at 130° C., is 0.2505. 
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Although the decrease in density difference is somewhat greater than would 
be expected by comparison with the results obtained by Tapp, the parallelism 
between the two sets of results is a convincing substantiation of their 


reliability. 


The change in density 
of liquid below the criti- 
cal temperature appears, 
from Fig. 2, to be depend- 
ent on the temperature 
of bomb D. The possible 
significance of this fact, 
together with the nature of 
the curves above the critical 
temperature will be fully 
discussed later in the paper. 


The results of the ad- 
sorption measurements are 
given in Tables III and IV 
and in Figs. 4 and 5. Ad- 
sorption is W, per gram 
of adsorbent. In Fig. 4, 
the apparent adsorption, 
plotted against tempera- 
ture of the arm AB, corre- 
sponds to the density curves 
of Fig. 2, while Fig. 5 cor- 
responds to Fig. 3. It 
will be noticed (Fig. 4) that 
the adsorption of liquid 
methyl ether up to 124° C. 
is zero, and that adsorp- 
tion, which begins at about 
this temperature, increases 
with increasing tempera- 
ture, reaching a maximum 
at a temperature approxi- 
mately 0.3° C. higher than 
the temperature of arm 
D. From this point the 
adsorption-temperature 
curve resembles the normal 
gas or vapor adsorption 
isobar. The second, or 


cooling part of the experi- 
ment, yields a complex 












































TABLE III 
RESULTS OF ADSORPTION DETERMINATIONS 
Temp. Temp. Temp. 
of AB, | Adsorp- ||} of AB, | Adsorp- || of AB, | Adsorp- 
"a tion "Se tion ae tion 
Arm D at 128° C. 
124.0 {| 0.0006 130.4 | 0.0251 128.3 | 0.0267 
125.4 | 0.0046 132.9 | 0.0201 127.7 | 0.0187 
126.8 | 0.0095 134.7 | 0.0173 126.3 | 0.0053 
127.6 | 0.0134 133.5 | 0.0210 125.1 | 0.0028 
128.3 | 0.0283 130.6 | 0.0255 
Arm D at 130° C. 
124.2 | 0.0018 130.9 | 0.0268 130.2 | 0.0226 
125.6 | 0.0057 132.0 | 0.0222 129.1 | 0.0191 
128.3 | 0.0125 132.7 | 0.0216 128.4 | 0.0130 
128.8 | 0.0138 134.7 | 0.0187 126.4 | 0.0065 
129.5 | 0.0188 132.5 | 0.0237 124.3 | 0.0007 
Arm D at 127° C. 
Temp. of Adsorption Temp. of | Adsorption 
AB, °C. 4B, CS. 
124.8 0.0000 128.0 0.0283 
125.9 0.0104 129.0 0.0255 
127.0 0.0180 
TABLE IV 
RESULTS OF ADSORPTION DETERMINATIONS 
Temp. Temp. Temp. 
of D, | Adsorp-|} of D, | Adsorp- |} of D, | Adsorp- 
se tion bo tion Nase tion 
Gel and float (Arm AB) at 128° C. 
123.6 | 0.0241 132.0 | 0.0108 127.2 | 0.0295 
125.7 | 0.0279 134.7 | 0.0058 126.4 | 0.0317 
127.0 | 0.0288 131.3 | 0.0127 125.6 | 0.0270 
128.3 | 0.0198 129.5 | 0.0167 
129.6 | 0.0149 128.3 | 0.0242 
Gel and float (Arm AB) at 130° C, 
123.5 | 0.0162 129.5 | 0.0215 129.5 | 0.0226 
125.7 | 0.0208 130.4 | 0.0180 127.7 | 0.0230 
127.3 | 0.0224 131.3 | 0.0166 126.4 | 0.0223 
128.2 | 0.0225 134.7 | 0.0085 125.1 | 0.0197 
128.3 | 0.0227 130.5 | 0.0194 
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type of hysteresis which cannot be satisfactorily explained on the basis of a 
lag in equilibrium. Thus from 134.7° C. to 131 or 132° the adsorption is 
greater than that obtained earlier, 
while the curve has a negative slope. 
Provided equilibrium had not been 
attained, the adsorption values at 
true equilibrium would result in a 
greater hysteresis than has been indi- 
cated. The second portion of the 
cooling curve, on the other hand, can 
be explained on such a supposition, 
but it appears unusual that such 
behavior should be apparent in a 
narrow region only: The hysteresis 





124 126 


Fic. 4. Variation with iaaperoters the 
adsorption from material in AB. Solid and 
broken curves have the same «poe as in 
Figs. 1 and 2. I. Arm D held at 127.0° C. 
II, Arm D held at 128.0° C. III. Arm D 
held at 130.0° C. 


at low temperature can also, of 
course, be explained by a lag in 
equilibrium. -An interpretation of 
these curves predicting their form 
and the hysteresis as found will be 


presented later. 

In Fig. 5 the situation is much 
clearer. To begin with, the curves 
are adsorption isotherms, the gel 
being held at 128 and 130° C. The 
low-temperature (of D) part of the 
curve merely represents the normal 
adsorption-pressure curve, the pres- 
sures being the vapor pressures of 
liquid methyl ether from 124° C., 
to the critical temperature. The 
change in slope of the curve above 
the critical temperature will be 
treated later, as also will the hyster- 
esis, which is not nearly as complex 
as in Fig. 4. 





Fic. 5. 
material in AB with temperature of arm D. 
Solid and broken curves have the same pee 


Variation of the adsorption from 


cance as in the previous figures. I. AB held 


128.0° C. II. AB held at 130.0°C. 


Density Discussion 


The primary object of this investigation was a study of adsorption at the 
critical temperature, but inasmuch as the density measurements recorded 
concurrently are of interest in themselves, some consideration may be given 
them at this stage. In the first place they confirm previous work carried 
out in this laboratory on the persistence of density differences above the 
critical temperature. In the second place density differences in the liquid 


below the critical temperature, depending on experimental conditions other 
than temperature, have been recorded for the first time. 
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For clarity of discussion, a hypothesis of the nature of the liquid state to 
explain these and related phenomena, and which will be elaborated more fully 
later on, may be stated here:— “There is an essential difference other than 
concentration between a liquid and a highly compressed gas. This difference 
may be ascribed to a structure in the liquid, whereas in a gas a completely 
chaotic distribution of molecules exists.” 

From the hysteresis obtained in the density values it is seen that the density 
of the material in the region previously occupied by the liquid depends on the 
extent to which the liquid is heated above the critical temperature (Fig. 2). 
Both in nature and magnitude of the effect this observation agrees with the 
results obtained in this laboratory by Tapp, Steacie and Maass (13), Winkler 
and Maass (16), and Morris (7). On the basis of the above hypothesis the 
density phenomena up to 14°C. above the critical temperature may be 
ascribed to the persistence, in part at least, of the structure of the liquid. 
This persistence was shown to be unimpaired by stirring of the material. On 
cooling from a temperature appreciably above the critical the re-appearance 
of the structure in the above sense does not occur until the critical temper- 
ature is passed; hence the hysteresis (Fig. 2, Curves II and III). 

Fig. 2 shows that density differences in the liquid are appreciable for a full 
degree below the critical temperature. The condition bringing about a 
decrease in density is the lower temperature of the arm D, the temperature 
of which governs the amount of material held in it. This immediately 
raises the question whether the additional weight of material above the point 
in the liquid where the density is recorded might bring about such a result. 
It was observed that the maximum difference in height of liquid was only 5%. 
It is true that, using the calculations by Gouy (4), a 10% difference in density 
would follow from an increase of 5% in height of liquid, provided the 
measurements ‘were made within 0.001° C. of the critical temperature, but, 
as Gouy observed, at 0.1° C. below the critical temperature, the gravitational 
effect should produce no appreciable density difference, whereas, as will be 
seen from Fig. 2, a 10% difference in density exists at 0.5° C. below the critical 
temperature. As this has been recorded for the first time and the apparatus 
was not suitably designed for very accurate measurements in this region, further 
discussion of this particular result will await more complete experimental 
evidence. The meniscus level was always above the top of the spiral, so that 
while a second order correction to the density values due to the buoyancy effect 
on the spiral should be made, the relative density values would be unaffected. 

The shift in the curves with increasing temperature of the arm D, and espe- 
cially the manner in which the inflection at the equalization of temperature 
decreases with increasing temperature of D are interesting. Since the tem- 
perature of D governs the distribution of material between the two arms, it 
was to be expected that the curves should be shifted in the manner illustrated, 
although it would be logical to expect no inflection in the curve. The meaning 
of this inflection may be closely linked with the sharp maxima obtained in the 
adsorption measurements from which a clearer picture of the true state of 
affairs can be deduced. These will be discussed later. 





366 CANADIAN JOURNAL OF RESEARCH 


Regarding Fig. 3, the experimental conditions involved nochange from liquid 
to gas, but rather from vapor to gas, so that no particular interest attaches 
itself to the form of these curves. They were necessary, of course, for the 
calculation of the adsorption, and otherwise merely show the distribution of 
material between the two arms. 


Adsorption 


The most interesting observations made on adsorption involved the transi- 
tion of liquid to gas, the results of which have been given in Table III and 
Fig. 4. A number of arguments may be put forward to account for the nature 
of these curves, and these will be the first object of discussion. In the first 
place the curves represent a distinct discontinuity in state in the region of 
the critical temperature. From the point of view of Le Chatelier’s principle, 
adsorption will always decrease with temperature since adsdrption is accom- 
panied by the evolution of heat. Actually a region is seen to exist in which a 
positive adsorption-temperature derivative may be obtained. Provided, 
then, that these curves are a true measure of the adsorption, definite state- 
ments may be made regarding the magnitude of the homopolar attractive 
forces between the molecules of sorbate in the vicinity of the critical tem- 
perature. In the liquid below 124° C. these forces exceed the attractive 
forces between sorbate molecules and surface of sorbent. For higher tem- 
peratures this is no longer true, so that adsorption becomes apparent. The 
obvious explanation lies in a change in the nature of the molecular attractive 
forces at the critical temperature. 


Following the region through which adsorption increases with increasing 
temperature there is seen to be a region in which adsorption falls off with 
temperature identically as does a normal adsorption isobar. The cooling 
part of the experiments yields complex hysteresis effects which, together with 
several other points of importance, will be dealt with later. 


The question may now be raised whether the curves in Fig. 4 represent the 
true adsorption. It must be remembered that the apparent adsorption, W., 
does not take into account the buoyant effect of the adsorbed phase itself, 
and that the true adsorption, W, is given by 

Wede 
"°t- 

The density of the adsorbed phase will depend to some extent on the density 
of the medium and furthermore will decrease slightly with increasing temper- 
ature, but the total change in d, should be very small over a narrow temper- 
ature range. Hence the validity of the curves in Fig. 4 may be tested by 
estimating W and calculating d, from the above equation. 


Let us assume that no discontinuity in true adsorption accompanies the 
transition of liquid to gas. To get an approximate value of W in the region 
of 124° C. the normal isobar part of the adsorption curve may be extrapolated 
to 124° C., and from the resulting values of W and the known values of W, 








DENSITY AND ADSORPTION NEAR CRITICAL TEMPERATURE 367 


and d,, the change in d, with temperature may be shown as in Fig. 6. It 
will be noted that a distinct discontinuity at the critical temperature is 
obtained. Thus, assuming that no discontinuity in adsorption exists, there 


results a discontinuity in the density 
of the adsorbed phase, a condition 
which cannot be logically expected 
under any circumstance. 

A more direct method of attack 
lies in the determination of the true 
volume of the adsorbed phase, and 
although this cannot be estimated 
directly, the pore volume of the gel is 
fairly accurately known. McGavack 
and Patrick (6) hold with the hy- 
pothesis that all of the sorbed 
material in a porous body is liquid 
condensed in the pores. Earlier 
work by the authors on the adsorp- 
tion of methyl ether below the 
critical temperature showed that, in 
accordance with Patrick’s results, 





Fic. 6. The necessary variation of the 
density of the adsorbed phase with temperature 
on the assumption that adsorption from liquid 
and from gas is continuous. Results cal- 
culated from Curve II, Fig. 4. 


the adsorption tended to be independent of temperature at equal relative 
pressures, p/p,, where ~, is the vapor pressure of the liquid. 

While the Kelvin equation may be used to obtain the radius of the pores, 
the pore volume may be estimated only by extrapolation of the adsorption 
to the saturation pressure. Inasmuch as the extrapolation was a small one, 
and the adsorption was fairly independent of temperature at equal relative 
pressures, Wecturation Could thus be estimated to well within 10%*. 

Assuming that the liquid condensed in the pores has a density equal to 
that of bulk liquid at that temperature, the volume of the pores may be 
obtained as v, = W,/d;. The value obtained was 0.09 cc. per gram of gel. 

If now we consider the form of the adsorption curves in Fig. 4 to be due 
to a lag in the critical temperature of material in the pores, the density of 
this material can be calculated from the known values of W,, dm, and 
(v, = volume of the pores). The first calculation is the true adsorption, and 
it is immediately seen that in this case W is obtained by means of Definition A 
of Coolidge (3) in which the volume of the gel is ¢:+¢2, 7.e., the volume of 
“skeleton” plus the pore volume. The resulting curve, corresponding to 
Curve II, Fig. 4, is given in Fig. 7. Assuming the material adsorbed to be 


all confined in the pores as liquid, the density of the material is given by 
W/0.09, so that a discontinuity in density of the liquid in the pores would 
have to be admitted. Considering that W, was consistently zero below 
124° C., the true adsorption according to Definition A of Coolidge must be 
considered as paralleling the liquid density curve in the region below 124° C. 


* The extrapolations were made for isotherms at 110, 98, 82, and 70° C., the final value of 
the pore volume being consistent to less than 5%. 
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Obviously liquid contained in the pores 
when the gel is surrounded by material 
of the same density cannot be considered 
as sorbed—the medium merely wets 
the gel. Even so, the form of the 
curve above the critical temperature 
cannot be explained by persistence of 
liquid in the pores above the critical 
temperature. That this statement is 
true can best be demonstrated by 
observing the necessary discontinuities 





as in the density of the adsorbed phase 
Pervert cn at 124° C. and 128° C. (Fig. 7). 
Fic. 7. Variation of adsorption with tem- , » : 
perature, taking the total volume of the gel as It may be definitely shown in another 


oe “Calculated from Fig. 4. ius tore way that the observed adsorption above 


the critical temperature cannot be due 
to the persistence of liquid in the pores. In the region below 124° C. the 
gel is merely wet by the liquid, so that the density of material in the pores 
is equal to the density of the bulk liquid. Assuming that above a certain 
temperature near the critical one, the change in density of material in the 
pores lags behind the rapidly decreasing liquid density, it must be admitted 
that the density of material in the pores still decreases with increasing 
temperature. By extrapolation of Fig. 4, Curve II, to 124° C. in the same 
manner as for Fig. 6, a value may be obtained at 128 or 130° C. which will 
lie near the density of the material held in the pores. In this way an 
approximate relation for the variation 
of the density of adsorbed material 
with temperature may be obtained 
for the range 124-135° C. Applying 
these values, the true adsorption may 
be calculated, in which case Fig. 8 
is obtained. From the values of W 
and d, it is seen that in the region 
near 124° C. the adsorbed material fem 
must occupy a volume 10 times that 
of the pores, or conversely, occupy- 
ing the pore volume, the material 
would have to have a density 10 times = 
the density of liquid at a lower tem- 
perature. This is an absurd value 
again, and we are led to the conclusion 
that the adsorption curves in Fig. 4 - ES. _ 
cannot be accounted for, except in small Fic. 8. 





Variation of adsorption with tem- 


art, b rsistence of liquid in the perature, assuming a continuous change in 
P y : q : the density of the adsorbed phase. Calculated 
capillaries by virtue of an increased from Curve II, Fig. 4. 
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critical temperature. Further experimental work has yielded interesting 
points concerning this problem, and in a later paper a full description will 
be given in the light of the conclusions of Patrick, Preston and Owens (10) 
that the critical temperature in fine capillary tubes is raised. 

It would appear then that, while W, is not a true measure of the adsorp- 
tion, an approximate idea of the changes taking place can be obtained from it. 
Considering the curves in Fig. 4 to be of the form of the true adsorption 
curve, the hysteresis can be completely accounted for on the basis of the 
hypothesis stated earlier. There is a persistence of the liquid structure above 
the critical temperature, and we may conveniently call this the ‘‘orientation 
effect”. An accurate measure of how this decreases with temperature has 
been obtained by Tapp, and Fig. 9 has been reproduced to show how the 
density difference above and below 
the region of meniscus disappearance 
decreases with temperature. This 
density difference may be used as a 
measure of the orientation remain- 
ing at one temperature. Hence 
when the material in the region 
previously occupied by liquid is 
~ heated to 135° C. the orientation has 
largely disappeared, and since, on 
cooling, it does not reappear until 
the critical temperature has been 
passed, the material may be regarded 
as gas-like in properties—i.e., the 
molecular motion is more chaotic on 
cooling than on heating. Since ad- 
sorption from the true liquid is zero, 
and that from the vapor or gas is 
large, we must expect a hysteresis in tine a—_, —. ee 
curve showing greater adsorption. meniscus disappearance and the temperature. 

A second consideration is the den- Reproduced from Tapp, Steacie, and Maass 
sity of the medium. The adsorption os 
will be proportional to the concentration of sorbate molecules, and in the inter- 
mediate temperature range of any of the curves in Fig. 4, this concentration 
is lower during cooling than upon heating, so that a decreased adsorption 
must be expected. The fact that the hysteresis curve lies above the solid 
curve shows the influence of the orientation effect, while the intermediate 
region of lower adsorption on cooling shows that the density effect is greater 
than the orientation effect. 

In the experiments from which the results shown in Fig. 5 were obtained, 
no liquid-to-gas change was brought about, and the hysteresis curve lies above 
the plain curve throughout the complete temperature range. This is due to 
the fact that the concentration of sorbate during the cooling part of the 
experiment was greater than it was before. 


DENSITY DIFFERENCE 
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A final comment may be made on the form of the curves (Figs. 2 and 4) in 
the region of temperature equalization between the two arms. From the 
inflection in the density curve and the abrupt increase in adsorption it appears 
that the decrease in liquid structure becomes very great in this region. Why 
this should be is difficult to explain, but it is highly probable that the expansion 
of material from one arm to the other, when the temperature of the first passes 
through that of the second, greatly influences the destruction of residual 
orientation existing in the material. 


Summary 


1. The difference in density above and below the region of disappearance of 
the meniscus observed by other investigators has been confirmed both in 
nature and magnitude by an entirely different experimental procedure. 


2. Adsorption studies in the region of the critical temperature have indicated 
that a discontinuity in adsorption accompanies the transition of liquid to gas. 


3. It has been demonstrated that the apparent changes in adsorption cannot 
be ascribed to changes in the density of the adsorbed phase or to persistence of 
liquid in the pores of the sorbent above the critical temperature, and further- 


more, there is no evidence of a lag in critical temperature existing for such 
material. 


4. The following hypothesis to account for the phenomena observed is 
proposed :— “‘In the liquid state a certain amount of molecular orientation 
exists. The nature of this orientation may be pictured by considering a 
sphere of radius equal to the average distance between the centres of molecules 
in any part of the liquid. The molecules in such a sphere have on the average 
a definite orientation relative to one another. If a radius twice as great is 
chosen the molecules may still have a resultant mean orientation but of 
smaller magnitude; thus with increased radius the average orientation rapidly 
diminishes. A liquid has therefore a tendency to form a structure into and 
out of which the molecules move. Increased temperature decreases the 
number momentarily in the structure, and at the critical temperature this 
structure no longer corresponds to the minimum potential energy of the system 
due to heat energy (translational, rotational, and vibrational) of the molecules. 
This structure has a tendency to persist above the critical temperature,and 
in a sense may be regarded as metastable, though efficient mechanical agitation 
has no appreciable effect. A discontinuity in properties other than that due 
to concentration exists at the critical temperature. Thus the total surface 
energy is not zero (15), the velocity of a chemié€al reaction may undergo a 
discontinuous change (12), and the adsorption from liquid and that from gas 
may be very different’’. 


The conception of the liquid state given above is not a new one (5, 14), 
though it has been held that the structure should become zero at the critical 
temperature. Contributions to this subject have been made by Stewart 
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(11) through a study of the X-ray diffraction caused by liquids, and the 
extension of this field to include critical conditions should lead to interesting 
results. 

Preliminary observations on meniscus behavior have indicated that the 
temperature and position of reappearance of a meniscus is dependent on the 
previous history of the material contained in the region previously occupied 
by liquid. Further experimentation along these lines has yet to be completed, 
but it appears that a meniscus reappears at a higher po-ition, and at a tem- 
perature less below the critical, when the contents of the bomb are not pre- 
viously heated appreciably above the critical temperature, so that the ease 
of formation of liquid is proportional to the residual orientation. 
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THE THERMAL CONDUCTIVITY OF DEUTERIUM! 


By A. B. Van CLEAVE? AND O. Maass? 


Abstract 


The thermal conductivities of deuterium and some mixtures of deuterium 
and hydrogen have been measured by a relative, “hot wire’ method. The 
results are consistent with the authors’ original conclusion that the deuterium 
molecule has the same molecular diameter as the hydrogen molecule. It follows 
also that the molecular heats of the hydrogen isotopes are the same. 


Introduction 


In a recent publication (7) the authors gave results for the coefficient of 
viscosity of deuterium and deuterium-hydrogen mixtures which indicated 
that the hydrogen isotopes have the same molecular diameter. The work 
described here was carried out to obtain additional support for the foregoing 
observation. 


Farkas and Farkas (2) have performed experiments on the thermal con- 
ductivity of deuterium-hydrogen mixtures, but they were primarily interested 
in devising a method of analysis for mixtures of the isotopes, and their results 
cannot be used to make calculations of the thermal conductivity of deuterium. 


Experimental 


Owing to the large number of somewhat uncertain corrections to be applied, 
measurements of the absolute thermal conductivity of gases are very incon- 
sistent and difficult to obtain. For this reason a relative method was em- 
ployed. 

The apparatus was modeled after that used by Palmer and Weaver (4). 
A platinum wire 0.07 mm. in diameter and about 10 cm. in length was 
stretched along the axis of a cylindrical brass tube of 4 mm. radius. The wire 
was kept under tension by means of a helical spring of copper wire. One end 
of the wire was electrically insulated from the metal chamber by attaching it 
to the terminal of an ordinary spark plug which had been fitted on to the end 
of the brass cylinder. The cell was provided with an inlet tube which could 
be connected to the deuterium container, air, hydrogen generator or pumps 
as desired. It was mounted vertically in an oil bath maintained at 25.00+ 
0.10° C. by means of a thermostat. 

The electrical arrangement took the form of a Wheatstone bridge. The 
thermal conductivity cell and an ammeter formed one arm of the bridge. 
The corresponding arm, which consisted of a coil of copper wire of about 
one ohm resistance, was also immersed in the thermostat. Two variable, 
high resistances made up the remaining arms of the bridge. One of these 

1 Manuscript received February 16, 1935. 

Contribution from the Physical Chemistry Laboratory, McGill University, Montreal, Canada. 
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resistance boxes could be varied to within 0.1 ohm. The current was supplied 
by a bank of dry cell batteries and could be varied by means of resistances of 
66 ohms and 1 ohm connected in series. 


Assuming that the loss of heat by radiation, by conduction through the 
connections to the ends of the wire, and by thermal convection is constant, 
the relative amounts of power expended in maintaining a given temperature 
difference between the wire and the walls of the chamber, when filled with 
various gases, affords a measure of the relative thermal conductivities of 
the gases. 


The heat loss, g, from an electrically heated wire is given by 
g = Jit, 


where J is a constant of proportionality, 7 the current, 7 the resistance, and: 
the time of flow. Since in any one experiment the temperature, and there- 
fore the resistance, of the wire is kept constant, the only quantity that need 
be measured is the current. The amount of heat given out to the thermostat 
in unit time is directly proportional to the thermal conductivity of the gas 
between the wire and the metal tube, and to the square of the current which 
must be used to keep the temperature of the platinum wire at the same value 
for all the gases concerned. 


If the thermal conductivity of hydrogen be taken as unity, that of any other 
gas may be found by measuring the currents which will maintain the platinum 
wire at the same temperature, first in hydrogen and then in the gas of unknown 
conductivity. That is, 


Pras , 
hen = ba 


where ky is the thermal conductivity of hydrogen, ig the current in hydrogen, 
teas the current and ky, the thermal conductivity of the gas under con- 
sideration. 


The cell was first filled with pure dry hydrogen. With the current and one 
of the variable resistances kept at definite values, the other variable resist- 
ance was adjusted until the bridge was balanced. After noting the values of 
the current and resistance for the balanced condition the first resistance was 
reset and the procedure repeated. The current was changed slightly and the 
process again repeated. The apparatus was then filled with deuterium (98% 
by weight) and, with the variable resistance arms at the same values as found 
for hydrogen, the current through the wire was varied until the bridge was 
balanced, thus indicating that the wire had assumed exactly the same resist- 
ance as it had in hydrogen. This procedure was repeated for each of the 
deuterium-hydrogen mixtures. These mixtures were prepared as described in 
the previous paper (7). ; 

The ammeter which measured the current flowing through the wire could 
be read to +1 milliampere, and was found to give relatively correct values 
over the range of current used. 


For all measurements the gas pressure was one atmosphere. 
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As a test of the accuracy obtainable the relative conductivities of hydrogen 
and air were measured and found to agree within 1% with values recently 
obtained by the above method. 


Results 


Table I is a summary of the currents necessary to maintain the platinum 
wire at the temperature corresponding to current readings of 0.2680, 0.2900, 
0.3175 and 0.3430 amp. in pure hydrogen. Each of these values is the 
mean of four independent settings of the variable resistance arms. 


TABLE I 
VARIATION OF HEATING CURRENT WITH CONCENTRATION 


Weight of hydrogen, % 


100 | 57.05 | 32.90 14.06 | 7.74 | 2.00 








Current, amp. 


0.2680 0.2571 0.2465 0.2355 0.2307 0.2280 





Table II shows the ratio between the thermal conductivities of various 
deuterium-hydrogen mixtures and hydrogen as determined from the relation 





Le 
ka 1h 
TABLE II 
YS 
VARIATION OF + WITH CONCENTRATION 
Weight of hydrogen, % 
100 | 57.05 32.90 | 14.06 | 7.74 | 2.00 
00/th 
1.0000 0.9225 0.8460 0.7722 
1.0000 0.9276 0.8477 0.7762 
1.0000 0.9108 0.8372 0.7722 
1.0000 0.9261 0.8434 0.7701 


Mean 1.0000 0.9216 0.7727 








Within the error of this determination the ratio ofthe thermal conductivities 
of deuterium-hydrogen mixtures and hydrogen is the same for the four tem- 
peratures at which measurements were made. These results indicate that 
the temperature coefficients of the thermal conductivity of deuterium and 
hydrogen are the same, as would be expected. The following calculations 
have been made on this basis. 


ee 
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Discussion 


Before discussing the foregoing results it is desirable to refer to some of the 
conclusions drawn in the previous paper (7) in regard to deuterium-hydrogen 
mixtures. It was pointed out that complete agreement between experimental 
results and calculations based on theoretical formulas for the viscosity of gas 
mixtures was not found, and that a complete discussion of this would be left 
until a purer sample of deuterium had been obtained. The necessity for this 
is much more apparent when it is considered that HD molecules are formed 
during the preparation of D. from 98% D.O, thus decreasing the molecular 
percentage of D,. A calculation* shows that the molecular percentages of 
D, and HD in the original sample were 92.54 and 7.06 respectively. This 
does not alter the conclusions reached with regard to the value of the viscosity 
of pure deuterium. 

If it is assumed that HD molecules have the same molecular diameter as 
Hz molecules, the viscosity of pure HD is 1.225 times that of pure H,. The 
original sample, which consisted almost entirely of D, with HD molecules 
rather than H, molecules as an impurity, was found to have a viscosity of 
1.400 referred to the viscosity of hydrogen as unity. A straight line extra- 
polation for these values of 100% HD and 7.06% HD gives the relative 
viscosity of 100% D,z as 1.414+0.003, which is exactly the same as that 
previously reported (7). 

This result may also be obtained by using Puluj’s relation (5) for the vis- 
cosity of a gas mixture as an interpolation formula, the value being 1.415+ 
0.003. 

It might be argued that since the HD molecule is asymmetrical it would 
tend to have an effectively larger collisiona area than Hz. Assuming this 
to be the case, it would of necessity follow from the extrapolations that the 
D, molecule has a smaller molecular volume than Hz. This does not seem 
possible. It follows, therefore, that the previous results (7) indicate that 
H., HD and D, molecules have the same molecular diameter to within 0.20%. 

The thermal conductivity measurements lead to the same conclusion. 
From the well known Maxwellian expression 

ke = enol, 

where k, is the thermal conductivity of a gas at 0° C., », is the corresponding 
coefficient of viscosity, C, is the specific heat at constant volume, and e« is 
a constant, it is possible to estimate k, for D, and HD. Reasons have been 
given above for believing that D., HD and He molecules have the same 
molecular diameter, hence their viscosities will be in the ratio 1.414 : 1.225 
: 1.000. According to Farkas and Farkas (2) the molecular heats of Ds, 
HD and H; are the same at temperatures above 0° C., hence their specific 
heats will be in the ratio 0.500 : 0.666 : 1.000. For isotopes the constant 
¢ should be the same. Hence, the ratio of the thermal conductivities of 
D., HD and Hz should be 0.707 : 0.817 : 1.000. 


* The authors are indebted to Dr. H. C. Torrey of Columbia University for bringing this 
calculation to their attention. 
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The experimental value for the relative thermal conductivity of the mixture 
containing 92.54 molecular per cent of D, and 7.06 molecular per cent of HD 
is 0.7165. On the hypothesis mentioned above, the thermal conductivity 
of HD is 0.817 referred to Hz as unity. A straight line extrapolation of 
these two values gives the relative thermal conductivity of pure D, as 
0.708+0.007, which checks exceedingly well with the theoretical value 
deduced above. Since the thermal conductivities of D, and HD are certainly 
not very different (2), the error in this extrapolation will be well within the 
experimental error. Since these results fit so nicely with the results on vis- 
cosity the authors believe that this is additional proof that D., HD and H, 
molecules have the same molecular diameters and mean free paths at 0° C. 


The absolute value of the thermal conductivity of hydrogen at 0° C. is 
still a somewhat uncertain quantity. The mean of what are, perhaps, four of 
the most careful determinations (1, 3, 6, 8) with the ‘‘hot wire’’ method gives 
k. for hydrogen as 0.000416 cal. cm.—' sec! deg! C. The maximum 
deviation from the mean is 0.72%. Relative to the above value of &, for 
hydrogen, the authors’ value of k, for deuterium becomes 0.000295 + 0.000003 
cal. cm.—! sec.-! deg. C. The calculated value of k, for HD is 0.000340 
cal. cm.—! sec.—! deg. —! C. 

Thermal conductivity measurements of other deuterium gas mixtures are 
contemplated. 
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MEASUREMENT OF DIELECTRIC CONSTANT OF 
ELECTROLYTES' 


By R. W. McKay? 


Abstract 


An apparatus for measuring the dielectric constant of conducting solutions 
by voltage resonance is described. The theory of the circuit and the sources of 
error are discussed. The design of the apparatus is such as to eliminate errors 
other than those due to inductance of cell leads, for which a correction is made. 

Measurements have been made on solutions of sodium chloride, hydrogen 
chloride, potassium sulphate, magnesium sulphate, copper sulphate and potas- 
sium ferrocyanide, at 23.0° C. and 2 X 10® cycles frequency. An increase in 
dielectric constant with concentration greater than that predicted by Debye 
and Falkenhagen has been found in all cases. The results are compared with 
those of other workers. Graphs are given of the results obtained. 


1. Introduction 


1.1 In recent years striking advances have been made in the theory of 
electrolytic solutions. Assuming complete dissociation for dilute solutions 
of strong electrolytes, Debye and others have calculated approximately, not 
only the thermodynamic functions for equilibrium states, but also such 
characteristic quantities for irreversible processes, as viscosity, conductivity 
and dielectric constant. The experimental data available are not sufficient 
to test the theory in all its phases. This is especially true of dielectric constant 
measurements. Although nearly 50 years have elapsed since Cohn and Arons 
made the first measurements of dielectric constant of salt solutions, great 
discrepancies still exist among the results obtained by many observers, and 
even the most recent measurements are uncertain. The experimental diffi- 
culties which have prevented accurate measurements have been discussed 
by many authors and complete summaries of the earlier work have been 
written by Walden and Ulich (19), Lichtenecker (13) and Blith (1). An 
extensive bibliography is given by Lattey (12), so that only the most recent 
work is referred to here. 


1.2 The technique of measurement of dielectric constants of insulating 
materials is relatively simple and several precise methods have been developed 
for work with such substances. When an attempt is made to extend these 
methods to conducting materials such as electrolytic solutions, many difficulties 
arise. Consequently the earlier papers in this field do not contain reliable 
results, although they give much information which is useful in developing 


1 Manuscript received October 18, 1934. 
Contribution from the Department of Physics, University of Toronto, Toronto, Canada. 
This paper was submitted in partial fulfilment of the requirements for the degree of Ph.D. at the 
University of Toronto. 


2 Holder of a studentship under the National Research Council of Canada. 











378 CANADIAN JOURNAL OF RESEARCH 


more accurate methods. The inconsistency of the results obtained is shown 
in Figs. 1 and 2, which are plotted from measurements for sodium chloride 
and copper sulphate by many authors. Some of the more recent papers, how- 
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CONCENTRATION  MULLIEQUIVALENTS PER LITRE 


Fic. 1. Measurements for univalent salts by various 
authors. (1) Millicka and Slama, 50 cycles. (2) Schmidt, 
500 cycles. (3) Smale, 50 cycles (KCl). (4) Pechold, 
50 cycles. (5S) Drude, 1 metre; Voigt, 40 cm.; Plétze, 
60 cm.; Deubner, 130 cm. (6) Malone, Case and 
Ferguson, 31 metres (1921). (7) Lattey. (8) Walden, 
Ulich and Werner. (9) Astin, 10° cycles. (10) Nernst, 
10° cycles. (11) Devoto. (12) Lattey and Davies (1931) 
4-7.5X10' cycles. (13) Sack, 2X10® cycles. (14) 
Skancke and Schreiner, 10° cycles. (15) Jezewski and 
Kamecki, 5 X10* cycles. 
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CONCENTRATION MIULLIEQUIVALENTS PER LITER 
. 2. ( Sree for divalent salts by various Drude, adapted for very 

authors. (1 midi, 500 cycles. (2) Smale, 50 cycles. i i 

(3) Drude, 1 metre; Plitse, 0.40 metre; Voigt, 0.60 high concentrations, finds 

metre. (4) Malone, Case and Ferguson, 3 metres. (5) 


Lattey and Davies, 10° cycles (1932). (6) Deubner, 
135 cm. (7) Wien, 7.5X10" — (MgS0,). (8) 
Wenk 1.1X10@ cycles. (9) Jesewski 
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ever, show much better 
agreement and appear more 
reliable. These researches 
are discussed briefly below. 

1.3 In the frequency 


yi —T_] 
Za ea ® ar 3 re 10° cycles per 
SES sec. the most consistent 
wy etal) a aan results have been obtain- 


ed. Using the original 
method of Drude, with a 
Barkhausen oscillator sub- 
stituted for the Hertzian 
resonator, Plétze (17) finds 
that there is no change 
in dielectric constant for 
copper sulphate up to 2 
millinormal within +0.3%, 
nor for potassium ferro- 
cyanide within +0.6% for 
the same range of con- 
centrations at 0.60 metres 
wave-length (1.67 xX 108 
cycles). Drake, Pierce and 
Dow (4) have modified the 
method of Drude by sub- 
stituting a concentric tube 
and wire for the parallel 
wire system. In agreement 
with Plétze they find no 
change for potassium chlo- 
ride up to 7 millinormal 
within +1.0%. The work 
of Malone, Case and Fer- 
guson (14, 15) is not in 
agreement with this, but 
does not appear to be re- 
liable. Stips (18), using 
the second method of 


that the dielectric con- 
stants for normal solutions 
of copper sulphate and 
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magnesium sulphate increase linearly with concentration up to 10% increase. 
This is not in disagreement with the results of Plétze or Drake, Pierce 
and Dow, as the value for small concentrations would be well within the 
limits of error stated by these authors. Using Wien’s ‘“‘Barretter’’ method at 
from 1.5 X 10* to 0.6 X 108 cycles per sec., Wenk (20) and Dobenecker (3) 
find a slight increase in dielectric constant for the salts they have measured. 
Wien (21) used the ‘‘Barretter’’ method at 3 X 107,1.5 X 10’ and 1.0 X 107 
cycles frequency and found an increase for all salts. The increase becomes 
larger the longer the wave-length. All these results, except those of Malone, 
Case, and Ferguson, are substantially in agreement with the theoretical 
predictions of Debye and Falkenhagen (2). 


1.4 For somewhat lower frequencies the measurements do not appear so 
certain. Lattey and Gatty (12) and Lattey and Davies (10, 11) working at 
10 ke., and Jezewski and Kamecki (9) at 5 kc., obtained widely divergent 
results although the work in both cases was carefully done. The work de- 
scribed in this paper gives results more nearly in agreement with those of 
Lattey. These are discussed in Section 5. 


Graffunder and Weber (6) have also made measurements in this frequency 
region but the method used is suited only to extremely low concentrations, 
and even for these it seems doubtful whether any weight can be given to 
their results. Their method depends on changing the frequency of an-oscilla- 
tor by placing the measuring cell in the oscillating circuit. This method is 
extremely sensitive but not very accurate. 


1.5 Many measurements have been made at very low frequencies using 
the Force method. Orthmann (16) gives a summary of work by this method 
up to 1931. In his own work Orthmann used a silver ellipsoid suspended 
in solutions of silver salts. In this way he believes he has avoided polarization 
effects. Even in such cases, however, ‘‘over-voltage” effects may cause 
large errors. A. R. Gordon* has shown in some unpublished work that for 
copper electrodes in copper sulphate solutions, overvoltage exists and is out 
ot phase with the applied voltage at low frequencies. If the same is true for 
silver in silver salts, errors from this cause would occur in Orthmann’s work. 


2. Apparatus, Design and Construction 


2.1 As the agreement which has been obtained among the results of various 
workers for the highest frequencies is fairly satisfactory, a somewhat lower 
frequency (2 X 10° cycles) was chosen for this research. At this frequency 
the choice of methods is practically limited to the bridge method or the 
voltage resonance method. Both have been used in this work, but no results 
of measurements with the bridge are recorded because the resonance method 
has been found to be much more satisfactory. 


*Prof. A. R. Gordon of the Department of Chemistry, University of Toronto. 








380 CANADIAN JOURNAL OF RESEARCH 


2.2 Bridge Method 


The advantages of a null method are well known. In the case of high- 
frequency capacity bridges, however, there are also many disadvantages, 
the most important of which are listed below :— 


1. Insensitivity. Since almost all alternating current detectors give read- 
ings proportional to the square of the input voltage, the sensitivity becomes 
small as zero reading is reached, making it difficult to balance the bridge 
accurately. By using an amplifier of sufficient sensitivity it is possible to 
overcome this difficulty. 


2. Difficulties in shielding. ‘These are common to all high-frequency 
measurements, but are somewhat accentuated in bridges because the circuit 
is more complex than, for instance, in resonance measurements. 


3. Calibration of phase-angle of variable resistances. Resistances can be 
made with very small phase angle, but this small remaining phase angle is 
still sufficient to produce errors in the measurement of capacity of cells con- 
taining electrolyte. Various “‘reactionless’’ resistances have been devised 
(see Schreiner, Sack, Weber) but some doubt exists as to the degree to which 
the reactance has been removed. 


These difficulties may be minimized by choosing a suitable bridge circuit. 
After a careful comparison of numerous bridge circuits, the design shown in 
Fig. 3 was chosen for preliminary work. It is essentially similar to the Scher- 

ing bridge used com- 
mercially for power factor 
determinations. It differs 
chiefly from the standard 
type in having no variable 
resistors in the measuring 
circuit. This change great- 
ly simplifies calibration 


aa of the bridge. The use 
fetes Se of a Wagner ground arm 
AM PUFIER ° . 

pant makes it possible to keep 

Fic. 3. Bridge circuit diagram. one corner of the bridge at 

ground potential without connecting it directly to ground. This permits the 

use of a three-electrode cell (Fig. 4A) so that the field, in the part of the cell 

which is measured, is entirely in the electrolyte and not through any sup- 
porting insulators. 


Condensers were used as far as possible for the various units of the bridge 
arms, since they are more nearly ideal for phase angle than resistances or 
inductances, and they are also more easily shielded. The theory of the bridge 
is well known. The value of the capacity in the measuring arm is given by 
the formula 


DE 
iter 
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where C is the total capacity of the measuring arm, E that of the comparison 
arm, and D and F the capacities of the condensers in the other two arms. 
C is independent of frequency, which is an advantage of this bridge in common 
with all those of the so-called ‘‘perfect balance’’ type. 


@ 








() 
Rubber rminal 
Greunded Rubber 
Shield 
Mevable 
Inside Eladirede Ourside Electrode Discs 
Cat A 
() 
@) 
Glass 
msc E> wil 
Celt C WIEN'S CELL 


Fic. 4. Cells. 


The variable condensers E and F each consist of fixed condensers in parallel 
with a small variable condenser. Their capacity is variable from Syuyf. to 150 
puf. and can be read at any point to within 0.02 uuf. They were mounted 
rigidly in a brass box and shielded from each other. The input from the 
generator may be varied from zero up to one watt. The detector consists 
of a three-stage screen-grid amplifier, and a rectifier. The sensitivity is 
sufficient to balance the bridge within 0.1% when the cell resistance is as low 
as 100 ohms. 


The only point causing trouble in calibration was the evaluation of the 
shunt capacity of the resistance Re. The measurement of any resistance of 
known phase angle suffices to evaluate the capacity of Rg. Various methods 
tried to find such a resistance are summarized below. 


1. Metal films sputtered on glass cylinders might be supposed to have no 
capacity apart from that of their terminals. This was shown definitely to be 
untrue for carbon and carborundum in the same geometrical form. A pencil 
mark on hard rubber was shown to have an effective capacity as high as 
10 ywuf. This is probably an anomalous behavior due to the structure of 
carbon, but one cannot be certain that the same effect is not found to a lesser 
degree in metal films also. This method cannot, therefore, be relied on. 
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2. Resistance can be measured accurately on the bridge. If a circuit consist- 
ing of resistance, capacity and inductance is adjusted in the bridge till its 
effective resistance becomes a minimum, it may be shown that for certain 
circuits the reactance is then zero. Such a circuit provides a method for 
calibrating Rg. The sensitivity of the method is not great. 

3. A shielded resistance was sent to the Bell Telephone Laboratories Inc. 
(New York) who were kind enough to offer to calibrate it. The calibrated 


value was 181.13 ohms and —1.3 + .5 pyf. at 400 kc. Using this resistor 
the bridge calibration was completed. 


The tests which were made with the bridge, the errors in measurement, 
and the reasons why the bridge was replaced by a resonance circuit are 
described below in Section 3 (1). 


2.3 Resonance Method ° 


2.31 Theory. The theory of the voltage resonance circuit has been given 
by Jezewski and by Astin for the case of an isolated resonator, without con- 
sidering the interaction with the generator and the detector. Lattey con- 
sidered the interaction with the detector only. Both these interactions 
may be large enough to affect the results. It is therefore worth while to con- 
consider the theory of the method in greater detail. 

After several other circuits had been tested, the one shown in Fig. 5 was 
adopted. A quartz crystal oscillator not closely connected to the rest of the 


= 
3 “2 rey 4135 


FIRST AMPLIFIER RESONANCE CIRCUIT 


SECOND AMPLIFIER 
Fic. 5. Resonance apparatus—circuit diagram. 


circuit supplied a high-frequency current of constant frequency. This current 
was used to drive a 32-type screen-grid valve, the plate circuit of which was 
coupled to the coil of the resonance circuit. The resonance circuit was con- 
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nected to a second 32-type valve which had an impedance of only 1000 ohms 
in the plate circuit, minimizing the possibility of feed back. These two 
valves, V2 and V3, effectively isolate the resonance circuit from the oscillator 
or detector, so that for the purpose of calculation the equivalent circuit Fig. 6 
may be substituted. The 

valves V2 and V3 are re- 2 Core 
placed in this diagram by 
their equivalent networks, t 
in the conventional way. 
The plate-filament and grid- 
filament capacities are in- 
cluded in the units /];, Z;, C 
and /,. It is necessary to 
consider the grid-plate capa- 
cities, so they are represent- 
ed as Cyy, and Cy,, in the dia- 
gram. The plate resistances Fic. 6. Resonance apparatus—equivalent circuit for 
and amplification constants  “““/aHons. 

of the tubes are represented by 1, p2 and si, we. It is necessary to assume 
that the voltage e is harmonic and that the frequency is constant. The experi- 
mental justification of these assumptions is given below. In Fig. 6, Cy repre- 
sents the capacity between the resonance coil L and the plate-circuit coil L; 
between which the mutual inductance is M. The cell capacity is C, and R 
includes the shunt resistance of the cell as well as the grid resistance of Vs. 





at 


Kirchoff’s equations for the steady state may be written at once in the 


following form. (The currents i, . . . . i, are all taken in the clockwise 
sense). , 
(jak + re atz)ie = whe =¢—méa (1) 
— pita + (pr + jwl:) in — jwo(Li + M) i + joMia = ply (2) 
— joLs + Mis + jo( Ls + L + 2M - 3h )ie - jl + Miz = 0 (3) 
—joMis — jo(L + Mic + (job tr - te) t+ Si wi (4) 
ee a , J jJwC. ” a 
Ziut(r-F i, — Riy =0 (5) 
—Ris + (R+ m— =o )ir - orn = -men (6) 
— priz + (a2 + jule + ro)ig ¥ = pals (7) 
Substitute 


tn = tars + jo) — 6 
Cn = (if — 4.)R 
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Divide Equation (1) by (1+ 4:) giving Equation (1') 
Equation (1) Xm 
Add Stes 
Add Equation (4) and Equation (5) giving Equation (4') 
Add Equation (6) and Equation (7) giving Equation (6') 
Replace the two variables 7, and i; by new variables 
(ty—t.) and (t,—ia). 


and (2) giving Equation (2') 


This gives 
1 
; Foleo, +"), — 
Fe at fa 


Tate it (jus +— >is — jolLi + Mic + joi, 

— joM(i, — ie)” 
—ju(Li + Mis + jo(Li + L + 2M - 37) + ju(L + Min 

+ ju(L + Mi. — 64) 

joMir — jo(L + Mic — Rliy — i.) + (jul + rie 

— (jul + rie — a) 
(R +3 =o) - te) +2 = te + (jul. + ro)tg 
(- an — pr) (iy — i.) — prise + (p2 + jul, + To)tg 
— Rliy = i) + elie — 50) 


Let A_ be the determinant of the coefficients of i., etc. 
Ax be the minor row ‘‘s”’ and column “k”’ of A 
Dj, be the minor of row “i’’ and column ‘“‘k’’ of Ajz 


~ 
Air po Du _ RDu oe - a 
nt 
Du = 0 identically . 
An _ 2 


So Zi = jwCR De 


where A+ represents sa conjugate complex a A 


An Ak Arn — An Oe , | Arn! 
‘I tpearn—— + Des 


e A 


- =Zy = 


= wor | 


which is quadratic in C. 


For any value C, of C another value C; will give the same |Z;7|? and there- 
fore the same output current. Since the expression is quadratic in C the 
mean of C,; and C; will be the same for any corresponding pair of values. 


~ ata. vey ax dR( An Ay — Az An) 
3 resonance capacity ZR? | And? 


oe An 


= the imaginary part of oRAx 
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Rearranging the columns slightly and multiplying the fifth row of each deter- 
minant by (p,jwC,,,) and adding the sixth row to it, the two determinants 
can be written :— 


-An _ 


al 
mt+i 


jwCor, 
m+i1 


° — jo(l, + M) 


jol; — jo(L, + M) 


eal 
mti 


jo(La +L+2M -a) 


° 
oO 
O° 


+ jwoM 
Oo 
oO 


— jo(L + M) 
oO 
° 


Right half of above-———— 
Oo 
joM 
— jo(L + M) 
jol +r oO 1 
pr + (jol, + 17.)(1 + joCor or) 
pr + jule + re m+ 


Mm — jwCo Ps 


confi 


m+ 
jel, + —“— 


mtil 
— jw(Ii + M) jolla + L + 2M — Se 
joM — jo(L + M) 
Oo 
O° 


— jo(li + M) 
1 


Right half of above——— 
oO 
joM 
— jo(L + M) 
joL +r ° 
oO pr + (jules + r)(1 + joCon, ps) 
—~ aa + jul, +1. 
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Let dig be the minor of the ith row and sixth column in either determinant 
Then 


a dig — (uz — jwCom pr) dss + (« + #) dee 
C = imaginary part of ———————@—_—$_ _ _<_ 
— wor dee 


— du + (us ~_ jwCov pr) ds6 


— ‘ — 
imaginary part o candies 


None of these terms are functions of R, so the important property of the 
circuit has been proved: 


The resonance curve is symmetrical and the resonance capacity is independent 
of the shunt resistance R. 


In the above it has been assumed that » and pare constants, and that r;, 1;, ro 
and |, do not depend on the voltage. These assumptions are only approximate, 
and it can be shown that in certain cases variation of these quantities may 
introduce errors in the results. To determine the magnitude of these errors 
the expression for C is worked out. 


E ie Oe’) .. imact bu 
C= Cup X (1 + 2) imaginary part of ob 


in which 
6= 


1 
ae 
dle ty > eg mnifite 


1+ im mt il 
4 1 ; ‘ . 
joli + 75 + =r. joli — jo(Li + M) jwoM 
. P 1 
o —ju(Li + M) jo L + Li + 2M -~37-) jell +) 
Oo jwoM — jw(L + M) jol +r 


and dy is the minor of the term in the lower right corner. 


Using the numerical values 
a“ =6xX 16 jol; + 7; = 10,000 ohms (non-reactive) 
wCon, 
p = 16 jol 500 ohms 


uw = 500 jeM = 100 ohms 


Thus jwl, + 1; is negligible in comparison with 
wCop, 

jwL, is negligible in comparison with 
wCop, 


If we neglect also Cy, we get simply 


“oa Ss 
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2.32.Construction. The construction of the resonance apparatus is shown 
diagrammatically in Fig. 7. The oscillator and first amplifier are enclosed in 
separate copper cases, and the battery leads and batteries are also shielded. 
The coil in the output 
circuit of the amplifier is 
wound on the same form 
as the resonance coil L. 
The theory shows that this 
close coupling is not harm- 
ful. On the contrary, it 
makes shielding much sim- 
pler and permits the use 
of an oscillator of lower 
power than is generally 
used. The capacity of the 
resonance circuit is made 
up of four parts: two 
uncalibrated variable con- 
densers, each 350 uuf. maximum, one cylindrical condenser 15 puf. maximum 
and calibrated to read to 0.003 uuf., and the cell containing the electrolyte. 
This part of the apparatus was built into a rigid box of brass and the leads 
to and from it were fixed in tubular shields by means of supporting insulators. 





In finding the resonance point, the condensers are adjusted so that the 
current in the galvanometer is zero for both positions of S (Fig. 5). The 
resonance capacity is then the average of the values when the switch is open 
and closed (since the curve has been shown to be symmetrical). The water 
in the cell is then replaced by solution and the small cylindrical condenser 
is varied until the “balance” is again obtained. The difference in capacity 
is then read on the dial of the small condenser. This method is much more 
sensitive than that of adjusting the condensers for maximum reading of the 
detector. 

The amplifier and detector circuits were enclosed in separate shields and 
their batteries were shielded separately from those of the oscillator. 

Voltmeters were permanently connected to the filament batteries and the 
rheostats kept adjusted to give the rated filament voltage. This was neces- 
sary since variations of filament currents produced small but noticeable 
errors in measurement. A change of 1% in voltage produced a capacity error 
less than 0.03 yyf. 


3. Preliminary Tests and Discussion of Errors 
3.1 Bridge Method 


After a series of tests with the bridge described it became apparent 
that though the shielding had been very carefully constructed, coupling 
between various parts ot the circuit still introduced errors for which no 
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correction could be made. The results of measurements on sodium chloride 
are roughly in accord with the resonance measurements given below, but are 
not considered sufficiently accurate to publish. The bridge is described, 
however, as it was found to be satisfactory for comparison measurements 
and might therefore be useful in later work. 


3.2 Resonance Method 


3.21 Harmonics and frequency variations. The theoretical discussion 
assumes a pure sine wave of constant frequency. A distorted wave form 
produces very small errors in resonance measurements with circuits having 
small dissipation. In the circuits used here, however, the effect of harmonics 

paling napa ongeona energie ys is much greater owing to 

ret the fact that the curves 

:, Pee a. ee are so flat that the curves 
Som | ak for the harmonics overlap 
on the curve for the funda- 

mental. Fig. 8 illustrates 

the effect of harmonics. 

By inserting the tuned cir- 

cuits (J; and 7; in Fig. 5) 

and taking care not to 

overload valve V2, the 

error due to harmonics is 

removed. The test for sym- 

metry of the curve, given 

below, shows that this is 


Pee. 6. Lies of harmenies. realized experimentally. 


To detect frequency changes a vacuum-tube wavemeter circuit was set up 
and coupled loosely to the oscillator. Operating switch S (Fig. 5) with the 
cell in or with the cell out produced no observable change of frequency. 
Neither did insertion or removal of the cell filled with either water or solution. 
It was concluded that under the conditions of operation frequency changes 
were always less than 10 cycles per second at 2 X 10* cycles per second. 


Small variations in frequency occurred immediately after the currents 
had first been switched on, but ceased within an hour. 


3.22 Coupling of oscillator and detector. Coupling between the detector 
and other circuits produces asymmetry of the resonance curve and hence 
errors in the results. This coupling was reduced to negligible magnitude by 
careful shielding. With the cell short-circuited no movement of the galvano- 
meter could be detected when the generator was turned on and off. The 
symmetry of the resonance curve is additional evidence that this error has 
been removed. 
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Symmetry test. The TABLE | 

capacity of a resistance Reading | Reading Diff 

Ci, with with ——— 
of 180 ohms was measured oak resistance | resistance parallel = Cagney 
for various values of Ci, out in oF senator 
the condenser in series —— 
with switch § (Fig. 5) 350 13.2 20.6 + probable erat 
corresponding to different . "5 = 0.2 div. 
heightsonthecurve. These 195 10.5 17.7 og Oe 


measurements are given in 
Table I. 
3.23 Changes in effective tube capacities. In section 2.31 the formula found 
for the resonance capacity is 
prs 1 ua To 
C 3 ee Coed 1 + 
aay 
The constants of the first valve do not occur in this approximate formula, 
so that their variations are unimportant. The constants of the second valve, 
however, are present in the formula and their effect may be regarded as a 
modification of the effective input capacity of the valve. The term r, includes 
the impedance, Z, of the coupling unit between V3; and Vyin Fig. 5, and the 
input impedance of the next valve, V4. The input impedance of this valve 
will depend on applied voltage, and therefore on R. For large values of 
R (i.e., for dilute solutions), 7. will be reduced so that the reading for C will 
be too large. This was verified experimentally with an impedance Z equal 
to 5 X 10° ohms. When Z was reduced to its present value, 10* ohms, no 
error could be detected. In this case the correcting term cannot be greater than 


Com( 1 ‘ ee) = 0,015 (1 aS ~ = 0.02 yf, 


which is less than 0.05% of the capacity of the smallest cell used. 


It is for this reason that the valve V3 is inserted in such a way as to have no 
value as an amplifier but merely serve to eliminate this error. 


Av. —7.1l units = 2.3 pupf. 








3.24 Temperature variation. Temperature variations may affect the results 
in two ways. The dielectric constant of electrolytes has a large temperature 
coefficient (0.4% per degree). Also the crystal frequency is dependent to 
a small extent on temperature. For these reasons the whole apparatus was 
enclosed in a box which could be kept at constant temperature. Temperature 
variations in the whole box were less than 0.1°. A Beckmann thermometer, 
immersed in the oil bath in which the cell was also placed, was used for measur- 
ing the cell temperature. All final readings were taken at 2.00 + .01° on the 
Beckmann, corresponding approximately to 23° C. 


3.3 Errors Due to Cell Construction 

A number of errors may occur if suitable precautions are not taken in cell 
construction. These are common to the cells used in both methods, except 
that in the three-electrode cells used with the bridge no errors due to in- 
sulators present in the field can occur. 
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3.31 Inductance of cell leads. Reference to Fig. 9 shows that the true 
capacity of the cell bears the following relation to the measured capacity :— 
Lg* 


Cus 2 c(1 + 1c — +f 


where L is the inductance of leads, and g, the conductance of 
cell. Terms in L? are neglected. 


This formula for inductance error has also been derived by 
Lattey, but the cells he used were so small that the correction 
was unimportant in his work. The use of small cells, how- 
ever, greatly decreases the sensitivity. Jezewski (9) neglects 

2 
the second term, =, which may have influenced his results as 
shown below. Most other observers have neglected both 
R terms. The importance of the correction is illustrated below. 


Fic. 9. Cir- For a loop of wire 1 cm. in radius and made of wire of 
cuit equivalent to 


cell and leads for 1 mm. diameter the approximate inductance is 5X10-8 
—— cor- henries. It is very difficult to make cell leads with lower 

‘ inductance than this. For this value and a cell, of 10 pyuf 
capacity when full of water, containing N/100 sodium chloride at two mega- 


cycles’ frequency, 
C = 10(1 + (5 X 10%) — 0.01) 
2 


: ; po L 
For this case the term w*LC is obviously negligible, whereas the term = 


is 1%. The first term becomes more conspicuous at higher frequencies and 
both terms become larger for cells of larger capacity. 


Several experimental tests were made which confirmed the preser — >f errors 
due to inductance. In the case of one cell, the impedance was 1. :asured 
with the cell short-circuited and was found to be about one ohm, corresponding 
to L=10-’. The impedance was reduced by placing a 0.05 yuyuf. condenser 
in series with the leads, indicating that the impedance had been inductive. 
This agrees in magnitude with the values of ZL found below. However, no 
direct method is sufficiently sensitive to evaluate so small an inductance, 
so the indirect method described below must be used. 


3.32 Insulators in the solution. Jezewski points out that when the insula- 
tors supporting the electrodes project into the solution the apparent capacity 
of the cell is increased when the concentration is increased. The same effect 
is produced by dirt on the electrodes. A small spot of shellac put on one 
electrode was shown to produce errors as large as 50% for fairly concentrated 
solutions. In designing the cells used here Care was taken to avoid these 
errors. The check obtained with two entirely different types of cell shows 
that this and similar errors have been eliminated. 


3.4 Corrections. The correction for inductance was made as follows. Cell 
B shown in Fig. 4B was constructed with movable plates so that the capacity 
could be changed without greatly varying the inductance. If readings are 
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taken with such a cell for various capacity settings, C,, the measured capacity 
for a given solution is 

C = Cz (ec + Ac) — LG? = C, (ce + Ac) — ay*C?? 


where 7 is the conductivity of the solution and G is the conductance of the 
cell. 


Then 


C—eC. Ae 
ee 


The term on the left is the measured difference between the capacity of 
the cell when it is filled with solution and that when it is filled with water. 
If this measured percentage change is plotted in a graph against C, the curves 
should be a family of straight lines, one for each concentration. The slopes 
of the lines should be proportional to the inductance and to the square of the 
conductivity. The intercepts on the axis of ordinates give the corrected 
percentage change in dielectric constant of each solution. 


Similarly the cell C (Fig. 4C) may be used with the plates separated by 
different distances. The inductance is not greatly changed in this way 
although the capacity may be changed by a factor of two. The level of the 
solution in this cell is not as important as it is in the silver cell, so that more 
accurate readings can be made. The level of the liquid in both cells was 
kept constant by filling with a measured volume of solution. The results of 
the measurements are given in the next section. 


4. Measurements 


4.1 Absolute measurements on sodium chloride. A silver-plated brass cell 
was constructed in such a way that its capacity could be varied (Fig. 4B). 
This cell was used to determine the inductance correction. Table II gives 
the averages of three sets of readings made with this cell. The capacity, C,, 


TABLE II 
AVERAGE READINGS WITH CELL B (SILVER PLATED) 


| fe | ac | 
Concentration, AC G { c G 
millinormal ——<$—$—— : 
C,. = 740 units C,. = 545 units 








2.8 1.31 
2.4 2.61 
—2.0 5.18 
~ii.5 4° 3.4 





C. = 393 units 
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of the cell when filled with water is corrected for the capacity of the leads by 
calibrating with benzene. The units in which it is expressed are scale divisions 
of the variable condenser. 
1 puf. = 3.20 units. 

The column marked AC gives the apparent increase in capacity when water is 
replaced by thesolution. The cell conductance, G, is calculated from tabulated 
conductivities and agrees within 5% with the value measured on a commutated- 
direct-current bridge. It is used in making the inductance correction. 


All measurements were made at 23.0° C. and 210° cycles per second. 


Assuming as above, AC = = C,—LG? (Section 3.4), and C,=AG, where 


i ; a ae 
A is a constant for a given solution, — G =A— -LG. 


The data in Table I are plotted in Fig. a For each solution there is no 


. : ee C—O 
systematic departure from linearity in the curves for | against G. The 


slope of the curve gives the 


Cea B lead inductance, LZ. It is to 
——— 
i 


be noted that any error other 
than one of the form (con- 
Thre O40 stant X G), would be detected 
ae because it would produce curva- 
ture in the lines in this graph. 
There is no reason to expect 
an error of the form (constant 
xX G). The lines for each con- 
centration should be parallel 
but not necessarily coincident. 
 siienheaneaedaand The points for 1.0 and 2.0 


* G+ Cele conoucrace (musimues) ~~ Mmillinormal concentrations are 
Fic. 10. Determination of correction for cell lead the most reliable and only 


industonce. these were used. 
Using the value of L (1.35 X 10-7 henries) determined from Fig. 10, the 





corrected values for =s C, are found from the formula. Dividing by the 


capacity of the cell the following values for the percentage increase in di- 
electric constant are found. 
TABLE III 
PERCENTAGE INCREASE IN DIELECTRIC CONSTANT 
(Sodium chloride in Cell B) 





Concentration C. = 740 C, = 545 C. = 462 C. = 393 Average 
0.5 0.4 0.6 0.5 
1.0 0.9 1.0 1a 0.9 1.0 
2.0 1.8 L7 2.0 > 1.95 
3.0 2.9 2.5 3.0 2.8 
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In the same way readings were taken with Cell C for two capacities. 


TABLE IV 
READINGS WITH CELL C (PLATINUM-GLASS) 


C = 398 C = 153 


Concentration, ©. |_———_—_<______ 
millinormal G G 


, , 
micromhos micromhos 


0.6 
1.0 
1.5 
2.0 
3.0 3.15 
4.0 3.5 
5.0 3 2.4 
6.0 1.5 
8.0 —1.9 
10.0 —7.0 


The values of Z calcu- TABLE V 
lated from the data for 
each concentration which Concentration | 2 | 3 | 4] 5 Average 
was measured in bothcells = 7, henries- | 1.6] 1.4] 1.3 | 1.4] 1.42 x 1077 
are shown in Table V. 


Using this value of L the corrected readings are as shown in Table VI. 


TABLE VI 


PERCENTAGES INCREASE IN DIELECTRIC CONSTANT 
(Sodium chloride in Cell C) 


Concen- Concen- 
tration, “te i 3 Average tration, Sts Be 3 
millinormal millinormal 
2.5 
3.3 
4.8 


5.5 





2 
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These values are plotted 
in Fig. 11. 


4.2 Comparison measure- 
ments on several electrolytes. 
The errors for which cor- 
rections were made in the 
last section are functions 
of cell conductance only. 
Different electrolytic solu- 
tions of the same conduc- 
tivity can therefore be- 
compared without any 


corrections. This was done 
CONCENTRATION fMIULLIEQUIVALENTS PER LITRE with several solutions to 
Fic. 11. Dielectric constant of sodium chloride solutions. test the method. The 
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high-frequency conductivity of the solutions used does not greatly differ 
from the d-c. conductivity, so the solutions could be matched on a 
commutated d-c. bridge. The value of R, the compensating resistance 
in the plate circuit of the detector, is a measure of the conductivity. As 
a check, the readings of this resistance were recorded and showed that 
when the d-c. conductivities were made equal the high-frequency conductiv- 
ity did not differ by more than 2%. The values of AC here are the differences 


TABLE VII 











Solutions matched with 
0.002 N common salt 


Solutions matched with 
0.001 N common salt 








Substance 

he be, 

€ € 

% % 
NaCl i 2.3 
HCl 1.2 2.4 
K:SO, 1.3 2.6 
CuSO, 2.0 o.2 
MgSO, 2.0 3.3 
K,Fe(CN)s 1.5 2.9 

Solutions matched with Solutions matched with 
0.005 N common salt 0.010 N common salt 

Conc., AC, Ae 

milli- |, ¢.' % 

C. = 153 % 
NaCl -- 9.5 
HCl 0.5 10.0 
K;SO, 1.4 10.9 
CuSO, 3.0 12.5 
MgSO, 2.9 ig. 
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TABLE VIII 
SuMMARY—% INCREASE IN DIELECTRIC CONSTANT 


a NaCl K:SOx CuSO, | MgSO, | K.Fe(CN)s 


between the capacity with the solution in the cell and that with common salt 


ca Ae. s 
of the same conductivity. <a then calculated by adding values from the 


graph for sodium chloride (Fig. 11.). 


These results are plotted in Fig. 12 
with conductivity as abscissa and in 
ii Fig. 13 with concentration as abscissa. ° 
* Het 
AK, Se 
KCuS04 
0 1795% 

U ky fe Cn, 


CONCENTRATION 


Fic. 12. Dielectric constant of electrolytes Fic. 13. Dielectric constant of electrolytes 
plotted against cell conductance. plotted against concentration. 


5. Discussion of Results 


5.1 The work described here differs in several respects from previous 
work with resonance methods. The method of finding the resonance point 
as described in Section 2.32 is not new but has not been used before in work 
on electrolytes. As is well known, it is possible by this method to increase the 
sensitivity greatly, for it is possible to adiust the size of the condensers so that 
the setting is made at the steepest part of the resonance curve instead of at 
the flat peak. This method is possible only when the resonance curve is 
symmetrical. It has been shown both theoretically and experimentally that 





PERCENTAGE INCRE, 








396 CANADIAN JOURNAL OF RESEARCH 


this is so. The circuit calculations, in addition to proving the symmetry of 
the curve, show that errors due to coupling are absent, thus permitting the 
use of a driving coil very close to the resonance coil. The power required 
in the oscillator is thereby decreased, thus greatly facilitating the shielding 
of the circuit. The shielding was such that the measurements were quite 
unaffected by the movements of the observer. 

The use of a quartz oscillator and of tuned filter circuits (TJ: and 73, Fig. 5) 
has been shown to give constant frequency and pure harmonic wave form. 
The necessity for these has been recognized before but the means of attaining 
them have not been used. Jezewski, however, corrected for errors due to 
frequency changes by comparison with a subsidiary quartz oscillator. 

The use of two cells of very different construction, both of which have 
several values of capacity, made possible the evaluation of error due to in- 
ductance of cell leads, and also assured the absence of other errors due to 
cell construction. 


5.2 Accuracy. The accuracy with which readings may be repeated is better 
than + 0.1% for solutions up to 3 millinormal and +0.5% for the strongest 
solutions measured. The accuracy of the final results, however, is not as 
great since the correction for inductance of the cell leads is large and com- 
paratively difficult to estimate. From the agreement obtained with the two 
cells, it seems likely that the results up.to 3 millinormal are correct within 


+0.5%. For the strongest solutions the error may be as large as +2%. 


5.3 Comparison with results of other observers and with theory. Debye and 
Falkenhagen (2) predict an increase in dielectric constant with increasing 
concentration due to the “ionic atmosphere.’’ This effect decreases with 
increasing frequency. Tables giving the predicted values in special cases are 
given by Falkenhagen and Vernon (5). As stated in the introduction, the 
experimental results of many 
observers for very high fre- 
quencies are in agreement 
with the theory. 

For lower frequencies, 
however, the agreement 
holds only in the case of 
Jezewski’s (9) work. As 
shown in Fig. 11, the 
author’s results, which agree 
approximately with those of 
Lattey and Davies (10, 11), 
show an increase in dielec- 
tric constant much larger 
than the predicted values for 
" CONCENTRATION * sodium chloride. For other 
salts the difference from the 
theory is of the sanie order. 


w 
= 
-& 
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Fic. 14. Deduction of approximate lead inductance 
from data of Jezewski and Kameckt. 
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The apparent agreement of Jezewski with theory seems to be accidental, 
as he has not properly corrected for the inductance of cell leads. In Table 
III on page 565, reference (9), Jezewski gives results for several salts. The 
results for concentration 0.0005 N were obtained with a cell of 20 uf. cap- 
acity, while the remainder were obtained with a smaller cell of 10 puf. capacity. 
Fig. 14, which is drawn from Jezewski’s table, shows a definite difference 
between the results with the large cell and the curve for the small cell. If, 
as seems likely, this difference is due to inductance, an approximate evalu- 
ation of the error can be made by the method outlined above (Section 3.4). 
This was done on the assumption that the inductance was the same for both 
cells. As this may not be so, the correction is not to be relied on, but is given 
merely to show a possible explanation for the discrepancy between the results 
of Jezewski and those of Lattey and of the author. 


The “‘corrected” values for sodium chloride, from Jezewski’s work, are 
plotted in Fig. 11. 


No attempt is made to account for the disagreement with theory as more 


data would first be necessary. It is hoped that these measurements will 
before long be extended to higher frequencies and to a wider range of salts. 
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THE TRANSMISSION OF SOUND BY A SERIES OF EQUIDISTANT 
PARTITIONS! 


By Donatp G. Hurst* 


Abstract 


The transmission of normally incident sound waves through a series of 
similar, equally spaced, partitions has been calculated for the cases of infinite 
partitions and of circular partitions. The transmission characteristics depend 
on the solution of simultaneous linear difference equations for the velocity 
potentials. The array of partitions possesses interesting properties as a filter, 
and to illustrate these the change in intensity levels due to glass partitions is 
calculated numerically for 1, 2, 5, and 10 partitions. 


Introduction 


Problems of acoustic insulation, to which much attention is now being 
given (8), suggest the calculation of transmission through a number of 
partitions of building material of various types. The question has recently 
been made the subject of an investigation by Constable at the National 
Physical Laboratory (2). 

The theoretical treatment in this paper is concerned principally with the 
transmission of normally incident sound waves through a number of equally 
spaced, infinite, rigid partitions of the same mass per unit area. In these 
circumstances the velocity potentials in the air spaces may be determined as 
solutions of simultaneous linear difference equations arising from the bound- 
ary conditions at each partition. 

This simplified problem shows that the array of partitions has interesting 
characteristics as an acoustic filter, 7.e., the arrangement is selective, readily 
transmitting certain frequency bands of the acoustic spectrum while remain- 
ing opaque to others. A numerical case has been worked out completely for 
systems of 1, 2, 5 and 10 glass partitions of mass 0.72 gm./cm.? and spacing 
5 cm. The regions of attenuation have also been found for a semi-infinite 
sequence of cellophane partitions of mass 0.0036 gm./cm.? and spacing 5 cm. 

The general procedure of this paper might be extended to equally spaced 
partitions through each of which the propagation of sound is taken into 
account. If desired, the dissipation of energy in each plate may be taken into 
account as well. 

When the partitions are considered to be of finite size and inserted in a 
rigid wall, the transmission characteristics are affected by diffraction phen- 
omena at the first and last partition. If we consider each partition to te a 
rigid circular piston the problem is still capable of rigorous solution in terms 
of cylindrical wave functions, making use of Hankel’s inversion theorem 
in the manner employed in a recent paper by King (6). It is beyond the 
scope of the present paper to consider any except rigid partitions, so that 
no allowance is made for other possible modes of vibration of actual panels. 


1 Manuscript received December 7, 1934. 
Contribution from the Department of Physics, McGill University, Montreal, Canada. 
2 Holder of a bursary under the National Research Council of Canada. 
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When the incident wave is inclined to the first partition, the latter does 
not move if it is infinite and rigid, since the resultant pressure variation on 
it is negligible in comparison with its mass. However, for finite panels we 
should expect a wave slightly inclined to the normal to act similarly to a 
normally inclined wave. The problem may be stated for circular partitions 
capable of oscillating under specified constraints, and the solution then 
depends on the corresponding problem of oblique diffraction by a piston sur- 
rounded by an infinite, rigid flange. Such problems are extremely difficult 
and are beyond the scope of the present paper. 


Infinite, Plane, Rigid Partitions 


(a). Let us consider the case of p 
equally spaced rigid plates or partitions 
of mass m, per unit area at a distance 
b apart (Fig. 1). 

A plane wave is normally incident 
on the first. The velocity potential ¢, 
apart from the time factor e* which a) 44 2 2a 
will be omitted for purposes of brevity, Fic. 1. Series of plane, infinite, rigid partitions. 
is determined from the equation 


d@? w 
i g=-—xK¢, (1) 





where c is the velocity of sound in air, # = 2 af where f is the frequency, x = 
w/c and z is measured normally to the plates in the direction of propagation 
of the incident wave. In the m” interspace the appropriate solution of 
Equation (1) is 
Gn = see + B,e**->) , (2) 
The first term corresponds to the forward wave, the second to the reflected 
wave. 
The particle velocity is given by 
— = — (d¢/d2) , (3) 
where £ is the displacement and é as usual denotes d¢/dt. At the (m+1)" 
partition, continuity of velocity requires 
bus = — (yn/d2)eand = — (Oe—n+i/92)2anb- (4) 
The variation of pressure is given by 
5p = poy ’ (5) 
where p, is the mean density of the air. 
It may be supposed that each partition is able to perform free oscillations 


of a simple harmonic type under an elastic restoring force uf per unit area, 
so that the dynamical equation for the (n+1)** partition is 


(Spn)sand = (5Pn+i)emnd = (mot + Be)emnd- 
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For periodic oscillations, § = —w*t, so that Equation (5) gives 
Po (gn = Qn+i)smnd = Me (1 bac wo?/w*) Esa ’ 


and consequently 


. a (Gn eset Pn+1 rend er O¢n ——— Dens 
fend (mo/ po) (1 — wo?/w*) on 02 a 02 ™ (6) 





where w/2zm is the impressed frequency, and w,/2z is the natural frequency of 
the partitions. From Equations (2) and (6) is derived the expression 


ie + B, dana Anye'*> — Bayie~**> 


=1 —_ =i ixdb —ixb 
(mo/ po) (1 aoe Wo?/w*) ix(An B,) tx(Anye Baye Ds 


On writing 

y = eo sy = § (xme/po)(1 — wo?/w) (7) 
we have for the A’s and B’s the simultaneous difference equations 

An — Anny = Ba — Bayry™ (8) 
and 

A,(1 — tv) — Angiy = — ivBy. (9) 


By elimination of B,, there is obtained for A, the linear difference equation 
Any: — 2Ancosha + Ann = 0, 


where 
cosh a = cos xb — vsin xb. (10) 


It is easily seen that the solution may be written 


A, = Asinh (na +e) = A’e™™ + AXe-"2. (11) 
Similarly 
B, = Bsinh (na + e) = B’e™® + B’e-"e, ' (12) 


where A, A’, A’’, B, B’, B”’, «, and & are independent of m. That vatue of a 
which is positive, or whose real part is positive, is chosen. 


(b). When the number of partitions is infinite it is seen that the appropriate 
solutions are 


A, = Ave: 3B, = Be, 
and hence there will be attenuation of amplitude if a has a real part, a condi- 
tion which requires that 


(i) cosha > 1 or (ii) cosha < —1. 


If we define u = }xb, the first condition becomes 


sin 2u (tanu + v) < 0. e 

Thus if 
sin2u > 0, we must have tanu < — » ) (13) 
sin 2u < 0, we must have tanu > — » f. 


The second condition becomes 


sin 2u (cotu — ») < 0. 
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Thus if 


sin 2u > 0, we must have cotu < 5 (14) 
sin 2u < 0, we must have cotu > » 

In studying Equations (13) and (14) as conditions of attenuation of ampli- 
tude, it will be sufficient as a preliminary step to neglect the effect of the 
natural frequency of the partitions and to write w, = 0. This makes 

vy = § xmo/po = mou/(bp.) = 1, 
where 

n = mo/(bpo). 
Since sin 2u = 2 tanucos*u, y 
sin 2% and tan u have the 
same sign, and since nu is al- 
ways positive we obtain from 
Equations (13) and (14) 


0>tanu > — TH Ce = 0) (15) 






0 <cotu < nu 


as the possible conditions. 
These conditions are illus- 
trated graphically in Figs. 2 
and 3. It will be noticed 
that the transmission bands 
become successively narrower 


with increasing frequency. Fic. 2. Tan u plotted against u; (i) isy = — ay 


The effect of the term ( = 2) . The regions of attenuation are marked by the 


i- wo'/ w is shown by the heavy lines on the axis of u. They are contained between 
Curves (ii) of both figures which 


represent y = Fv respective- 

ly. For w > -w, there isa y 
small reduction in the width 
of the attenuation bands. . 

However, for w < w, the aa 
effect may be considerable, 
and as is seen in Fig. 2, a 
band of attenuation corre- 
sponding to the condition 
0 < tan u < —vis introduced 
near the origin. Inasmuch, 
however, as the natural fre- 
quency of an actual panel is 
usually low, relative to audi- 
ble frequencies (if we neglect ; 
possible radial and transverse Fic. 3. Cot u plotied against u; (i) is y =  . Here 
modes of vibration), the dis- the regions of attenuation are bounded by the roots of 
turbance due to w, does not coy = =. and u = A a: The combined system 


extend appreciably into the of attenuation bands is shown at the bottom of the figure. 
audible frequencies. 


the roots of tanu = —Zuand u = x, 2x, 3r. ‘ 


ee el 
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Transmission Characteristics of a Finite Number of Partitions 
The existence of alternate regions of attenuation and perfect transmission 


in the semi-infinite case leads one to expect wide variations in the character- 
istics of a finite number, p, of partitions. 


The sound wave generated by the last partition is a progressive wave 
traveling in the direction of the incident wave, and it will’ be supposed that it 
is not reflected, so that, referring to Equation (2), B, = 0. According to 
Equation (12) therefore 

B, = — Bsinh (p — n)a. (16) 
On giving m the value n = p—1 in the difference equations (8) and (9), and 
remembering that B, = 0, 


Ay-1 — YAp = Bp = — (tv) {Apia (1 — iv) — Ap} 


is derived. 
Let B = pa + «. Then from the above equation 
cosh 8 = cos xb + (1/») sin xb (17) 
and 
sinh 8 = (i/y) sinh a. (18) 
Finally 
A, _ sinh (pa +4) _ sinh 8 (19) 
A. sinh « sinh (8 — pa) 


The intensity level of the incident sound (1) is 
a = 10 logio (I,/10-*) decibels, 
where J,/10-"* is the ratio of the incident sound intensity in watts/cm? to the 


standard reference intensity (10—* watt/cm?). Similarly the intensity level 
of the emerging sound is 


a» = 10 logy (Ip/10-*) db. 
The difference in levels due to the array of partitions will be denoted by 
dpe, and is given by 
Spp = a — a, = 10 logio (I./I>). 
Since the waves are plane, 
(I./Ip) = | Ae/A,l*, 


so that 
bpp = 20 log | Ao/A» | 
sinh (8 — pa) 
= 20 logio ss . (20) 


It may be shown, as indeed is to be expected from energy considerations, 
that 


|} At? =| BI +] Ap PP 


1.¢., the intensity of the incident wave is equal to the sum of the intensities 
of the reflected wave and the transmitted wave. 
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Single and double parti- 
tions, corresponding to p = 1 
and p = 2, have been exam- 
ined theoretically and experi- 
mentally in recent papers by 
Davis (3) and Constable (2). 
The formulas obtained, trans- 
formed to the present nota- 
tion are 

| A/4i kt? =1+ 
(single panel) 
and 
| Ao/As |? = 1 + 40° cosh? a 
(two equal panels) 

to which Equation (19) re- 
duces for the values p = 1 
and p = 2. 


Numerical Calculation 


To illustrate the results of 
the preceding sections some 
numerical examples may be 
considered. We take as par- 
titions panes of glass for 
which m, = 0.72 gm./cm.?, 
and which are separated by 
a distance b = 5 cm. Ata 
temperature of 20° C. and 
normal pressure, pp = 1.20 
10-* gm./cm’, c = 3.44 X 
104 cm./sec. 


Hence 7 = m,/(bp.) = 120. 
The natural frequency is 
taken as 18.3 cycles per 
second, thus making » = 
300x — (300x)—' and greatly 
simplifying the arithmetic. 

The method of computation 
is to choose x, then compute 
v, cos xb and sin xb, and finally 
cosh a, cosh 8. From these 
a and 8 are readily found. 
With the aid of Kennelly’s 
chart (5), it is a simple matter 
to obtain | sinh 8 | and 
| sinh( B— pa) |. One calcula- 
tion of a and @ suffices for any 
number of partitions. dps, is 
calculated from Equation (20). 
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Fic. 4. Reductions in intensity levels (measured in 
decibels) due to various numbers of partitions, plotted 
against a logarithmic scale of x (=2xf/c). The standard 
partition is taken as glass, having a mass per unit area 
mo = 0.72 gm./zm2 and natural frequency f. = 18.3 
cycles per second. The spacing between partitions is 5 cm. 

The separate curves are: A, one partition; »/ B, two par- 
titions; C, five partitions; D, one partition having m, = 5 
X 0.72 gm./cm.2 and f, = 18.3 cycles per second; +.e., 
its weight equals the combined weight of the five partitions 
in C; E, ten partitions. The transition to the sem1-infinite 
case can by seen taking place as the number of partitions 
is increased. In all curves the first zero is at the resonant 


frequency of a single partition; the other zeros (p — 1 in 


number in the first transmission region ) occur at the points 
where pa = im X (integer) and 8 # O, since then 
| A./Ap| = J and dnp = 0. 





EOE eEore— 
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The computations were performed for 1, 2, 5, and 10 partitions and for a 
single partition of mass per unit area equal to that of the five partitions com- 
bined. The results are plotted in Fig. 4. 


Although the higher frequency ends of the attenuation bands are deter- 
mined solely by the spacing, the width of the bands decreases with decreasing 
weight of the partition. The following table gives a comparison for the glass 

TABLE I partitions considered above, 
REGIONS OF ATTENUATION and similar partitions (¢.g., 
of cellophane) of mass per 
unit area m, = 3.6 X 107 
gm./cm.? so that 7 = 0.60. 
Using Equation (15) it is 

2230 3440 
4850 : f 2 6880 found (4, pp. 36 and 37) that 
7840 < f < 10,320 the regions of attenuation are 

7 A 

pee ne ee oe given by the data in Table I. 
Thus, while with glass the 
regions of transmission are 
extremely narrow, they are 





Glass 






Cellophane 





200 <f < 3440 
3446 <f 
Further transmission regions 
are too narrow to be indi- 
cated and all frequencies 
above 200 can be considered 
as attenuated. 





TABLE Il much wider with cellophane 

and might possibly be ob- 

or med Spe (f = 3300 cycles per second) = served experimentally. Ow- 
partitions Glass | Cellophane ing to the smaller values of 
vy for cellophane it is to be 

1 45 2.6 expected that the reduction 

; Too ae a a will not be as great as for the 

10 practical 44 heavier partitions. As an 


significance example, the reductions for 
glass and cellophane at 3300 
cycles per second are compared (Table II). 
It should however be recalled that 10 partitions of the cellophane have a 
mass only 0.05 times that of a single pane of glass. 


Transmission Through a Series of Circular Panels Set Into a Rigid 
Wall, Taking into Account Diffraction Effects 


The cross section of the series of panels 
is shown in Fig. 5. 
In terms of cylindrical co-ordinates (p, 2) 
the velocity potential may be written (6) 
in the form 

© us 
bo = 2Agcos ut +, [face (*) Jno) aa, 
(21) 

where p2? = 2 — 2, and 2 is measured 
from the first plate in the direction of 





GEL UY 


Fic. 5. Cross section of series of circular = propagation of the incident wave. The 
first term on the right represents a plane 


panels se into a rigid wall. 
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wave which undergoes perfect reflection at the wall, and the second term 
is a solution of the wave equation representing reflection by the circular 
panel. The function f,(d) is determined by the boundary condition 


8 Po &, O<p <a 
-G )_- i a@a<p<o (22) 


and £, is constant since the partitions are rigid. 
Ser) =f? fonrrooian. 


According to Hankel’s inversion theorem, if 

Flo) = f° f0)rJa(r0)dr 
then 

fo) = [° Flo)edJo(rodo . 
In the present instance 

fin) = [°(— &) 0 Jodo) do = — & ($) 102), (23) 
so that 

= thane = th [+ Je(Xp) Ju(da) 2. (24) 
The dynamical equation for the first partition is 

[200 - SPdee Axode = atm, (1 — w/t) be. 


Equations (5) and (24) give 
— (bo) f° ode f° 7.000) 1000) 2 - 4 (dee = O/ 0b. 
° ° 7 


On performing the integration with respect to p and writing 
M, = ii, / * (na) 2 (25) 
° Au 


we finally obtain 


a 2Ao — (1)s—0 — 2Ae Ae — Meiers - 
> “704 6 -ae ~ -(% e) (26) 


where 
VW = 2 (vo/x + M, — 1M) (27) 
The integrals M, and M, have been evaluated by King (6), who finds 
M, = Hy (2«a)/(2a) ; M2 = {1 — Ji (2a)/(xa)} /2n. 


In the notation introduced by McLachlan (9), 
Gi = 1 — Ji(2na)/(xa) ; Ge = Hi(2«a)/(«a) 
so that 
My = Ga/(2x) ; Me = Gi/2k. 


The functions G, and G2 are plotted on page 54 of McLachlan’s treatise. 


It is seen that the solutions (Equation (2) ) derived for infinite partitions 
satisfy the boundary conditions in the interspaces, so that we can write for 
the velocity potential y, in the ”** interspace 

gn = Anemixe-md) + Bye nh, (m9 = 1,2,..,p—l). , (28) 
It is noticed further that Equations (8), (9), (10), (11), and (12) also apply. 





a 
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The velocity potential for the wave generated by the last panel may be 
represented by 
& r 
ve = [Pe T0r0) 2 f0)0d, (29) 


the origin being at the centre of the last panel and f,() being such that 


=) {fi 0<p <a 
08 Jeno a 6@<29 <<, 


As before we obtain from this boundary condition 
er = ok, fe J.09) La) &, (30) 
° u 


and the dynamical equation for the last partition gives 
Vin = — WY (8971/02) = (Ye) 0 - (31) 
From Equation (9) it follows that 
{(1 — t») A cosh e, — yA cosh a + «& + i» B cosh e2 }sinh na 
+ {(1 — iv) A sinhe — yA sinha + « + i» B sinh eg }cosh na = 0, 
and the requirement that this must hold for all values of m leads to 
Be=Ayj;a=a-— 8, (32) 
where 8 is defined by Equation (17). 
Returning to Equation (31) and writing 


f=p-larta, (33) 

it is found that 
sinh 8 

taah f = See Tl — ac) + ee (34) 
so that 

1—ixv | sinh (8 — £) 

ftw. “~aT (35) 
On making use of Equations (34) and (35), we may write Equation (26) in 
the form 

As _ (1 + ix¥) 7 sinh B sinh (2 ¢ — B—p— 2 a) (36) 

A 2 sinh ¢ , 


and it follows from Equation (6) that 
és 


(Ep-s)ongmie = — (0 ¢-1/02)s=5-0 37 
ix (Ap. — Bp) = 2icA(1 + ia®) sinh ¢. (37) 


As shown by King, the radiation output is given by 


(dW/dt} = a*puo| ky |*Mz, (38) 
and the rate of transmission of energy in the incident wave across the area 
xa? is 

[dW./dt] = wa*hp,| Az |?xw . ¢ (39) 
Hence a “measure of reduction”’ for transmission through the series of panels 
is given by 
[¢W./dt} . £ | Af ms 1 (1 + ix)? sinh 8 sinh (2 — B — p — 2 a) [ (40) 
[dW/dt}~ 2Mz\ & 32 Max sinh *f¢ ; 








As before Kennelly’s charts aid in carrying out numerical calculation. 


TRANSMISSION OF SOUND BY A SERIES OF PARTITIONS 407 


Summary and Conclusion 


1. In a semi-infinite sequence of equally spaced, infinitely extended par- 
titions having the same mass per unit area, there are alternate ranges of the 
frequency corresponding to complete transmission and complete attenuation. 
The upper frequency ends of the ranges of attenuation are determined by 
spacing alone, viz., by the equation: 

Integral number of wave-lengths = 2 X spacing; 


but the widths of the regions decrease with decreasing mass per unit area of 
the partitions. The regions of attenuation are regions of perfect reflection. 


2. In a finite sequence, the transmission is alternately greater and less 
than that of a single partition of the same total weight, the alternate regions 
of frequency corresponding to those of the limiting case of a semi-infinite 
sequence. The effects of viscosity and mechanical damping would probably 
bring the curve nearer to that for the single partition. 


3. The equations obtained for circular panels indicate a behavior somewhat 
similar to that for infinite partitions although no numerical calculations were 
made in this case. 
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REVIEWS AND NOTES 


Note on the §- and 7-Rays of Radium D 


It has been found in this laboratory that of every 100 atoms of Radium D 
disintegrating, approximately 3.5 emit a primary y-ray and 25,Z radiation, pre- 
sumably characteristic of atomic number 83. Hevesy and Lay (1) have shown 
that for this atomic number, the fluorescent yield of L radiation is about 40% 
and of M radiation 24%. If it is assumed that this is true for the internal 
absorption of y-rays, it is found that the total number of §-rays ejected from 
L levels is 25/0.40, or 62% As the number of §-rays ejected from the M 
and N levels is approximately 60% of that from L levels, the number of 
B-rays ejected from M and N levels will be 0.6062, or 37%. As 62+37+ 
3.5=102.5, there is strong evidence that every atom of Radium D on dis- 
integrating emits a y-ray, a point of some importance. The internal con- 
version coefficient of the y-rays is 0.62 for the L levels and 0.37 for the M 
and N levels, coefficients that are much higher than those usually accepted. 
The number of tertiary f-rays will be 62—25, or 37, with energies corre- 
sponding to differences between the L and M and L and N levels, at least 
0.76X62 and 0.76X37, or 75, with energies corresponding to differences 
between M and N and lower levels. Some of these tertiary rays will have 
such small energies that they will be very difficult to detect. However, the 
total number of §-rays emitted by Radium D in 100 disintegrations should 
be 100 primary, 96 secondary and of the order of 100 tertiary. It is not easy 
to find experimentally the number of atoms emitting M characteristic radiation 
as these are so soft. The writer has observed the presence of this type of 
radiation but has not been able to form a reliable estimate of the number 
of atoms emitting it. 


1. Hevesy, G. and Lay, H. Nature, 134 : 98-99. 1934. 
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Note on the Scattering of X-rays at Small Angles 


Three methods have been used in this laboratory for investigating the scatter- 
ing of X-rays at small angles, less than 5° of arc. In the first, the Soller slit 
method, the primary beam is defined by a set of parallel slits formed by thin 
metal sheets and analyzed by a similar set of slits which are attached to a 
graduated circle, so that the intensity distribution of the scattered radiation 
may be determined. In the second, the primary beam is defined by two very 
narrow slits and analyzed by a third slit which can. be moved perpendicularly 
to the direction of the beam. By this method it can be ascertained in the 
cases of intense scattering whether the beam is broadened or whether there 
is a scattering of the rays out of the beam. This is called the movable slit 
method. In the third method, the same three slits are used, but are first 
put in line and then the third is widened. As this is done, more and more 








SCATTERING OF X-RAYS AT SMALL ANGLES 409 


scattered radiation enters the ionization chamber placed behind the third 
slit. This is called the absorption method. The second and third methods 
are suitable only when the scattering at small angles is very intense indeed. 

The Soller slit method was introduced here by H. M. Cave, and used by 
Gray and Cave (2) to show that at 2° of arc the intensity of the scattered 
radiation from gas molecules was proportional to the square of the number 
of electrons in the molecule. It was also used by W. H. Zinn in an examination 
of the radiation scattered by water and aqueous solutions. He showed that 
with the more dilute solutions the solute behaved as a gas, the effective mole- 
cule being approximately the solute molecule minus the displaced water. 
Miss M. F. McDonald also observed the very important result that artificial 
graphite scattered X-rays intensely at small angles (it was expected that the 
scattering would diminish to zero as zero angle was approached). Gray 
and Zinn (3) continued this work and showed that charcoals in particular 
displayed this effect. Silicon powder obtained by forming a spark under 
water between silicon electrodes also showed it. 

It appeared that this scattering was somewhat similar to that shown by 
gases, and later work by W. E. Bennett has shown this to be the case. Using 
the second and third methods, Bennett has been able to show that there are 
two types of scattering shown by powders, one similar to that referred to above, 
and the other shown by microscopic crystalline powders where the scattering 
is really due to refraction in, and reflection from, the surface of the particles 
forming the powder. A theory of the refraction in and reflection from, 
spherical particles has been worked out by von Nardroff (4). He has deve- 
loped a theory following the method used by Debye (1) to account for the 
charcoal type of scattering, which explains the facts fairly well. The charcoal 
is assumed to consist of sub-microscopic particles of diameter of the order 
10-*to10-’cm. The charcoal is then treated asa gas. One of the interesting 
things that has been found is that the larger particles cannot be entirely solid, 
for when the charcoal is wetted by gasoline, the scattering from them is prac- 
tically zero. It is assumed that they consist of chains of smaller particles. 
The total intensity of the scattered radiation should be proportional to the 
square of the wave-'ength of the primary rays, and Bennett has shown this 
to be true. Some idea of the intensity of this scattering may be obtained 
from the fact that with a primary wave-length of 2.0 A, one of the specimens 
of charcoal used gave a mass scattering coefficient of 30 for angles of scattering 
between 3’ and 3°, the photoelectric absorption coefficient being 12.0. 
Artificial graphite, of course, is not a true crystal, but presumably a mixture 
of graphite crystals and charcoal. 
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